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Preface 


Conventional computing or hard computing often depends on a precisely stated 
analytical model and many times involves a lot of complexity during computation 
along with high computation time. Nowadays, there are a lot of real-world engi- 
neering problems, which cannot be studied precisely due to the presence of im- 
preciseness and uncertainties. To tackle these types of problems, soft computing 
techniques have found wide applications and have been proved to be a powerful 
problem-solving methodology because of their strong learning and cognitive abil- 
ity and good tolerance of uncertainty and imprecision. Basically, they are an ex- 
tension of heuristics and help in solving typical problems that are hard to model 
mathematically. Moreover, they are easy to accommodate with changed scenario 
and can be executed with parallel computing. Needless to say, soft computing 
techniques are an emerging approach, which includes techniques such as fuzzy 
logic, evolutionary computing, artificial neural network, and applied statistics. 
An interesting fact to be noted down is that they are actually distinct applied tech- 
niques to solve a wide range of problems but when applied in collaboration help a lot 
in solving complicated problems easily or relatively easily. They have wide applica- 
tions in fields such as machine intelligence, computer vision, VLSI design, medical 
diagnosis, pattern recognition, network optimization, and weather forecasting. 

Keeping aforementioned facts in view, we have tried to edit a book comprising 
some of the contemporary researches in different real-world problems where soft 
computing techniques have been applied. The book consists of ten chapters in 
which the following studies are conducted: 

— Chapter 1 focuses on the application of neural network models for freezing of 
gait detection and prediction in Parkinson's disease in which authors develop 
a user-independent freezing of gait detection and prediction model that would 
go along with nonpharmacological treatments. 

— Chapter 2 examines a new fuzzy de novo programming approach for optimal 
system design, which proposes a new approach for solutions of de novo pro- 
gramming problems in fuzzy environment. This approach is built on the ap- 
proach by Li and Lee (1993), which uses positive and negative ideal solutions. 

— Chapter 3 discusses a probabilistic bilevel programming in Stackelberg game 
under fuzzy environment in which researchers develop a computational algo- 
rithm to solve a stochastic bilevel programming in Stackelberg game using fuzzy 
optimization technique. 

— Chapter 4 studies intuitionistic fuzzy trigonometric distance and similarity mea- 
sure and their properties. In this chapter, the concept of intuitionistic fuzzy sets 
is introduced. Similarity and distance measures between intuitionistic fuzzy sets 
are also explained and extended these measures to intuitionistic fuzzy sets. 
Some new trigonometric similarity and distance measures of intuitionistic fuzzy 
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sets are studied and it is shown that an intuitionistic fuzzy distance measure sat- 
isfies the required identities of intuitionistic fuzzy similarity measures. 

Chapter 5 develops a mathematical model through transmutation of activation 
energy function. This work incorporates the different perspective of modeling 
of activation energies through ingression of an additional parameter of distri- 
bution function to increase the flexibility and controlling ability of modeling by 
the linear mixing. 

Chapter 6 addresses forecasting of air quality parameters using soft computing 
techniques. In this chapter, authors studied the application of soft computing 
techniques of artificial neural networks and genetic programming for forecast- 
ing of air quality parameters a few time steps in advance. 

Chapter 7 examines arithmetic operations on generalized semielliptic intuitionis- 
tic fuzzy numbers (IFNs) and its application in multicriteria decision making. In 
this chapter, authors discuss the generalized semielliptic IFN and performed its 
arithmetic operations with the help of (a, 8) cut method. In this study, they com- 
pared the proposed IFN with normal triangular, trapezoidal, semielliptic IFNs. 
Chapter 8 discusses the method for solving intuitionistic fuzzy assignment 
problem. In this research, authors use a newly proposed centroid concept rank- 
ing method for intuitionistic fuzzy numbers and solve assignment problem 
where assignment costs are taken as triangular intuitionistic fuzzy numbers. 
Chapter 9 focuses on optimization of electric discharge machining process through 
evolutionary computing and fuzzy multicriteria decision-making techniques. This 
study intended to explore the best possible set of process parameters, namely, cur- 
rent, voltage, and pulse on time during electric discharge machining process in 
order to determine the changes in performance characteristics like material re- 
moval rate and electrode wear rate for machining of high carbon high chromium 
die steel workpiece by applying the titanium nitride-coated copper electrode. 
Chapter 10 studies the fuzzy reliability of systems using different types of level 
(A, 1) interval-valued fuzzy numbers. In this chapter, authors obtain the fuzzy 
reliability of series, parallel, parallel-series, and serie parallel systems as- 
suming that all components of the system follow different types of level (A, 1) 
interval-valued fuzzy numbers (both triangular and trapezoidal). 


The book will be useful to scientists, engineers, managers, senior and postgraduate 
students, as well as research scholars. 
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Bharatendra Rai, Amruta Meshram 

1 Application of neural network to detect 
freezing of gait in patients with 
Parkinson’s disease 


Abstract: Freezing of gait (FOG) consistently reoccurs in the later phases of a pa- 
tient suffering from Parkinson’s disease (PD). Although it is treated with pharma- 
cological treatment, the impact of the medication fades with increasing duration 
of the disease and thus diminishing the mobility of a patient. This chapter aims at 
developing a neural network-based classification model that helps to detect FOG 
episodes in a patient at early stages so that lethal mishaps can be avoided. In this 
application example, we build user-independent FOG recognition system that 
would work along in conjunction with nonpharmacological medications. The 
structured system of developing a neural network-based classification model can 
be organized into three different stages. The process starts with extraction of suit- 
able features from the dataset. In the subsequent stage, patients are additionally 
grouped into two clusters depending on the FOG episodes. In the final stage, two 
neural network models are developed using feedforward network on the two clus- 
ters that were formed. The accuracy of the model is computed using sensitivity 
and specificity. 


Keywords: healthcare, Parkinson’s disease, freezing of gait, feature extraction, clus- 
tering, neural network 


1.1 Introduction 


Parkinson’s disease (PD), caused by Parkinsonism, was first discovered by James 
Parkinson in 1817. This symptom is observed due to a decrease in the neuromelanin 
neuron called dopamine causing unusual activity in the cerebrum [1]. Main features 
of PD are frequently categorized into four groups: trembling, rigidity, akinesia, and 
posture imbalance [2]. The Hoehn-Yahr score (H&Y) is used to measure different 
stages of PD in a patient, which are classified into five distinct stages indicating rel- 
ative disability levels as follows [3, 4]: 

Stage 1: Affects only one side of the body with minimal or no functional impairment. 

Stage 2: Affects both sides of the body without the loss of balance. 
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Stage3: Impairment is found; functionality is restricted. In this stage, the patient is 
physically capable of doing work independently, and the disability is mild 
or moderate. 

Stage 4: A higher level of disability, but the patient can still walk or stand without 
being assisted with a conspicuous concentration difficulty. 

Stage 5: Bounded to a bed or a wheelchair, but needs an assistant. 


Freezing of gait (FOG) is the short, inconstant loss or an explicit deterioration in the 
movement of the feet, despite inclination or willingness to walk [5, 6]. Some com- 
mon symptoms of a FOG episode involve difficulty in executing efficient stepping 
even when the patient makes an effort to do so, a feeling of difficulty to move legs, 
and also shivering of legs [6]. About 70% of FOG is observed in a patient during 
advanced stages of PD, but 2696 of FOG is reported in the beginning stages for those 
patients who have not undergone a levodopa therapy yet [5]. FOG experienced by a 
patient is mostly a brief episode that can sometimes last more than 1 min [7]. 

As the mechanism of the freezing experienced by a patient is multifactorial and 
not well understood, the treatment to find a solution is very challenging [8]. Levodopa 
is the most efficient drug used for PD [9]. Some other drugs that are available are dopa- 
mine agonists, catechol-o-methyl-transferase inhibitors, and nondopaminergic agents, 
but the effect of these drugs does not stay for a longer duration in patients [9]. 

Although the complete cure of PD is not available at this time, if FOG events 
can be predicted even a few seconds before they occur for a patient, nonpharma- 
ceutical treatments can be very helpful to avoid fatal accidents. These nonphar- 
maceutical treatments may take the form of visual, auditory, and tactile cueing 
provided to the patient [10]. 

The objective of this chapter is to build a user-independent classification model 
by application of neural network techniques to enhance FOG identification in PD. 
To accomplish this objective, the data related to patients are partitioned into two 
clusters depending on the features and then we develop a feedforward neural net- 
work model for each of the two groups. 


1.2 Data and approach 


In this chapter, we will use Daphnet FOG dataset which is freely available [4]. This 
dataset is based on experiments that were performed at the Laboratory for Gait and 
Neurodynamics at the Department of Neurology of the Tel Aviv Sourasky Medical 
Centers (TASMC). The dataset consists of FOG events captured with the help of 
wearable accelerometer sensors. 

Figure 1.1 provides a representation of the sensor data that will be used in this 
chapter. In this study, three sensors are utilized to quantify acceleration in three 
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Figure 1.1: Sensor usage flowchart. 


dimensions: first appended at the shank, second at the thigh, and third at the 
lower back. Every sensor collects reading from movements related to three axes: 
X-axis, Y-axis, and Z-axis. The horizontal forward direction reading is collected 
using X-axis, the vertical direction reading is collected using Y-axis, and the hori- 
zontal lateral information is collected using Z-axis. The dependent variable de- 
scribes the presence or absence of the FOG. In total, the database consists of ten 
variables: one dependent and nine independent variables. 

Figure 1.2 provides a flowchart depicting the data collection phase. From the flow- 
chart we can observe that in all there are ten patients: three female patients and seven 
male patients. These ten patients have the disease span varying from 2 to 30 years. 
The age for ten patients ranges from 59 to 75 years. The patients are further classified 
as belonging to the *ON" phase of medication and *OFF" phase of medication. At the 
time of collecting data, each patient completes diverse walking patterns that involve 
an arbitrary walk including 360? turn, walking in straight line including 180? turn and 
walking while performing simulation tasks. All these events were recorded using a dig- 
ital video camera. A physiotherapist is used for noting events during the experiments 
to determine start time, duration, and end time for the FOG event. In addition, the cur- 
rent activity of the patients is also recorded. 

The data collection period involved recording for 8 h and 20 min with an inter- 
val of 15-20 ms. Out of the ten patients, eight patients experienced FOG episodes. 
Overall, 237 FOG episodes were detected during the experiment. Patient 5 encoun- 
tered as many as 66 FOG episodes, whereas patient 6 experienced 10 FOG episodes 
which were least among all the patients with FOG. Figure 1.3 represents a boxplot 
of the duration of FOG events for each patient. Data of patient 4 and patient 10 are 
excluded from the plot as they did not experience any FOG event. Although patient 
6 (disease duration: 22 years) experienced only 10 FOG events, we can observe from 
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Figure 1.2: The patient data collection process. 
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Figure 1.3: Boxplot for duration of FOG events in each patient (Patients 4 and 10 are excluded as 
they did not experience FOG events). 
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the boxplot that a period of FOG events was the highest variability (sd = 14.32 s) for 
this patient. Patient 8 (disease duration: 18 years) also indicates relatively high vari- 
ability in length of FOG events. 

On the other hand, patient-7 (disease duration two years) experienced the least 
variability (sd = 2.73 s) among all patients. Patient-8 was also observed to have the 
highest average FOG event duration (mean - 14.34 s), and patient 7 had the least 
average (mean - 3.38 s). It is to be noted that patient 8 was on the *ON" stage and 
was the only patient with the highest H&Y score of 4. 


1.3 Methodology 


This section provides a detailed overview of the three stages of developing a classi- 

fication model using neural network. 

Stage 1: (Extracting essential feature for model building): In this analysis, use of a 
5096 overlapping window is utilized to derive 11 statistical features. The sta- 
tistical features are derived from the nine variables and consist of mean: 
mean1 (lower 10% of the values were trimmed), mean2 (lower 20% of the 
values were trimmed), minimum, maximum, median, standard deviation, 
kurtosis, variance, skewness, and mode. Thus, in all we develop 99 features 
from the data. 

Stage2: (Clustering analysis for model building): In this phase, data with FOG epi- 
sodes are considered for the analysis and the patients are grouped into 
clusters depending on their FOG episodes. FOG events using all 99 features 
were utilized to perform k-means cluster analysis. This analysis resulted in 
the formation of two distinct clusters: patients 1, 2, 3, and 5 are classified 
as cluster 1, while patients 6, 7, 8, and 9 are grouped as cluster 2. 

Stage 3: (Model building): In this stage, we build a neural network model to clas- 
sify and predict FOG events. In this dataset, the overall percentage of FOG 
events is 9.796, which means the percentage of no-FOG event is 90.28%. If 
we consider the worst scenario where without developing any classifica- 
tion model we make a prediction that all events belong to no-FOG cate- 
gory, then we will be right 90.2896 of the time. We can use this number as 
a benchmark and develop a classification model that has a better perfor- 
mance. In this study, we will make use of sensitivity and specificity to as- 
sess the performance of the classification model. Sensitivity is defined as 
the ratio of FOG episodes that are correctly classified out of all the FOG 
episodes observed in the data. Whereas specificity is defined as the ratio 
of no-FOG episodes that are correctly classified out of all the no-FOG epi- 
sodes observed in the data. 
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To perform classification of the two clusters, feedforward networks (nntraintool) 
in MATLAB are utilized to classify inputs according to the target classes. While 
developing a classification model on the two clusters formed, 70% of the data are 
used as training data, 15% of the data are used as validation data, and the re- 
maining 15% of the data are used as test data. The representation of cross- 
entropy versus performance using training, validation, and test data for cluster 1 
is depicted in Figure 1.4. 


Best validation performance is 0.14834 at epoch 58 
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Figure 1.4: Cross-entropy versus performance using training, validation, and test data for cluster 1. 


For developing a good classification model we minimize cross-entropy. As shown in 
the plot, the cross-entropy reduces with the increasing epochs for training, validation, 
and test data. The dotted green line indicates the best epoch, which in this case is 58 
epochs. The performance of the network is measured using four lines that capture 
train, validation, test, and best situations. If any of the three lines (training, validation, 
and testing) meet or pass near the best (dotted) line, it symbolizes convergence. 

To assess the classification performance of the mode, we develop a confusion 
matrix. Figure 1.5 provides the confusion matrix for cluster 1. 

From the confusion matrix, we can observe that there is consistency in results for 
training, validation, and test data. For training data, the model correctly classified 
8,855 events as no-FOG and 457 events as FOG. The sensitivity and specificity values 
for the training data are 66.496 and 96.196, respectively. The sensitivity and specificity 
values for the validation data are 63.1% and 96.3%, respectively. Similarly, sensitivity 
and specificity values for test data are 64.496 and 9796, respectively. This consistency 
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Figure 1.5: The confusion matrix for cluster 1. 


in results across training, validation, and test data indicates the absence of over- 
fitting problem. Another observation that we can make from this confusion matrix is 
about the events that are misclassified. It is observed that the events that belong to 
FOG category are misclassified as no-FOG at a slightly higher rate compared to a situ- 
ation where events belonging to no-FOG category are misclassified as FOG. This pat- 
tern occurs consistently for the training, validation, and also the test data. 

Figure 1.6 represents the receiver operating characteristics (ROC) curve of the 
feedforward network at epoch 58 for cluster 1. 

From the ROC we observe that the classification performance for the training, 
validation, and test data are again consistently high. The diagonal line in the plot 
depicts the classification performance with 50% accuracy. And closer the ROC 
curve to the upper left corner of the plot, better the classification performance. 

The representation of cross-entropy versus performance using training, valida- 
tion, and test data for cluster 2 is depicted in Figure 1.7. 

As observed, the cross-entropy reduces with the increasing epochs in training, 
validation, and test data. The dotted green line indicates the best epoch, which is 
75 epochs. 
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Figure 1.6: The ROC curve for a feedforward network at 58 splits for cluster 1. 


Figure 1.8 provides the confusion matrix for classification model with cluster 2. 

From the confusion matrix we can observe that the sensitivity and specificity 
values for the training data are 63.7% and 91%, respectively. The sensitivity and 
specificity values for the validation data are 71.6% and 91.9%, respectively. 
Similarly, we observe that the sensitivity and specificity for test data are 74.6% and 
90%, respectively. Once again we observe high consistency in classification perfor- 
mance for training, validation, and test data. We also note that the classification 
performance for cluster 2 is better compared to results obtained with cluster 1. In 
addition, it is also observed that misclassification of FOG event as non-FOG event is 
significantly higher compared to situations where non-FOG event is misclassified as 
FOG event. 

Figure 1.9 represents the ROC curve of the feedforward network at epoch 75 for 
cluster 2. 
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Best validation performance is 0.21957 at epoch 75 
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Figure 1.7: Cross-entropy versus performance using training, validation, and test data for cluster 2. 


The ROC curves based on training, validation, and test data show consistently 
high performance compared to the diagonal line representing 5096 accuracy. 

In this chapter, we classified all patients into two clusters to enable us to de- 
velop patient-independent model. The classification model where both training and 
testing data belonging to the same patient is a user-dependent model, whereas the 
model where training is performed on a set of patients and testing is performed on 
another set of patients is a user-independent model [11]. When classification models 
are made for an individual patient, this model is likely to have better performance. 
However, such a model is of limited use to another patient. Developing a user- 
independent model is cheaper and efficient compared to a user-dependent model 
but accuracy of both models is less. Reasons for such lower accuracy are inconsis- 
tency in data due to different walking styles, small size of the dataset, and a limited 
number of patients. 


1.4 Conclusion 


A user-independent methodology for detection and prediction of FOG events using 
combination of clustering and neural network model is illustrated in this chapter. 
This chapter is divided into three stages: feature extraction, clustering, and model 
development/evaluation. The dataset contains data from three different sensors 
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Figure 1.8: The confusion matrix for cluster 2. 


that are used to derive 99 variables with the help of 11 statistical numeric summary 
values. An overlapping window of 50% is utilized to extract the features that were fur- 
ther used to group patients into two clusters using k-means clustering. Feedforward 
neural network classification model is built for each of the two clusters. 

In this study, data from all sensors were utilized to perform the experiment. We 
used neural network-based classification model in this chapter. Further research 
may involve use of various other machine learning and deep learning methods to im- 
prove classification performance. In addition, this dataset used has a significantly 
lower number of FOG events (about 10%) compared to no-FOG events (about 90%). 
Such a severe imbalance in the two categories is likely to cause bias in results where 
no-FOG events are classified with a much higher accuracy compared to accuracy in 
classifying FOG events. This bias is observed in the classification model that was de- 
veloped for cluster 1 and cluster 2. In both classification models, accuracy in correctly 
classifying no-FOG events is higher than 90%. Whereas accuracy in correctly classify- 
ing FOG events for classification models for both cluster 1 and cluster 2 has a percent- 
age lower than 75%. This disparity in accuracy for two classes observed due to class 
imbalance is an interesting problem and needs further investigation. 
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Figure 1.9: The ROC curve for a feedforward network at 75 splits for cluster 2. 
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Nurullah Umarusman 
2 Anew fuzzy de novo programming 
approach for optimal system design 


Abstract: De novo programming, which evaluates the concept of optimality from a dif- 
ferent point of view, makes it possible to rearrange the number of resources for con- 
straint functions depending on the budget constraint. While the resource amounts of 
the constraints are used at full capacity, this arrangement improves the performance 
level of the objective functions. With the development of fuzzy set theory, mathemati- 
cal programming models can be examined in a fuzzy environment as in other scientific 
fields. Mathematical programming models are handled in different ways as fuzzy deci- 
sion variable, fuzzy goal, and fuzzy parameter. The fuzzy goal and fuzzy parameter 
models are used in the literature for de novo programming, also known as optimal sys- 
tem design. In this study, a new fuzzy de novo programming approach was proposed 
by using positive and negative ideal solutions. In the proposed approach, using fuzzy 
unit prices of resources and fuzzy resource amounts of constraints, the fuzzy budget 
was created. The solution steps of the developed approach are given step by step on an 
illustrative example accepted in the literature. 


Keywords: Multiobjective de novo programming, multiobjective linear program- 
ming, fuzzy de novo programming, optimal system design 


2.1 Introduction 


In classical linear programming (LP) problems, there is a goal to be achieved based 
on certain decision variables and constraints. The solution of such problems deter- 
mines the decision variable values as well as the objective function values based on 
them. On the other hand, it is possible to define the resource amount in each con- 
straint using the values of decision variables. Instead of focusing on objective func- 
tion values after solving the problem, it may answer significant problems to 
concentrate on the utilization amount of constraint resources, such as “Are resource 
amounts utilized optimally?” The answer to this problem would result in the interpre- 
tation of all the functions in the problem. When a problem is analyzed in terms of the 
amount of resource utilization, it mostly leads to an excess amount of resources and/ 
or a need for more. Either of them is a reason that hinders the optimal utilization of 
resources for production and also causes a downfall in the performance level of 
objective function. All constraints must be active for classical LP to determine an 
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optimal result. In other words, all constraint resources must be utilized at full capacity. 
Then, such a problem is named as an optimal model in terms of resource utilization. In 
multiobjective linear programming (MOLP) problems, the conditions described above 
are taken to a different dimension with “the increase in the number of goal function." 
It is because each goal function is achieved at different values of different variables in 
the solution of MOLP problems. That would lead to differences in resource utilization 
amounts of constraints in terms of each objective function. 

It is an expected case in MODM problems that constraints have limited amounts 
and goal functions are not homogeneous. Therefore, there is no concept of “optimiza- 
tion" in MODM (multiple objective decision-making) problems. That is why it is almost 
impossible for all goal functions to be realized optimally and simultaneously [1]. A sig- 
nificant feature of MODM methods is the concept of “optimality,” which is defined for 
single objective optimization problems but not valid for MODM problems. Instead of 
the concept of “optimality,” decision makers therefore use the concept of pareto- 
optimality (alternatively, efficient, nondominated, or noninferior), proposed by Italian 
economist Pareto, for MODM problems [2]. The concept of “Pareto-optimality” is 
named as “pareto-optimal” solutions for MODM problems as defined by Yu [3]. The set 
number of nondominated solutions is quite high in MOLP problems. In a MODM prob- 
lem with p-number of goal functions, the concept of *optimality" for any goal function 
may not be applicable to other p-1 goal functions. The fundamental process in evaluat- 
ing nondominated solutions is how to determine the proximity to ideal solution points 
(Cohon, 1978). The methods in MODM classification are used in the solution of MOLP 
problems. It is possible to reach a satisfactory solution among nondominated solutions 
by using the interactive methods such as step method, Geoffrion- Dyer-Freinberg 
method, and Zoints and Wallenius method [4]. Another process in evaluating nondomi- 
nated solutions is the analysis based on positive and negative ideal solutions. The util- 
ity approach, goal programming, interactive approaches, and fuzzy approach models 
can be applied to reach a compromise solution [5]. Additionally, Tabocanon (1988) de- 
scribed solutions of MOLP problems using global criterion method, interactive meth- 
ods, and compromise programming methods. 

The studies by Zeleny [6, 7] provided a different viewpoint in solving single objec- 
tive linear programming (SOLP)/MOLP problems. Named as “optimal system design,” 
this approach enables the reorganization of resource amounts using “budget.” Also 
known as de novo programming, “optimal system design" can be regarded as a spe- 
cial extension of LP. The most significant feature of de novo programming is that it 
can reorganize resource amounts of constraints based on the budget constraint. 
Resource amounts of constraints are utilized at full capacity with this reorganization. 
It brings an optimal model in terms of resource utilization. Therefore, de novo pro- 
gramming approach may be named as a redesigned SOLP/MOLP model. *De novo" 
hypothesis may be easily applied to MOLP problems as well. However, “optimal sys- 
tem design" application for MOLP inhabits a long process. 
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This study proposes a new approach for solutions of de novo programming 
problems in fuzzy environment. This approach is built on the approach by Li and 
Lee [8], which uses positive and negative ideal solutions. The proposed approach is 
analyzed with an illustrative example. 


2.2 An overview of the multiobjective de novo 
programming process 


The first study about de novo programming, which designs an optimal system in- 
stead of optimizing a given system, was carried out by Zeleny [6]. Today’s produc- 
tion and management systems are necessarily becoming more and more flexible. 
They need to be designed and redesigned rapidly, disassembled and reassembled 
again, which requires continuous reorganization of resources for ensuring opti- 
mality [9]. As all systems are created within their limits, they lack alternatives, op- 
tions, and design variables in their creation environment. Therefore, when a 
system is to be redesigned, reconfigured, or optimized, its limitations and bound- 
aries have to be done so as well. It is not sufficient to reconfigure it based on a 
given system with its own priorities and options. Therefore, system design re- 
quires creations of alternatives instead of selection [10]. In contrary to optimizing 
a given system as standard approaches do, de novo enables constraint variables 
and works to reach a solution that is more reassuring than those within the sys- 
tem’s limitations. This approach must be implemented prior to the onset of pro- 
duction because optimal production is necessarily dependent on optimal 
utilization of raw materials which is only possible with an optimal production 
plan (Babic and Pavic, 1996). 


2.2.1 De novo programming literature survey 


It is possible to classify the scientific articles on de novo programming based on 
their theoretical and application aspects. This subtitle summarizes the studies that 
theoretically contribute to de novo programming problems. 

As stated above, Zeleny [6] was the first to propose de novo programming ap- 
proach, which designs an optimal system. Zeleny [7] applied to a classical LP 
problem. These two studies make up the foundation of de novo programming 
problem. Zeleny [11] summarized the approach and methodology of de novo pro- 
gramming, which is related to extensive LP and LP-related problems. Implications 
for solving traditional LP problems are also given in his study. It also provided 
new knowledge on the “new concept” of optimality and of MCDM. Hessel and 
Zeleny [12] used the concept of metaoptimum to explain how the system design 
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can be carried out for multiobjective de Novo programming. Zeleny (1989) proposed an 
ERA (External Reconstruction Approach) method with very simple philosophy. The 
right side constants of constraints are organized using its specific algorithm. This orga- 
nization enables the full utilization of constraint constants. It can be regarded as a dif- 
ferent aspect of de novo approach. This approach does not use the budget constraint. 
Zeleny [10] suggested the *optimal path method" to create the design of an optimal 
system to conduct the solution of a de novo problem through an ideal system design. 
Shi [13] proposed in his study that several optimum path ratios could be used as alter- 
natives to select the optimal system design to solve multiobjective de novo program- 
ming problems. Apart from “optimum path ratios” proposed by Zeleny [10] and Shi 
[13], there is no certain method to be used in solving multiobjective de novo program- 
ming problems. Some studies make use of MODM methods to solve multiobjective de 
novo programming problems. Tabocanon (1988) solved multiobjective de novo pro- 
gramming problems by using global criterion method, compromise constraint method, 
and step method. Umarusman [14] demonstrated how multiobjective de novo program- 
ming problems could be solved by using minmax goal programming. Umarusman and 
Türkmen [14] suggested the global criteria method, which relies on positive ideal solu- 
tions to solve multicriteria de novo programming problems. Zhuang and Hocine [15] 
proposed the meta goal programming for the solutions of multicriteria de novo pro- 
gramming problems. Banik and Bhattacharya (2018) suggested weighted goal program- 
ming method to solve multiobjective de novo programming problems. 

The short literature survey above lists some of the articles that investigate the solu- 
tion of multiobjective de novo programming problems in classical environments. In ad- 
dition to classical decision-making process, there are significant studies on de novo 
programming problems in fuzzy environments as well. Below are the theoretical stud- 
ies on the de novo programming in fuzzy environments. First, Lai and Hwang [16] 
made use of Chanas' (1983) nonsymmetric approach to solve a single criterion de novo 
programming problem. Li and Lee [17] were the first to propose a fuzzy approach to 
solve a multicriteria de novo programming problem using “the possibility concept of 
fuzzy set theory." In their study as a continuation of the former, Li and Lee [18] sug- 
gested a “two-step fuzzy approach” based on positive ideal solutions and negative 
ideal solutions. They proposed that their approach is very efficient and can be applied 
to most multicriteria de novo programming problems. In their next article which is con- 
sidered to be a continuation of their first study where they analyzed multicriteria de 
novo programming in a fuzzy logic frame, Li and Lee [8] studied fuzzy goals and fuzzy 
coefficients simultaneously and proposed a different approach. Bhattacharya and 
Chakraborty [19] suggested an alternative method for general multiobjective de novo 
programming problems under fuzzy environment, which uses Luhandjula's compensa- 
tory H O-operator. Banik and Bhattacharya (2018) put forward a one-step method to 
solve multiobjective de novo programming problem, which makes use of min-max 
goal programming technique where all parameters are fuzzy numbers. Umarusman 
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[20] suggested for solving multiobjective de novo programming problem in fuzzy set- 
ting using Wu and Guu (2001) approach. 


2.2.2 Multiobjective de novo programming formulation 


Similar to MOLP problems, multiple objective de novo programming problems are 
handled with more than one and conflicting objective functions based on the same 
constraint functions. Multiple objective de novo programming is demonstrated 
mathematically as follows [10]: 


n 
Max Zk: 5 Cup 
IEN 


n 
Min Ws: 5 Lat 
j=l 


Subject to (2.1) 
m 
3 AXGB 
jel 
xj20, k=1,2,...,1 s=1,2,...,7, and j=1,2,...,n 


where 

Aj - Y i-,piaj is the total budget with B utilization, and P; is the unit price of 
ith resource. ) symbol means that the budget constraint is “<” or “=” type. This 
decision rests on the decision maker. “=” allows the total use of budget constraint, 
whereas “<” means that it should be less than the current budget. 

De novo formulation provides not only the best mixture of outputs but also 
the best mixture of the inputs (Tabucanon, 1988). One of the advantages of this 
approach is that it enables the adjustment of resource constraints in multicrite- 
ria decision-making problems and that it makes the initial solution more conve- 
nient at the same or less cost. Resource amounts of constraints are utilized at 
full capacity with de novo solution. It enables slack, surplus, and artificial varia- 
bles to be equal to zero. When the solution of de novo programming problems is 
investigated from MOLP perspective, positive and negative ideal solutions [21] 
must be determined for each objective functions. The positive ideal solution set 
is indicated as I = {Z}, Z}... ZiWi, W3» ..., W;) and negative ideal solution set 
as I^ - (Z,,Z;,..., Zi :IWLLW; , ..., W; ). Then the solution of eq. (2.1) may be 
done in terms of *optimal path ratio" as proposed by Zeleny [10] and Shi [13], 
based on positive and negative ideal solutions. 
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2.3 An approach proposal for fuzzy multiobjective 
de novo programming 


Conventional optimization methods evaluate parameters based on precise data. 
However, decision makers want to know which changes in parameters would result 
in which changes in the solution of the problem. It is because the precise parame- 
ters used by decision makers are only valid for “that moment." It becomes harder to 
determine optimal values of goals when short- and long-term changes in parame- 
ters are accounted for. The development of fuzzy set theory by Zadeh [22] and fur- 
ther improvements made it possible to acquire more flexible and informative 
results that present solution by forming given parameters structured as fuzzy 
ranges. Fuzzy set theory, which helps decision makers in fuzzy decision environ- 
ments, may be easily applied to multiobjective de novo programming problems. 
This chapter presents a new methodology based on the proposal by Li and Lee [8], 
which accounts for both goals and coefficients at the same time. 

The model by Li and Lee [8], which presents fuzzy goals and fuzzy coefficients 
(Gig, Čsj» à, pi, and B) formed by using positive ideal and negative ideal solution in- 
formation is here solved as a two-step fuzzy approach. Consider the fuzzy multicri- 
teria de novo programming problem: 


n 
Max Zut? = 5 CigXj 
Zei 


n 
Min W,(x) = 5 CgXj 
j=l 


Subject to (2.2) 


xj20, i212, ...,m, j=1,2,...,n,k=1,2,...,1 and s=1,2,...,r. 


Parameters Cu, Co, Oo, Pi, and Bin eq. (2.11) are fuzzy numbers. The fuzzy parameter 
model proposed by Li and Lee [8] may be summarized as follows: different ranges are 
used for maximization or minimization goals in eq. (2.11). For a maximization prob- 
lem, if’, f?) type ranges may be used for the fuzzy coefficients in the goal function. 
On the other hand, fuzzy coefficients in the goal function must be set in (5 f! range 
for the minimization problem. Then, the fuzzy coefficients aj, p; are (jo. f! type, and 
the fuzzy budget B is (f°, f'] type. Therefore, an optimal system design at the safety 
level a is formed according to the highest profit, lower cost, highest amount of 
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available, lowest resources price, and smallest units of materials to be used for each 
product. Membership functions may be defined for these ranges using linear form. 
Next, the compromise solution is conducted for goal functions based on the defined 
a — cut value. According to Li and Lee [8], (x)g be a solution of (2.2) where a € [0; 1] 
denotes the possibility to which all fuzzy parameters satisfy, and f € [0; 1] denotes 
the compromise to which the solution satisfies the fuzzy goals at given possibility a. 
À denotes the overall satisfaction for the solution (is, When a = =A the optimal solu- 
tion is obtained. If the results do not satisfy a= B, choose another a and solve phase 
I. The two-step approach [8] is described as follows. In the first step, 


Max B 
Subject to (2.3) 
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Analyzed by Li and Lee [8], in fuzzy environment based on de novo programming, 
model parameters and investment budget were made more flexible according to the 
features given below for the product mixture model: 

a. When the constraint resources (b;) were constant, the fuzzy investment budget 
(B) was determined using the fuzzy unit prices of resources (p;). In an investment 
budget like this one, the fuzzy investment budget was formed based only on the 
greatest tolerance value of unit prices. 

b. It was estimated that the resource values that were needed to produce each 
product would decrease based on the smallest tolerance value, and the unit 
price for the resource amount of each product was evaluated based on the 
smallest acceptable unit price. In other words, the cost of resource amount was 
lowered in terms of resources and the lowest acceptable level of unit price of 
resource amount. 

c. Finally, based on the order of goal function, if the goal function was maximiza- 
tion, goal coefficients were increased based on the tolerance value acceptable 
by the decision maker; and if goal function was minimization, those coeffi- 
cients were decreased. 

d. The model that was formed based on the factors above was solved according to 
the *1— a" risk protection (trust) level. When factors of risky and risk-free pro- 
duction are considered together, it must be determined in the beginning for 
which factors the decision maker would take risks. Then, the problem is to be 
modeled this way and solved according to an acceptable “1-a” risk protection 
level. 


Rather than modeling the present production models at an acceptable risk level, 

the model should be formed to produce more efficiently by lowering risks in the be- 

ginning of product mixture problems. Based on de novo hypothesis, the increase or 
decrease of the budget will directly affect production as resource amounts are deter- 
mined accordingly. Therefore, the increase or decrease of resource unit prices is sig- 
nificant for the budget. Additionally, the definition of resource amounts in a range 
will reform the budget. Consequently, the appropriate determination of both re- 
source amounts and their unit prices results in a reasonable budget realization. 

The risk factors in the studies model and assumptions on their analyses in 
fuzzy environment are listed as follows: 

1. When allocating a budget for investment or resources, the constraint resources 
should be analyzed in a range based on the de novo programming. The budget is 
formed based on the constraint resource amounts. As the budget cannot exceed 
the acceptable tolerance value of constraint resources, there is less risk in terms 
of the budget, which is formed according to the expectations of the decision 
maker instead of surpassing. Nevertheless, businesses may only increase their 
budget up to a certain level. 
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2. In the fuzzy de novo model proposed by Lee and Li [8], utilization amounts for 
products were decreased and unit prices were evaluated according to the small- 
est acceptable value and both of which increased production. On the other 
hand, the risk factors in production are lowered when the evaluation is based 
solely on unit prices as resource utilization amounts for each product are fixed. 
In other words, only the acceptable minimum value of unit prices should be 
taken into consideration when forming constraints in de novo models. 


The proposed approach uses ei eil type range for fuzzy coefficients in maximization- 
type goals and Gag f! range for fuzzy coefficients in minimization-type goals. aj coeffi- 
cient is used as the constant number. The fuzzy coefficients p; and b; are of (f° ; f 
type, and the fuzzy budget B is of if D f°) type. It means that the optimal system design 
provides lowest cost, highest profit, lowest resource unit price, and highest budget uti- 
lization at a trust level. Consider the fuzzy de novo programming model based on the 
aforementioned information: 


Subject to (2.5) 


xj20, i=1,2,...,m, j=1,2,...,n, k=1,2, ...,l ands=1,2,...,r 


Solutions are determined as positive ideal or negative ideal for each goal function 
based on model (2.5): 
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Based on the “a-cut” value determined by the decision maker, the compromise solu- 
tion for goal functions is expressed as B = min{ ys, (Zz), Ho, (Walk: El 
S-1,2,...,r). The proposed model for these defined formulations is given below. It 
makes up the first step of the two-step approach: 
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Max f 
Subject to (2.10) 
Bx Bee): Xj (Zr), 
(Z9. - (1). 


j=12,...,.nands=1,2,...,r 
After getting the optimal value of a and f, this solution is tested whether it is effi- 


cient or not. If eq. (2.11) is rewritten based on the membership functions of the pa- 
rameters, the second step of the model is formed as follows: 


= 1 l r 
Max B= SR Kä yA.) 
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k=1,2,...,1,s=1,2,...,r, andi- 12, ...,m 


The optimal system design is achieved at A = a = B, where À represents the minimum 


level of possibilities under consideration of both the fuzzy goal and the fuzzy 
coefficients. 
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2.4 Illustrative example 


The multiobjective de novo programming problem prepared by Zeleny [23] was 
used to compare the results of the proposed approach and the approach by Li and 
Lee [8]. This problem consists of three maximization-type objective functions, six 
constraint functions, and three variables as follows: 


Max Z, = 50x; + 100x2 + 17.5x3 
Max Zp = 92x, + 75x» + 50X3 
Max Z; = 25x, + 100x2 + 75X3 
Subject to (P1) 
12x, + 17x? € 1400 
3X1 + 9X2 + 8x4 < 1, 000 
10x; + 13x; + 15x3 <1, 750 
6X4 + 16x3 x 1325 
12x; + 7x3 < 900 
9.5X1 + 9.5X2 + 9.5x3 x 1, 075 


X1,X2,X32 0 


The resource unit prices for (P1) are $0.75, $0.6, $0.35, $0.5, $1.15, $0.65, respec- 
tively. The fuzzy parameter model for (P1) is as follows: 


Max Z, = 50x; + 100x + 17.5x3 
Max Z) = 92x, + 75x; + 50X3 
Max Z3 = 25x; + 100x + 75x3 
Subject to (P2) 
125 + 17x; < bi 
3X1 + 9X2 + 8x4 < 1, 000 
10x; + 13x2 + 15x3 x bs 
6x1 + 16x3 x 1325 
12x; + 7x3 < 900 
9.5xX1 + 9.5X2 + 9.5x3 < 1,075 


0.75b; + 0.6b; + 0.35b; + WICH + 115b; + 0.65D¢ < B 
Xp X2, X3, bi, bo, bs, ba, bs, be = 0 
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Defined as fuzzy triangular number, goal function coefficients are chosen for each goal 
function. They are 50 = (40,50, 80), 75=(70,75,90), and 25 = (10, 25,35). a - Cut in- 
formation for the left side and right side corresponding to these coefficients are: 


(50) * = (40+ 10a), (ai = (80 — 30a) 


a Qa 
SVT ~\R 
(75) = (70+ 5a), (75) = (90 - 15a) 
a a 
UA LR 
(25) - (10 4 15a), (5) - (35 - 10a) 
Q a 
a-Cut information for the left side and the right side are organized as follows for a 
total of five fuzzy numbers such as two triangular fuzzy numbers accepted for the 
selected right side constants (b;), the two trapezoidal fuzzy numbers accepted for 
the selected unit prices (p;), and one triangular fuzzy number: 


b, = 1400 = (1, 300, 1, 400, 1, 500); 
Lg € R 
(1400) - 1, 300 + 100a), (14000) =, (1500 - 100a) 


a a 


b; = 1750 = (1,500, 1, 750, 2, 000); 

— NL —~—\R 
(1750) = (1,500 + 100a), (1750) = (2,000 - 250a) 
Qa 


a 


B, = 0.75 = (0.6, 0.7.0.8, 1.15); 
GL —~—\R 
(075) = (0.6 0.10), (075) = (1.15 - 0.35a) 


[4 a 


P = 0.35 = (0.3, 0.35, 0.5); 
ENGL —~—\R 
(035) = (0.3 + 0.054), (035) = (0.5 — 0.154) 


a a 
ps =1.15= (1,1.15, 1.2, 1.4); 
AE coy R 
(15) = (14 0.150), (115) - (14 - 02a) 
a a 
Using this formulation, the calculation of the left side and the right side of the bud- 
get is as follows. For the left side of the constraint: 
x E 
(075) * (12x, + 17x2) + 0.6* (3x, + 9x; + 8x3) + (0:35) * (10x; + 13x; + 15x) 
a a 
aes 
+0.5*(6x; + 16x3) + (115) *(12x5 + 7x3) + 0.65* (9.5, + 9.5x; + 9.5X3 
a 


Or 


26 —— Nurullah Umarusman 


= (0.6 + 0.1a)* (12x; + 17x2) + 0.6* (3x3 + 9X2 + 8x3) 
+ (0.3 + 0.05a)* (10x, + 13x5 + 15x3) 
+ 0.5* (6x + 16x3) + (1+ 0.15a)* (12x2 + 7x3) + 0.65*(9.5x1 + 9.5xX2 + 9.5xs 
is acquired. The final version of the calculation is 
(21.175 + 1.7a)x; + (37.675 + 4, 15a) x + (30.475 +1, 8a)x3 


For the right side of the budget constraint: 


B= (1,500 -100a)* (1.15 — 0.35a) + 1, 600*0.6 + (2, 000 — 250a)*(0.5 — 0.15a) 


+ 1, 325*0.5 + 900*(1.4 — 0.2a) + 1, 075*0.65; B - 6,306.25 - 1, 245a + 40.50? 


Objective functions are organized as 
Max Z, = (80 - 30a)x; + 100x + 17.5x3 
Max Ž = 92x, + (90 - 15a)x; + 50x3 
Max Z; = (75 - 50a)x, + 100x; + 75x3 
(P1) is formed as below after these processes: 
Max Z, = (80 - 30a)x; + 100x; + 17.5x3 
Max Z» = 92x, + (90 — 15a) x2 + 50x3 
Max Z3 = (75 - 50a)x, + 100x + 75x3 
Subject to (P3) 
(21.175 + 1.7a)x; + (37.675 + 4, 15a) x2 + (30.475 +1, 8a)x3 
< 6, 306.25 — 1, 245a + 40.50? 
X1, X5, X320,a € [0,1] 
First, if a parametric solution is made for a = 1 (P3), it is obtained as follows: 
Max Z; = 50x; + 100x + 17.5x3 
Max Zp = 92x, + 75X2 + 50x3 
Max Z3 = 25x, + 100x + 75x3 
(P3.1) 
Subject to 
22.875X, + 41.825x; + 32.275x3 € 5,101.75 
X1, X2, X3 20 


acquired. 
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Positive ideal solutions for each objective function from (P3.1) are calculated as 
(Z),., = 13,528.99; (Z2), ., 2 22,757.68; (Z3),_, = 14,621.45 


and negative ideal solutions are calculated as 


The first phase is written for (P3.1) using eq. (2.10) as follows: 


Max f 
Subject to (P3.2) 


B " [50x1 RE 100x2 + 17.5x3 -0 
Sd 13, 528.995 - 0 
Be [92x4 + 75x? + 50x3 - O 
^t  22757.68-0 
[25x; + 100x + 75x3 - O 
p< : 
"t 1462145-0 


22.875x4 + 41.825x» + 32.2/5x3 € 5,101.75 


X1, X2, X3 2 0,ß € [0; 1] 
Values obtained from the solution of (P3.2) are given in Table 2.1. 


Table 2.1: Decision variables of Phase | for (P3.1). 


Decision variables Z3 Za Z3 
X 101.208 101.208 101.208 
X) 66.625 66.625 66.625 
Xs 0 0 0 


Objective function value 11,722.9 14,308.01 10,457.8 


Information in Table 2.1 is the result of the solution of (P3.1)'s first phase. The objec- 
tive function value of (P3.1) is B = 0.6287 that is for a- 1, fj - 0.6287. According to 
this result, the solution of (P3) will take place in the range of [0.6287; 1]. Because of 
a + p, the solution should be continued by selecting a new a. The result obtained 
for the selected a = 0.71999 is optimal. First, if (P3) is taken as a= 0.71999 
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Max Zi; = 58.4x, + 100x2 + 17.5x3 
Max Z2 = 92x, + Ki e SR 50x3 
Max Z3 = 27.8x1 + 100x + 75X3 
(P3.3) 
Subject to 
22.399x, + 40.633x; + 28.521x3 x 5, 430.845 
X1, X2,X3 20 


is obtained. First, positive and negative ideal solutions for (P3.3) are given as 
follows: 


D 


(Z1)a-071999 = 14 159.62; (2), sie = 22, 306.25; (23). age = 14, 281.17 
and 
(Za) a-071999 =0; (22) amas 70; OR 071999 "H 
Phase I is arranged using these solutions as follows: 
Phase I 
Max f 
Subject to 
58.4x + 100x; + 17.5x3 — 14, 159.628 > 0 
92x1 + 79.2x; + 50x3 — 22, 306.258 > 0 
27.8X; + 100x + 75x3 — 14, 281.178 = 0 
Subject to (P3.4) 


22.399, + 40.633X2 + 28.521x3 < 5, 430.845 


X31, X2, X3 2 0, B € [0; 1] 


f — 0.71999, xı =113.057869, x)= 71.332710, x3 - 0 are obtained from (P3.4)’ solution. 
In the first phase, f =0.71999 is obtained for a= 0.71999. It is now tested in 
phase 2 whether the solution of (P3.4) is efficient: 
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Max p= 5 (B, +B, +B) 


Subject to 


B<B3< 


[58.4x4 + 100x2 + 17.5x3 - O 


14,159.62- 0 


22, 306.25 - 0 
[27.8x1 + 100x2 + 75x3 — O 


[92x, + 79.2X2 + 50x3 — °] 


14, 281.17 - 0 


| 


22.399x, + 40.633x + 28.521x3 < 5, 430.845 
B, = 0.71999 
f,2 0.71999 
B3 > 0.71999 


Xy, X2 X3 20,8 € [031] 
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(P3.5) 


Variable values obtained from (P3.5)'s solution created for phase 2 are specified as 


B - 0.80334, f, — 0.970072, B, = 0.71999, B, = 0.71999 


X1 2113.057861, x2 = 71.332710, and x32 0 


This result shows that values obtained from phase 1 are efficient. Depending on the 
different values of a, variable and objective function values are given in Table 2.2. 


Table 2.2: The results of (P3). 


a B 
1.0 0.6973268 
0.75 0.7169044 
0.73 0.7191670 
0.725 0.7188874 
0.71999 0.71999 


Zi 


13002,15 
12664.65 
12743.9 

12761.55 
12786.05 


Z 


15869.4 

15584.99 
15673.59 
15718.16 
15751.14 


Z3 


11598.99 
11268.66 
11341.23 
11351,69 
11372,84 


Budget 


5101.75 

5395.281 
5418,982 
5424.913 
5430.845 


The optimal result is obtained at a = f = 0.71999. At this point, resource utilization 
amounts according to decision variable values are 2569.328, 981.159, 2,057.886, 


678.342, 855.984, and 1,751.696, respectively. 
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2.5 Conclusions 


As the main objective in traditional production processes is maximum profit and 
minimum production process, conciliatory or satisfactory solution is sought instead 
of optimal solution due to the multiple objective functions. In multiobjective de 
novo programming problems, it is possible to specify the optimal solutions at the 
end of similar processes. For example, if multiobjective de novo programming prob- 
lem is solved through compromise programming, the result can be named as com- 
promise optimal system design, or if it is solved through goal programming, the 
result can be named as satisfactory optimal system design. In this study, multiob- 
jective de novo programming problem is solved because the solution of this prob- 
lem is made in terms of fuzzy set theory, and is named as fuzzy optimal system 
design. However, the concept of compromise solution is used in some studies in- 
stead of fuzzy solution. 

In this study, the model offered for fuzzy optimal system design comes into 
prominence with some features. First of them is although it is possible to define aj 
parameters as a fuzzy number, it is wrong to fuzzify this parameter in the produc- 
tion process. This is because the change in the parameter aj causes the amount 
used to produce the product during the production process and this disrupts the 
generated model. Also, p; and b; fuzzy parameters were used while setting a fuzzy 
budget. This situation caused the upper limit of the fuzzy b; determined by the deci- 
sion maker to directly affect the budget while creating the fuzzy budget. The pro- 
posed approach can be accepted as an extension of Li and Lee’s [8] approach 
because the solution procedure of the proposed approach is the same as that of Li 
and Lee’s approach. Although it is possible to compare the results of Li and Lee’s 
[8] approach with the result of the proposed approach, it is not recommended be- 
cause of the differences explained earlier. Fuzzy parameters are different from each 
other, even though the solution of both approaches is the same. 
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3 A probabilistic two-level programming 
in noncooperative game under fuzzy 
environment: comparative study 


Abstract: This chapter develops a computational algorithm to solve a stochastic two- 
level programming in a noncooperative game using a fuzzy optimization technique. It 
also includes all parameters such as fuzzy variables except the constraints on the right 
side that follow an extreme value distribution. Mellin transformation is considered to 
convert the fuzzy numbers into the crisp value. Based on stochastic programming, 
probabilistic constraints are converted into the deterministic form. Furthermore, as the 
fuzzy programming is used for the degree of acceptance, we present a comparative 
study of nonlinear and linear membership functions for the acceptance corresponding 
to LINGO 15.0 iterative scheme and genetic algorithm to see its impact on optimization. 
A numerical experiment is considered in order to present the applicability and feasibil- 
ity of the proposed algorithm. Finally, conclusion about the findings and outlook are 
presented. 


Keywords: Bilevel programming, fuzzy number, GA, membership function, 
Stackelberg game, stochastic programming, satisfactory solution 


3.1 Introduction 


Two-level programming problem (TLPP) consists of two levels in any hierarchical 
decision process, and large numbers of research works have been investigated on 
TLPP in the past few decades ([19],[20],[24]). In TLPP, each decision maker (DM) 
optimizes his/her own objectives without considering the other one. Stackelberg 
game was first introduced by an economist H. von Stackelberg to solve noncoop- 
erative decision problems. In Stackelberg game, the leader has the capability to 
enforce his/her decision on the followers. Many researchers have carried out in- 
vestigations on the Stackelberg game (cf [1—4]) Ref. [17]. In real-life-based applica- 
tions, when the input data of any mathematical model are inexistent or scarce, 
data are considered as uncertain. To tackle the uncertainty problem, two principal 
resources say, fuzziness and randomness, are considered. Fuzzy programming (cf 
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[5, 6]) developed by characterizing the input data as fuzzy sets, and stochastic 
programming is also developed to express the input data as random variables. 
Using the Mellin transformation, all fuzzy numbers (triangular or trapezoidal) are 
reduced to a crisp form. The constraint parameters on the right side follow an ex- 
treme value distribution Ref. [22]. Then, using a stochastic programming, these 
constraints are converted to the deterministic form and the equivalent crisp prob- 
lems are solved by LINGO 15.0 iterative scheme and the results are compared. This 
chapter includes the following: 
- Mellin transformation is used for the fuzzy numbers (triangular or trapezoidal) 
into crisp value. 
- The constraint parameters on the right side are considered to be an extreme 
value distribution. 
— A computational algorithm is developed to solve a noncooperative game. 


We develop a computational algorithm for solving a stochastic TLPP in a noncoop- 
erative game involving an extreme value distribution by fuzzy programming. Our 
aim is to study the impacts of different types of membership functions. For that in- 
stance, we consider a linear function and a nonlinear function. The research works 
of several authors related to this area are shown in Table 3.1. 


Table 3.1: Research works of several authors related to this area. 


Author(s) Real number Fuzzy number Fuzzy random Multichoice 


Anandalingam and Apprey [1] Yes 


Anandalingam and White [7] Yes 
Hamidi and Mishmast [8] Yes 

Maity and Roy [2] Yes Yes 
Roy and Maiti [4] Yes Yes 
Sakawa and Matsui [9] Yes 

Singh and Chakraborty [10] Yes 

Yeh et al. [11] Yes 

Youness et al. [12] Yes 

Zheng et al. [13] Yes 

Our chapter Yes Yes 


The sequence of the chapter is given as follows: Section 3.2 introduces the basic 
knowledge of fuzzy set. In Section 3.3, model formulation of stochastic TLPP in a 
noncooperative game involving an extreme value distribution is presented. The 
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computation algorithm of fuzzy optimization is discussed in Section 3.4. A numeri- 
cal example is presented in Section 3.5. Results are reflected in Section 3.6. The con- 
clusion and an outlook of future research directions are described in Section 3.7. 


3.2 Preliminaries 


Here, we recall the useful definitions of fuzzy set and related theorems that are to 
be required for our subsequent developments. 


Definition 3.1 ([11]). Let A be a fuzzy set defined on X, the universal set. Then 
A= {(x,yq(x)):x € X], where pj: X — [0,1] is the membership function such that 
Hy (X) =Oif xA, pa(x)=1ifx € A. 


Definition 3.2 ([18]). Let A bea fuzzy set defined on R, the set of real numbers, 
which is said to be a fuzzy number if its membership function satisfies the following 
conditions: 

(a) q(x): R — [0,1] is continuous, 

(b) uj(x) 0 Vx € (- es, c]u [b, ee), 

(c) Hat) is strictly increasing on [c, d] and strictly decreasing on [a, b], (c « d « a « b), 
(d) wa(x)=p Vx € [d a], where O«p «1. 


Definition 3.3 ([14]). Let B= (bi, bj, b3) be a fuzzy number. Then B is said to be a 
triangular fuzzy number if ug(x) is given by 


(E) if bi «x b; 


ug(x)-7 (Zz): if b; «x € b; (3.1) 
0, if otherwise 


Definition 3.4 ([14]). Let C= (C1, C2, C3, C4) be a fuzzy number. Then C is said to be 
a trapezoidal fuzzy number if u¿(x) is given by 


(£2). ife; «xc; 


1, if@<x<c3 (3.2) 


ax), ifc £x x c4 


O, if otherwise 
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Definition 3.5 ([15]). Let g(x) be a probability density function (pdf). Then Mellin 
transform My (r) of g(x) is given by 


provided the integral exists. 


3.3 Model formulation 


This chapter discusses on a stochastic TLPP for a noncooperative game in which 
cost coefficients of both level objectives are either triangular or trapezoidal fuzzy 
numbers. Generally, the constraint coefficients in a TLPP are considered as crisp. 
But in real-life situations, the constraint parameters are imprecise. So, we consider 
that the constraint parameters are of two different natures. The left side of the con- 
straints is a fuzzy number, whereas the right side follows an extreme value distribu- 
tion. The following assumptions are considered to design the model: 

— Gu: fuzzy number (triangular or trapezoidal) for leader objectives, j - 1,2, ...,n; 

- Cy: fuzzy number (triangular or trapezoidal) for follower objectives, j= 1,2, ...,n; 

- ay: a mxn matrix with fuzzy numbers; 

— bi: random variable follows an extreme value distribution. 


Therefore, the model on a stochastic TLPP for a noncooperative game involving ex- 
treme value distribution is defined as follows: 


P1 
n 
maximize Z(x)= 3" Gu 


n 
maximize Z)(Xx) = 3. Can 
j=1 (3.3) 


n ~ = 
subject to pr ys <B,) 21-7 i=1,2,...,m 
ja 


xj20 Vj and O<y,<1 Vi 


3.3.1 Deterministic form of a fuzzy number 


Here, we present a defuzzification procedure of a triangular or trapezoidal fuzzy 
number. 
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3.3.1.1 Triangular fuzzy number 


Let a= (a1, 5, a3) be the triangular fuzzy number. Then the crisp value is obtained by 
finding the expected value using the pdf corresponding to the membership function. 

Therefore, the proportional pdf corresponding to a is defined as follows: 
fa(x) = du, (x), where p; (x) is given by eq. (3.1) . 

To find d such that TT fa(x) dx=1, we get d= SE 

Then the proportional pdf corresponding to à is given by 


2(x-a1) 


(aa) (45-41) > if ay < x < a2 
fao) = 2(a3 -x) ; 
a EEEN if a2 <X< a3 
0, if otherwise 


According to Definition 3.5, we obtain 


2(x- a1) r4 2(a3 - x) 
EE il (a; — a5)(as - a1) 


co a» 
Ms(r)- Je ps dx = E dr 
0 ai 


On integration, we obtain 


mo- A. [205 =) _ alas ~ a 
T (as-aj)r(r-1) | a-a a-a; 


Therefore, the mean (E[a]) and variance (07) corresponding to à are obtained as 


a 
follows: 
= dı + A2 + A3 
3 
dj + a3 + d$ — A102 — Q203 — A303 
18 


Then the crisp value of a is given by 


, dı + d2 + A3 


Crisp(ã) = a'- 2 (3.4) 


3.3.1.2 Trapezoidal fuzzy number 


Let b- (bi, b», b3) be the trapezoidal fuzzy number, and the proportional pdf corre- 
sponding to b is defined as follows: f;(x) = cu; (x), where uj (x) is given by eq. (3.2) 
and c is given by relation 
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| f(x) dx=1 (3.5) 


From eq. (3.5), we obtain c= i poe 


Therefore, the proportional pdf corresponding to B is given by 


2(x — by) : 
(bz - by) (ba + b3 -b2 - bi)" HUM SxS) 
2 : 
(by, +b3 - b; - bi)" if by Eé b; 
d (04-2) 
2 by -X : 
(ba - b3) (b4 *b3 - bo - b)" if b3<x<by 
O, if otherwise 


According to Definition 3.5, we obtain 


M;(r)= Ra 3 (x) dx 
0 
by b> 
= 2 x-b,) | a 2 
= | xo! dx rt jg 
| (b; — by) (bg + b; - b2 - b1) ; (by + b; - b; - bi) 
bi b) 
by 
E 2(b, - x) 
x dx 
y | (ba — 3) (ba + bz- b; - by) 
3 


On integration, we obtain 


2 


Mi(r)- (by + b3 - b; - byyr(r 1) 


(bit-t) (EIN 
b, — b; by - by 


Then, the mean (E DN and variance (03) corresponding to b are obtained as 
follows: 


EN =M,(2) = 5 [Oih + b+ by) + (biba — bsba) | 


(by 2 b3 = b; = bi) 
(bi + b2) (b2  b2) - (bs + ba) (b2 ed 
(by bs b3 = by = bi) 


dw Dat s |t + geen, 


Then the crisp value of b is given by 


(3.6) 


Crisp(b) = b= 5 | (bi baba + by) + (rb: "bb, | 


(by + bs - b-b) 
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Using eq. (3.4) (or eq. (3.6)), model 1 is transformed into the following model with 
the probabilistic constraint as follows: 


P2 


n 
maximize Z;(x)= M Cu 
j=l 


n 
maximize | Zi(x) = 3 Can 
j=l 


n 
subject to (3e) 21-7 i=1,2,...,m 


j=l 


xjizO Vj and O«y;«1 Vi 


3.3.2 Deterministic form of probabilistic constraint 


Here, B; Vi are independent random variables that follow an extreme value distri- 
bution with 6; and f; as location and scale parameters, respectively, where the aspi- 
ration level is y,, O < y; < 1. 

Now, the pdf of B; Viis defined as follows: 


fes gee(- Ce) (CC CR) on 


-oo«Bj«oo and 6;>0, B;>0 Vi 


Therefore, eq. (3.3) can be written as follows: 


| oe: G8) 
3. Am 
Jzd 


Using eq. (3.7), the above eq. (3.8) can be written as 
| ie tel Aal Pollen on 
i Bi Bi 


3. Au 
j=l 
Bi - 6; 
Let us take exp{| —{—~—]]=q. Then eq. (3.9) reduces to 
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On integration, we get 


St zii (zm) a, 


Taking logarithm of both sides two times and rearranging, we get 


Ye x6; - B; In( - In(y,))] (3.10) 


Then we obtained the following deterministic model: 


P3 


maximize —Zij(x 2 Cu 
maximize Z2(x 53 Ca 


subject to ya iG € < ôi - fi [In(- In(y,))] 


130 Vj and O<y,<1, 6;>0, Bj»0 Vi 


3.4 Computational algorithm 


After getting the deterministic form of P1, we develop the following two algorithms 
corresponding to the satisfactory solution using fuzzy approach. The following 
steps of fuzzy programming are given: 


Algorithm 3.1 (For linear membership function). 


Step 4.1: 


Step 4.2: 


Independently solve the leader and follower objectives with respect to 
the constraints and assign the best (Z2) and worst (Zl) solutions for 
each objective function, s - 1, 2. 

A linear membership function is considered for each objective function 
as follows: 


1, if Z;(x)Z2 
-z1 ] 
AZ. 9))- KE if Z? < Z,(x) < ZQ0, 
0, if Z;(x) « Zl, 


Again, considering a linear membership function for the upper-level deci- 
sion variable (xı) with f tolerance value is defined by 
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xy xb +f 8 
(5 p—) ifxi-fsxsx} 


Ha n) = (=), ifx| «x xxi f 
0, if otherwise 
where xj is the solution set of the top-level DM. 
Step 4.3: According to Zimmerman approach [16], solve the following problem: 


P4 


maximize A 


subject to p,(Zs(x))=A Vs 


x2 0 and À c [0,1] 


If the top-level DM is satisfied with the obtained optimal result, then the sol- 

utions become a satisfactory solution and go to Step 4.4. Otherwise, top- 

level DM updated his/her membership function and then go to Step 4.2. 
Step 4.4: Stop. 


Algorithm 3.2 (For nonlinear membership function). 
Repeat Step 4.1 and follow the steps. 


Step 4.2: Consider a nonlinear membership for each objective, and the top-level 
decision variable (xı) with tolerance value f is given by 


1, if Z;(x) > Z2 
Zei) -Z; 
dc v^ E 
us(Zsx)) = S iret |, if Z$<Z;(x) <Z? 
0, if Z;(x) < Zl 
(st) 
£O iff «xd 
Hy (x1) = E gl 
Been dixisxsxb4f 
0, if otherwise 


where x. is the solution of the top-level DM. After that repeat Step 4.3 and 
Step 4.4. 
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3.5 Numerical experiment 


Here, we present two numerical examples corresponding to triangular and trapezoi- 
dal fuzzy numbers and one practical application. 


3.5.1 Example 1 
Let us consider the following example: 
P5 


maximize Zj(x)- (200, 302, 388)x; + (234, 280, 347)x2 + (187, 199, 223)x3 + (207, 289, 301)x, 
maximize Z2(x) = (89,107, 112)x; + (106, 127, 132)x2 + (67, 80, 93)x3 + (73, 79, 104) x4 
subjectto (34, 42, 49) x, + (27, 37, 52)x2 + (30, 47, 53)x3 + (18, 22, 28)x4 < bi 

(32, 37, 42)x, + (18, 22, 25)x; + (25, 27, 29)x3 + (30, 32, 36)x, < b2 

(26, 36, 42)x, + (16, 20, 27)x2 + (41, 44, 49)x3 + (51, 55, 59) x, > b; 

(38, 43, 47)x1 + (52, 56, 57)x; + (61, 65, 68)x3 + (43, 45, 49) xy > by 

xii 0 vj 

The location and scale parameters and specified probability levels of bi, b2, bs, by, 
are given in Table 3.2. 


Table 3.2: Input parameters for b;, i— 1, 2, 3, 4. 


Random Location parameter (6) ^ Scale parameter (6) Specified probability levels (y;) 
parameters b; 

bı 800 200 0.5 
b 600 100 0.6 
bs 150 75 0.7 


b, 86 25 0.8 


Using eqs. (3.4) and (3.10), P5 can be written in the following model: 
P6 
maximize Zi(x) = 296.67x; + 287X2 + 203x3 + 265.67x;, 
maximize Zo(X) = 102.67X, + 121.67X2 + 80x3 + 85.33X4 
subject to 41.67X, + 38.67x» + 43.33x3 + 22.67x, x 873.3 (3.11) 
37X1 + 21.67x» + 27X3 + 32.67x4 x 667.17 (3.12) 
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34.67x + 21x? + 44.67x3 + 55X4 2 136.08 (3.13) 
42.67x1 + 55X2 + 64.67X3 + 45.67x, 2 74.10 (3.14) 
x;z0 Vj 


Using Algorithm 3.1 of computational algorithm, P6 is converted to a linear pro- 
gramming problem, that is, P7. 


P7 


maximize A 


^ Zi (x) - 618.407 
subject to (ae) BA 


Za) — 211.122 
( Xexs ) ZÀ 


X1 +X2 - 17.363 +f 
a) 


constraints (3.11) - (3.14) 
xji2 0 Vj 


Now, P7 is solved for f = 0.1 by using LINGO 15.0 iterative scheme and genetic algo- 
rithm (GA), and the optimum results are presented in Table 3.3. 


Table 3.3: Optimal results for Algorithm 3.1. 


Method A Zi LZ X1 X2 X3 Xy 
LINGO 0.998 7,335.433 2,867.078 0.127 17.297 0.0 8.783 


GA 0.899 7,338.571 2,848.670 0.216 16.549 0.0 9.504 


Also, using Algorithm 3.2 of computational algorithm, P6 is transformed into the 
following model: 


P8 
maximize H 
(A (x) = 618.407) 
H i 7 6730487 . set 
subjectto s747 Je >A 
X (2 (x) 207 
e (TXE j o-r A 
Let 
SR (s *X)-17.363 +i) 
f -r 
e -e 
Let 2A 


constraints (3.11) - (3.14) 
x20 Vj 
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Using LINGO 15.0 iterative scheme and GA, P8 is solved for r=0.5 and f=0.1, and 
the results are presented in Table 3.4. 


Table 3.4: Optimal results for Algorithm 3.2. 


Method A A LZ X4 X2 X3 X4 


LINGO 0.99999 2,717.578 1,035.560 0.00001 3.70111 3.61948 3.46520 


GA 0.99956 2,718.753 1,037.992 0.00005 3.7811 3.5989 3.3989 


3.5.2 Example 2 


Again, let us take another example: 
P9 


maximize  Zj(x)- (200,302, 388)x; + (234, 280, 347)x2 + (187, 199, 223)x3 + (207, 289, 301) x4 
maximize Z)(x) = (89,107, 112)x; + (106, 127, 132)x2 + (67, 80, 93)x3 + (73, 79, 104)x4 
subjectto (34, 42, 49)» + (27, 37, 52)x2 + (30, 47, 53)x3 + (18, 22, 28)x4 < bı 

(32, 37, 42)x, + (18, 22, 25)x; + (25, 27, 29)x3 + (30, 32, 36)x4 < b; 

(26, 36, 42)x1 + (16, 20, 27)x2 + (41, 44, 49)x3 + (51, 55, 59)x4 > bs 

(38, 43, 47)x1 + (52, 56, 57)x + (61, 65, 68)x3 + (43, 45, 49) x4 > Dy 

xjzO0 Vj 


The location and scale parameters and specified probability levels of bi, b», bs, ba 
are given in Table 3.5. 


Table 3.5: Input parameters for b;, /=1, 2, 3, 4 corresponding to Example 2. 


Random Location Scale parameter (;) Specified probability levels (y) 
parameters b; parameter (0j) 

bi 1,200 100 0.02 
b 1,400 200 0.03 
bs 200 140 0.05 


b, 500 120 0.07 
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Using eqs. (3.6) and (3.10), P9 is reduced to the following model: 
P10 


maximize Zi(x) = 15.29x,+ 24x? + 30.55x3+ 44x, 
maximize Zo(X) = Sint 15.19x; *. 6.42X3 + 5.43X;, 


subject to 16.29x; + 34.09x» + 7.8x3 + 34.21x, < 1, 063.59 


14x, + 24.4X5 + 57.22x3+ 54.43x,< 1,149.07 
25.62X, + 29.335 + 49.64x3 + 54.3x,> 615.83 


4x, + 25.81%) + 67.29x3+ 26.81x,2 814.78 
x20 Vj 
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Using Algorithms 3.1 and 3.2 of computational algorithm, optimal results are shown 


in Tables 3.6 and 3.7, respectively. 


Table 3.6: Optimal results for Algorithm 3.1. 


Method A F Zi Z2 X1 X2 
LINGO 0.998 4.4 1,045.933 549.124 40.598 10.503 
GA 0.977 4.4 1,045.983 559.914 39.985 11.897 


Table 3.7: Optimal results for Algorithm 3.2. 


Method A F r Eé Z; X X2 


LINGO 0.999 4.4 0.3 379.686 80.292 0.091 0.0 
GA 0.995 4.4 0.3 380.318 80.329 0.074 0.0 


x3 


5.667 
4.88 


X3 


12.382 
12.412 


3.5.3 Practical application (production planning problem) 


X, 


0.0 
0.0 


X, 


0.0 
0.0 


India bought crude oil from different cities of the World. That crude oil cannot be 
directly used in our daily life. So it should be refined. Under different circumstan- 
ces, the crude oil produces various by-products from different process units. India 
has eight oil refineries in different locations such as Haldia, Gujarat, Panipath, 
Paradip, Mathura, Guwahati, Barauni, Bongaigaon. Each oil refinery has different 


process units and are discussed as follows: 


1. FOB (fuel oil block): It has two atmospheric distribution units: CDU-I and II 
along with the naphtha pretreating unit, catalytic reforming unit, and kero hy- 


drodesulfurization unit. 
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2. DHDS (diesel-hydrodesulfurization): It produces Euro-III quality eco-friendly 
diesel, which reduced sulfur dioxide in vehicular emissions and as a result pol- 


lution-free environment. 


3. FCCU (fluidized catalytic cracking unit): It processes heavier feedstock, includ- 
ing a proportion of residue, and produces value-added products and contrib- 
utes to improve the distillate yield from the refinery. 
MSQIU (motor sprit quality improvement unit): It produces EURO-III petrol. 

5. OHCU (once-through hydrocracker): It comprises hydrogen generation unit, 
sulfur recovery units, hydrogen storage facility, air compressor, heat recovery 
steam generation, effluent treatment plant, and feed preparation unit. 

6. LOB (lube oil block): It consists of vacuum distillation unit, propane de- 
asphalting unit, lube hydrofinishing unit, vis-breaking unit and catalytic ISO 


dewaxing unit. 


7. DHDT (diesel hydrotreating unit): It produces EURO-IV-grade MS and HSD 


(High Speed Diesel). 


8. GT (gas turbines): It generates electricity by using gas turbine. 


All the process units are not available in all refineries. Due to fluctuation of global 
market and plant equipments, the capacity of the process units in each refinery is 
not fixed. So considering these capacities are fuzzy numbers (triangular or trapezoi- 
dal). We mainly discuss here six types of oil refineries and their capacities that are 
supplied in Tables 3.8 and 3.9, respectively. 


Table 3.8: Input data for FOB, DHDS, FCCU, and MSQIU units. 


Refineries Process units 

FOB DHDS FCCU MSQIU 
Haldia (5,5.5,7) (6,6.3,6.5) (10,12,12.3) (11,11.2,11.4) 
Gujarat (1.1,1.4,1.5) (2.1,2.3,2.6) (2.4,2.41,2.43) (2.13,2.16,2.17) 
Panipath (1.9,2.3,2.7) = (1.2,1.3,1.7) = 
Paradip (0.47,0.49,1) (0.57,1,1.1) = (0.39,0.45,0.5) 
Mathura (5.6,5.86,5.89) (6.2,6.5,7.2) (6.15,6.19,6.34) (5.1,5.15,5.24) 
Guwahati (2.12,2.15,2.17) zs (2.18,2.19,2.21) (1.89,1.92,1.95) 


Let us consider that xi, X», Xs, X4, Xs, X6, X7, and xg be the decision variables (quan- 
tity) corresponding to the process units FOB, DHDS, FCCU, MSQIU, OHCU, LOB, 
DHDT, and GT, respectively. So the problem gives a TLPP where the top-level DM 
maximizes the profit and the lower one also maximizes the products. The cost 
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Table 3.9: Input data for OHCU, LOB, DHDT, and GT units. 


Refineries Process units 

OHCU LOB DHDT GT 
Haldia (13,13.2,13.5) (12,12.6,14) (7.2,7.5,7.9) (5.5,5.9,6.3) 
Gujarat (1.8,2.1,2.3) (1.12,1.14,1.17) (1.43,1.45,1.48) (2.6,2.9,3.1) 
Panipath (2.1,3.3,3.4) (3.21,3.42,3.46) = (3.42,3.45,4) 
Paradip - (0.31,0.38,0.47) (0.67,0.69,0.72) - 
Mathura (4.3,4.4,4.49) = (4.18,4.24,4.38) (4.67,4.75,4.86) 
Guwahati = (2.15,2.19,2.27) (1.73,1.79,1.83) za 


coefficients of both level objective functions and location, shape, and specified 
probability levels of oil refineries are shown in Tables 3.10 and 3.11, respectively. 
Using Tables 3.8—3.11, model P3 is transformed into the following model: 


Table 3.10: Input data for cost coefficients. 


C11 = (800, 820, 850) 


€12 = (380, 420, 435) 


C13 = (305, 310, 325) 


C14 = (210, 225, 240) 


C15 = (180, 197, 213) 


C16 = (170, 175, 180) 


C17 = (150, 156, 159) 


Cig = (140, 148, 153) 


C21 = (80, 85, 89) 


C22 = (52, 53, 55) 


C23 = (46, 48, 49) 


C24 = (32, 35, 38) 


C25 = (27, 27.6, 28) 


C26— (17; 17.6, 17.9) 


C27 = (12, 12.9, 13) 


Cog = (8, 8.7, 8.9) 


Table 3.11: Input parameters for location, shape, and probability levels. 


Refineries Location Scale Specified 

parameter (6;) parameter (8;) probability levels (y;) 
Haldia 200 20 0.05 
Gujarat 64 10 0.06 
Panipath 19 9 0.07 
Paradip 6 3 0.07 
Mathura 5 3 0.08 
Guwahati 7 2 0.09 
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P11 


maximize Z,(x) = 823.33x, + 411.67x2 + 313.33x3 + 225Xx4 + 196.67X5 + 175X6 + 155x; + 147Xg 
maximize Z5(x) = 84.67x; + 53.33X2 + 47.67X3 + 35X4 + 27.53X5 + 17.5X6 + 12.63x; + 8.53xg 
subjectto 5.83x, + 6.27x2 + 11.43x3 + 11.2%, + 13.23X5 + 12.87X6 + 7.53X; + 5.9Xg x 178.06 

133x,4 2.33x2 + 2.41x3+ 245x44 2.07x5+ 1.14x6+ 1.45x7 + 2.87Xg < 53.66 

2.3X, + 1.4X3 + 2.93x5+ 3.36X6 + 3.62xg< 10.20 

0.65x;  0.89x; 4 0.45X4 + 0.39xg + 0.69x7> 13.87 

5.78X, + 6.63X2 + 6.23x3+ 5.14X4 + 4.4X5 + 4.27X7 + 4.76Xg > 12.45 

245x,4- 2.19x3+ 1.92x4 + 2.2x6+ 1.78x7> 11.72 

xj vj 


Using Algorithms 3.1 and 3.2 of computational algorithm, optimal solutions are ob- 
tained from P11 by LINGO 15.0 iterative scheme and GA and put into Tables 3.12 and 
3.13, respectively. 


Table 3.12: Optimal results for Algorithm 3.1. 


Method ZA Eé Xi X2 X3 Lé Lé Xe X; Xg 


LINGO 11,936.45  1,448.545 4.435 19.574 0.0 064 0.0 0.0 0.537 0.0 


GA 12,141.87  1,448.597 5.568 17.774 0.0 0.58 0.0 0.0 0.71 0.0 


Table 3.13: Optimal results for Algorithm 3.2. 


Method ZA Z5 X X2 X3 Xu X5 Xe X; Xg 


LINGO 3,447.215 304.995 0.0 0.0 00 1.542 0.0 0.322 18.914 0.764 


GA 3,447.324 306.866 0.0 0.0 0.0 1.627 0.0 0.413 18.728 0.722 


3.6 Discussion 


In our chapter, we solved a stochastic TLPP for a noncooperative game under fuzzy 
environment using linear and nonlinear membership functions. For Example 1, from 
Tables 3.3 and 3.4, it has been observed that the optimal solutions have been derived 
using LINGO 15.0 iterative scheme and GA, respectively. In case of linear membership 
function under GA, we have the best solution for A= 0.899 and d = 0.1 corresponding 
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to both level DMs Z; = 7338.571, Z; = 2848.670. However, for the nonlinear membership 
function under GA, the results are Z, = 2718.753, Z = 1037.992 with A = 0.99998, t=0.5 
and d=0.1, respectively. For Example 2, from Tables 3.6 and 3.7, it is concluded that 
the optimal solutions are reached only for the linear membership function under GA 
Z, = 1045.983, Zz = 559.914 with A=0.977 and d=4.4. Again, for the production plan- 
ning problem, from Tables 3.12 and 3.13, we have observed that the optimal results 
have been obtained for both levels Z, = 12141.87, Z; = 1448.597, d = 4.4, À = 0.999 corre- 
sponding to the linear membership function, whereas Z, = 3447.324, Z; = 306.866, d = 
1.1, t= 0.2, and À = 0.9567 corresponding to the nonlinear membership function. So the 
best solution occurred corresponds to the linear membership function under GA. 
Finally, we have concluded that the three numerical experiments that attained the best 
solution correspond to only consider the linear membership function under GA. 


3.7 Concluding remarks and future research 


In this chapter we considered a TLPP for a noncooperative game where the cost pa- 
rameters and the constraint parameters except on the right-hand side are either tri- 
angular or trapezoidal fuzzy numbers, whereas the constraints on the right side 
follow an extreme value distribution. Using Mellin transformation, fuzzy numbers 
are transformed into crisp value. Thereafter, the constraints in probabilistic nature 
have been converted into a crisp form by using the stochastic programming proce- 
dure. So it becomes a nonlinear programming problem. Then based on the compu- 
tational algorithm, reduced problem is solved by LINGO 15.0 iterative scheme and 
GA, respectively and then compared the results. 
The following aspects are treated as new features of our chapter. 
(a) For the first time, we have proposed a probabilistic TLPP in Stackelberg game 
under a fuzzy environment using linear and nonlinear membership functions. 
(b) Based on Mellin transformation, all types of fuzzy numbers (triangular or trape- 
zoidal) are converted into crisp values. 
(c) Probabilistic constraints follow an extreme value distribution to be converted 
into a deterministic form using stochastic programming. 
(d) Comparative study about the results by LINGO 15.0 iterative scheme and GA 
correspond to nonlinear and linear memberships. 


In future studies, this method can be considered in another type of nonlinear mem- 
berships such as hyperbolic, parabolic, and logarithm. This research work can be 
extended for multiobjective TLPP in noncooperative game, including managerial 
decision-making problem and transportation problem. 
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4 Intuitionistic fuzzy trigonometric 
distance and similarity measure and their 
properties 


Abstract: In this chapter, we introduce the concept of intuitionistic fuzzy sets (IFS). 
Similarity and distance measures between IFS are explained and also extended these 
measures to IFSs. Some new trigonometric similarity and distance measures of IFSs 
are studied, and it is shown that intuitionistic fuzzy distance measures satisfy the re- 
quired identities of intuitionistic fuzzy similarity measures. Also, basic properties of 
trigonometric similarity and distance measures of IFS (IF,) are described. 


Keywords: Fuzzy set (Fs), intuitionistic fuzzy set (IFs), distance measures (DM), sim- 
ilarity measures (SM), intuitionistic fuzzy distance measures (IFDM), intuitionistic 
fuzzy similarity measures (IFSM) 


Subject classification codes (2010): 94D05, 94A15, 94A17. 


4.1 Introduction 


The concept of fuzzy set was proposed by Zadeh [1]. Fuzzy sets conjecture has to 
transact business with uncertainty and dreaminess anywhere the notch of belonging- 
ness of any part to a backdrop which is represented by a attachment function. 
However, in the frank mean condition, it is not constantly doable that the scale of 
nonmembership meeting of a facet to the fuzzy set is barely peer to 1 minus the level 
of membership. It implies that there is some kind of degree of ambiguity and has 
been called hesitation degree. 

The classical concept of fuzzy sets [1] is a strong tool to deal with uncertainty, 
vagueness, and imprecise information. There are various generalizations of higher 
order fuzzy sets for different specialized purposes. Intuitionistic fuzzy sets (IFS) in- 
troduced by Atanassov [2] are quite useful and applicable, and fuzzy sets are de- 
scribed by two functions: a membership function and a nonmembership function. 
Atanassov [2] introduced the concept of IFS as simplification of fuzzy sets, which is 
competent of capturing the order that includes several grades of timid and has been 
practical in countless return of research. 
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Thus, another degree of freedom has been introduced to describe IFSs. Such a 
generalization of fuzzy sets equips us with an additional possibility to represent in- 
complete knowledge, which leads to describe many real problems in a more adequate 
way. Further, Attanssov and Gargov [3] also introduced the concept of interval-valued 
IFS and studied their properties. 

Atanassov [4] provided an example, where fuzzy sets are IFS; however, con- 
verse is not automatically rightful. Burillo and Bustince [5] suggested an algorithm 
for IFS from fuzzy sets with the help of Atanassov’s [6] operator. Each of the issues 
that manage fuzzy set hypothesis can also manage IFS hypothesis; however, there 
are circumstances where fuzzy set hypothesis cannot be appropriate to make a rea- 
sonable investigation. The membership values cannot be always determined, due 
to inefficiency of available information, besides the presence of vagueness in the 
information. Similarly, while determining the nonmembership, the same situation 
arises. A part of such estimation remains of indeterministic nature. This indetermin- 
istic part has not been considered in the case of fuzzy set theory, so we have IFS. 
For instance, suppose that a human being is to express his/her opinion about the 
beauty of a flower. In classical theory, he/she has only two choices 0 or 1, that is, a 
flower is either beautiful or ugly. In fuzzy sets, he/she has the liberty to express 
his/her view that the flower is beautiful with a grade of membership 0.8, which de- 
termines that the grade of nonmembership is 0.2, but in real life, a flower may have 
the beauty with a grade of membership 0.8 but not ugly at all, that is, its grade of 
nonmembership has the value 0. Therefore, one needs a most appropriate ap- 
proach, which would help to deal with such a problem, whose situation depends 
on human perception. Therefore, IFS is a very useful mathematical tool in many sit- 
uations having imprecise, vague, and inexact data. 

It is noted that IFSs can describe the fuzzy characters of the things more compre- 
hensively and thus are powerful and effective tools in dealing with fuzzy information. 

An intuitionistic fuzzy conventional is characterized by two functions. Former 
is the point of link and subsequently is the level of nonmembership function. Their 
addition is each time fewer than one and the same to one. 


Definition 4.1 (Fuzzy Set). 
Fuzzy set E in a universal set Y is defined as 


E={<y,ag(y)>:y € E} (4.1) 
where a@g(y;) is a membership function defined as follows: 


0, y does not belong to E and there is no ambiguity 
üg(y) za 1, y belongs to Eand there is no ambiguity 


0.5, there is maximum ambiguity whether y belongs to E or not 
(4.2) 
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In fact ag(y) associates with each y, a grade of membership, to the set E. 


Definition 4.2 (Intuitionistic Fuzzy Set). 
An IFS E in a universal set Y is defined as 


E- {<y,ag(y),Bg(y) >: y € E) (4.3) 
where the functions 


E- («y.ag(y).Bg(y) >: y € E) (4.4) 


are the degree of membership to the element and the degree of nonmembership to 
the element, respectively, with the condition 


0 xaz(y) * Be(y) <1 (4.5) 
the intuitionistic index of an element y € Y is defined as follows: 
eil =1-az(y) Bei) (4.6) 


mal is called a degree of membership of an element in E and it is evident that 
Ox rte(y) x 1. 
Note that for an IFS E, if ag(y)=0 then ag(y)+f8,(y) and if ag(y) -1 then 


7tg(y) =0. 
We present the following basic operations on IFS, which is necessary for our 
discussion: 


Definition 4.3 (Basic Operations of Intuitionistic Fuzzy Set). If E and F are two IFS 
of the universal set Y, then 

i) iff E C Fthenag(y) <ar(y) and gly) > Beil Vy € Y 

ii) iff E- F then ag(y) 2 ar(y) and Bei) - Bet) 

iii) EC = {< y, Bi (y), az(y) >: y € Y}, where E% is complement of E 


iv) En F- (xy, min(ag(y),ar(y)), max(ag(y),ar(y)) >: y € Y} 


v) E UF= {< y, max(az(y), ar(y)), min(az(y),az(y)) >: y € Y] 


Burillo and Bustince [5] have also given some results on IFS. In the subsequent dis- 
cussion, the following notations are used: 

R=[0, 2%), Y= {y1, Y2 .. yn) is the universal set; FOI is the class of all fuzzy 
sets of Y; IFS(Y) is the class of all IFSs of Y. 

Liu [7] given new measures without saying definitions of aloofness actions and 
similarity events of fuzzy sets and examined the connections between them. These 
axioms are spontaneous and possess laboring in the fuzzy literature. 
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Definition 4.4 (Fuzzy Distance Measure). Let D be a real-valued function defined as 
D:F, xF; — R*. A distance measure D is called fuzzy distance measure if it satisfies 
the following conditions: 


(d1) D(E,F)=D(F,E) VE,F € F, 

(d2) D(E,E)=0 VE € Fs 

(d3) Did, d^) = maxz rer, D(E, F), if dis the crisp set 

(d4) VE,F,G € F,, if E C F C G then D(E, F) < D(E, G) and D(E, F) < D(F, G) 


Definition 4.5 (Fuzzy Similarity Measure). Suppose that S be a real-valued function 
T:F,xF,— R*. A similarity measure T is called fuzzy similarity measure if it satisfies 
the following conditions: 


S1) T(E,F) - T(F,E) VE,F € Fs 
S2)T 
S3 
S4) VE,F,G € Fs, ifEC F C Gthen T(E, F) > T(E, G) and T(E, F) > T(F, G) 


( 
(S2) T(d, d^) =0, if dis the crisp set 

(S3) T(K, K) = maxz zer, D(E,F) VK F 

( 

If we standardize d and T in Definitions 4.3 and 4.4, we get Ox D(E,F) x1 and 
Ox T(E, F) «1 for E, F € F,. The relation between standardized d and S is D=1-T. 
Hence, from this relation, it can be seen that the distance and similarity measures 
are dual concepts, and we can give the distance and similarity measure axiomatic 
definitions for IFS. 


Definition 4.6 (Intuitionistic Fuzzy Distance Measure). Suppose D is a real-valued 
function D:IFS x IFS — R*. Then distance measure D is called intuitionistic fuzzy 
distance measure if it satisfies the following conditions: 


(ISI) T(E, F)=T(F,E) VE,F € IF, 

(IS2) T(d, d^) =0, if is the crisp set 

(IS3) T(K, K) = maxz rer, D(E,F) VK € IF, 

(IS4) VE, F, G € IF,, ifE C F C G then T(E, F) 2 T(E, G) and T(E, F) > T(F, G) 


Definition 4.7 (Intuitionistic Fuzzy Similarity Measure). Let T be a real-valued 
function T:IF, x IF; — RL Then similarity measure T is called fuzzy similarity mea- 
sure if it satisfies the following conditions: 
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(IS1) T(E, F) = T(F,E) VE,F € IES 

(IS2) T(d, d^) =0, if dis the crisp set 

(IS3) T(K, K) = mass rar, DIE PI VK € IFs 

(IS4) VE, F,G € IFs, ifE C F c Gthen T(E, F) > T(E, G) and T (E, F) > T(F, G) 


4.2 Sine trigonometric intuitionistic fuzzy 
distance measures 


Definition 4.8. Let E and F be two IFS on universal set Y as Y= {y1,y2, .... Yn}, 
whose membership values ag(yi), i=1,2, ...,nandar(yi), i=1,2,...,n, respectively, 
are the measures of sine intuitionistic fuzzy distance Dir, (E, F) between Dyr, (E, F) fuzzy 
sets E and F, as 


Dire, (E, F) = 53 (sin (hl Soe sin (Nd 5091) (4.7) 


i=1 2 2 


Dal BF) = AE nl DER VEDAL) ns sin (Le vB) 
(4.8) 


Theorem 4.1. Identities (4.7) and (4.8) prove all the conditions of intuitionistic fuzzy 
distance measures. 


Proof. Given that 


Here, we prove all four conditions of intuitionistic fuzzy distance measures that are 
mentioned in Definition 4.6: 


(Id1) Dig, (E, F) = x X E (Eran Sen) m+ sin (Ern fett) d 


iz1 2 2 


M 27 E S Je sin (Pro? 2 d 
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(Id2) Dir, (E, E) = = S E (Econ zer) m+ sin (Ego Bi d 


-0 VE, c IFs 


(Id3) Let 

E - d and F - d^, Le, ag(yi) = aa(yi), ar (yi) = age (yi) and Bg(yi) = Bali) 5 
Bei) = Bac (Vi) 

Then we get 


De, (d, d^) -1 VE,F € IFs 


that is 
Dir, (d, d^) = maxz, rer, Dig; (EF) VE,F € IFs. 
(Id4) Forall E,F, G € IFs(Y), ifE C FCG, thenaz(yi) <ar(yi) x ac(yi) 
or Br (Yi) 2 Be(yi) 2 Boy) 


(En un Jn (Eoo SA Jes e SE )n 


p (Bebo Aon), 


and the nature of sine function, we get 


(till, sin (M) OL). 
2 2 


ac (E? Pal) | )n 


> Dr, (E, F) < Dir, (E, G) 


sin 


(E SE | )n 


Similarly, 


(E - ar(yi)] Je (en -Br(yi) es Gua - ag(yi)] Jn 


2 2 2 


EEN 


4 Intuitionistic fuzzy trigonometric distance and similarity measure —— 59 


sin (E99 cee) m+ sin (UN Bi m< sin (Seb) ces) T 
2 2 2 
+ sin joe T 
> Dr, (E, F) < Dir; (F, G) 
Now we will prove all the required conditions using identity given in eq. (4.8): 


i=1 


~ 2n4 


Léi E ( deht Vay] Je ge ( Gran : VA Ja 


= Dis, (F.E) 


(Id2) Di, (EE)- LY; z : GE væ) Je ol : EU V Beo) Ja 
=0 VEEIFs 

(Id3) Let 

E-dandF- d^, i.e., ag(yi) =aa(yi), ar (yi) = Gar (yi) and Bp (Vi) » Ba(yi); Br (yi) = Bac (Vi) 

Then we get, 


Di, (dd) -1 VE,F € IFs 


that is, Dig, (d, d^) = Dar, FelFs Dr, (E, F) VE,F € IFs. 


(Id4) For all E,F, G € IFs(Y), ifE C F C G, then ag(y;) xar(yi) xac(yi) or 
Br (Vi) 2 Be (yi) 2 Bg (Vi) 


Ga - Vary] y (Bn fei) y (Sp - V'ac(yi) ) 
2 2 R H 


" (2w SEH ) 
2 


and the nature of sine function must be defined as 
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sa (Ier rbl y ol VB EN) ) 


< sin 


(ie - Vasy] y sa (Iro - VB On) y 
2 2 


=> Di, (E, F) € Dra (E, G) 
Similarly, 


(m - vVar] y 
2 


felt) - v/Beon)| Jrs (xem - V/ac() ) 
2 B 2 


+ 


WEIN APO 


IVBO) - VBW] ) 
2 


nl NS as) in nl VT AEN) ) 


Ne 


< sin m+ sin 


(aw? - Vac] [^A - VBson)| ) 
2 2 f 


Se et 


> Dig, (E, F) < Dir, (F, G) 


4.3 Properties of sine trigonometric intuitionistic 
fuzzy distance measures 


Let D be a distance measure on IFS and for E,F € IF,. If they satisfy that for all 
yi € Y, either E(y;) C F(yi) or E(yi) 2 F(yi), then 


(i) Dig, (EUF, En F) = Dig, (E, F) (4.9) 

and 
Dg, (EUF,EnF)-Dig, (EF) VE,F,G € IFs (4.10) 
(ii) Dig, (E UF, G) < Di, (E, G) + Digg, (F, G) (4.11) 

and 


Dif (E UF, G) < Dir,, (E, G) Ke Dir,, (F, G) VE,F,G € IFs (4.12) 
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(iii) Dir (E nF, G) < Dir, (E, G) + Dira (F, G) (4.13) 


and Dig, (E nF, G) < Dig, (E, G) + Dip? (F, G) VE,F,G € IFs (4.14) 
then the above equations satisfy the properties. 


Proof. 
(i) Using Definition 4.6, we get 


1X i) - i 
Dig, (E UF, EnF) = oD [ein (eee sea )n 


i=1 


+ sin 


gu E [sm (ero eil je sin (20 Beo) Ja 


is sm (ert ml je sn (20 = Bly y] 


Yi CH: 


(Sp Bene Jz 


= Dire, (E, F) 


Similarly, we can prove for eq. (4.10). 
(ii) Let us consider for 


Dir, (E, G) $ Dir, (F, G) Se Dir, (E UF, G) 


f =a E = aa) Ier inf BE yi) - B«(Y?)| ja 
" x bis (0 2:500 Jas ol EH) At. yi) - Bet) d 


- LY bel Seite te), sin (Bevel) Bull) y] 


= L55 | (ee 9697), in (0 = Bol) 


(ler Seul ) (Eoo zd ) d 


1 n 


+ sin m+ sin 
2n n 
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SÉ E sin (10 75000), Bo) = Bol J 
á 2 nl D yi) SE | Je ol Vi ra )n| 
af bs al" d zept ye sn Sr y) Aan! y 
n ve 
n F M [s (Eon M Ion sin (Pro? ZEN d 2 
mech 


> Dir, (E, G) + Dir, (F, G) - Dir (E UF, G) 20 


Similarly, we can prove 


Dr, E,G)+ Dr, F,G)- Dr, (E UF, G) >0 
(iii) 
Dir, E,G)+ Dir, FG - Dig (E nF, G) >0 


and 


Dig, E,G)+ Dig, FG - Dig, (E nF, G)20 VE,F,G € IFs 


4.4 Cosine trigonometric intuitionistic fuzzy 
similarity measures 


Definition 4.9. Let E and F be two IFS on universal set Y as Y = {y1, y2, .... ya]. 
whose membership values are ag(yi), i2 1,2, ...,n and ag(yi), i=1,2, ...,n, respec- 
tively. Then the measures of cosine intuitionistic fuzzy similarity Sig, (E, F) between 
Sig, (E, F) fuzzy sets E and F are defined as 


Sir, (E, F) = » 5 feos (See wL) T+ COS (NH Beer) d (4.15) 


i=1 2 2 


WE 


E 


i=1 


(4.16) 


4 Intuitionistic fuzzy trigonometric distance and similarity measure —— 63 


Theorem 4.2. Identities (4.15) and (4.16) prove all the conditions of intuitionistic fuzzy 
similarity measures. 


Proof. Given 


Si, (E, F) = > : feos (eman T+ COS (fm tees) d 


i=1 


Here we prove all four conditions of intuitionistic fuzzy distance measures that are 
mentioned in Definition 4.7: 


(IS1) Si, (E, F) = aps cos (Ero ee 7-4 COS (Bos cnn d 


EIS STEE 


(IS2) S(d, d^) =0, if dis the crisp set 
(IS3) S(K, K) = MaXg, re]Fs D(E, F) VK € IFs 


(IS4) VE, F, G € IFs, ifEC F C G then ag(y;) x ar(yi) «ag(xi) and 


pon shes pie 2 (E pope zs 


Aldi pit, «(Eo hon, 


and the nature of cosine function is given by 


ape ax(yi) 2609) IN a cos (Ie yi) Patil d 
«x Ze! az (Vi) ECN) cos (FeO) - Br). 


Sing, (E, F) 2 Sig; (E, G) 


Similarly, 


2n 2 2 


> 5, | (Mees c(yi) Pss „(Ew FDN 


a(o (Wat = Bob) 
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and the nature of cosine function is given by 


Rol? ee) - ao) secos (I y) AD, 


A Rol ar(y:) ODN on (BeBe) 


Sigg, (E, F) E Sigg, (F, G) 


Similarly, we can prove the identity given by eq. (4.16). 


4.5 Properties of cosine trigonometric intuitionistic 
fuzzy similarity measures 


Let T be a similarity measure on IFS and for E,F c IFs. If they satisfy that for all 
yi € Y, either E(y;) C F(y;) or E(yi) 2 F(y;), then 


(i) Tir, (E UF, ENF) = Trg, (E, F) (4.17) 

and 
Trg, (E UF,EnF) - T, (E,F) VE,F,G € IFs (4.18) 

and 
(ii) Tir, (E UF, G) < Tr (E, G) + Tif (F, G) (4.19) 
Tis, (E UF, G) < Tr, (E, G) t Ti, (F, G) VE,F, Gc IFs (4.20) 
(iii) Trg, (E nF, G) < Tr (E, G) + Tr (F, G) (4.21) 

and 
Ti, (E nF, G) < Tir; (E, G) + Ti, (F, G) VEFGc IFs (4.22) 

Proof. 


(i) Using Definition 4.7 


n . — . 
Ti, (EU F,EnF)- = E (Een aenF(¥i)| )n 
1 


2 


i- 


He (P 


- Bea rQYi)| ) 


2 
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as bs E Ge = ast) ID e (Eco Pro) Jz 


ch 


n > E [tm eem. ies) mt Jr cob (Em Bun JW 


A= 6) 


"x [cos (Met Prt yi) SAM Jn cos (Aa = j| 


= Trig, (E UF,EnF 


Similarly we can prove eq. (4.18). 
(ii) Let us consider for 


Tir (E,G)* Ties, (F,G)- Tir, (EUF,G) 
i E J feos (=w) T+ COS [En d 


" I5 [cos (Ero? 2d Jr: cos (Pe SC )n| 


i=1 


ioe ? [cos (Mev - ac(yi)] Je: cos( Best -Ao d 


H 


SE feos (20 -ac(y)| Jas cos (Bet Bol) VW 


ke 2 


" Se feos (2% e Jr: cos (Be SC )n| 


i » cos (#20 2500) ID cos (He eil y 


5 feos (Me SE Jn cos (Bet? %1) d 


y1€Y2 


Similarly, we can prove Tig; (EU F,G) < Tg; (E,F) + Tr (F, G). VE,F,G € IFs 
(iii) Similarly, we can prove 


Tr (E nF, G) < Tir; (E, G) + Tr (F, G) 
and 


Tr, (E nF, G) < Tr, (E, G) * Tr, (F, G) VE,F,G € IFs 
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4.6 Conclusion 


First, we defined IFS. Then, we introduced a general intuitionistic inclusion on a uni- 
verse of discourse X and intuitionistic similarity measure between two IFS on X. 
Definitions have been given to illustrate the concept of IFS, intuitionistic similarity 
between two IFS, and intuitionistic fuzzy relation based on intuitionistic similarity 
between two IFS, intuitionistic fuzzy distance measure, sine and cosine trigonometric 
intuitionistic fuzzy distance measures. Properties of sine and cosine trigonometric in- 
tuitionistic fuzzy distance measures are also given. 
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Alok Dhaundiyal, Suraj B. Singh 

5 Distributed Activation Energy Modeling by 
the transmutation of different density 
functions 


Abstract: This chapter incorporates the different perspective of modeling of activa- 
tion energies through ingression of an additional parameter of distribution function 
to increase the flexibility and controlling ability of modeling by the linear mixing. 
The transmuted functions are adopted to bring flexibility in the numerical solution 
of multireaction models. A stochastic model is used to represent the activation en- 
ergy function. Thus, the different activation energies can be represented by a density 
function to describe the pyrolysis process. The qualitative, as well as quantitative ef- 
fects of the transmuted family of distribution on biomass pyrolysis, are carefully ex- 
amined by the proposed mathematical solution. The nonlinear thermal history is 
implemented to demonstrate the approach. 


Keywords: Biomass pyrolysis, transmuted function, kinetic parameters, thermal pa- 
rameters, nonlinear ramping 


Notation 

SES Transmuted form Dimensionless 

W Weight parameter for volatile content Dimensionless 

T Temperature K 

a Conversion Dimensionless 

A Frequency factor min“? 

E Activation energy kJ/mol 

p(E) Density function of activation energy Dimensionless 

T Shape parameter Dimensionless 

R Universal gas constant kJ/K mol 

0 Shape parameter for the Lindley Dimensionless 
distribution 

y Shape parameter for the Lindley Dimensionless 


distribution 
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y Location parameter for the Lindley 
distribution 

B Shape parameter for the Weibull 
distribution 

n Scale parameter for the Weibull 
distribution 

y Location parameter of the Weibull 
distribution 

p Shape parameter of the Frechet 
distribution 

ð Scale parameter of the Frechet distribution 

Hu Location parameter of the Frechet 
distribution 

Eo Mean value of the Gaussian distribution 

oO Standard deviation of the Gaussian 
distribution 

A Location parameter of the Rayleigh 
distribution 

X Scale parameter of the Rayleigh 
distribution 

E Stochastic nature of regression model 
(error) 
Rate of change of heating rate 

E Instantaneous value of time interval 

m Heating rate 

To Initial temperature 

T The maximum temperature 

Es Central value of activation energy 

Ew Step width 

W(k) Lambert W function 

K Time rescaling factor 

yt Maxima point 

ys Energy rescaling corresponding central 
value 

Yw Energy rescaling factor for step width 

[^] Dimensionless number 

Ez; Upper bound of activation energy 

fiand fo The components of the Lindley 
distribution 
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Modeling of a system is one of the most challenging tasks when you are dealing 
with the unit where multivariables are acting simultaneously within a system. One 
of such systems is a biological activity of constituents that have different chemical, 
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thermal, and physical properties, and it emerges out to be a reason of anomalies in 
the measurands. There are numerous ways of tackling such complicated problems, 
which need myriads of steps to evaluate them properly. Therefore, some assump- 
tions, approximations, and advanced calculus are required to save three M’s (men, 
machine, and money) from overexploitation. 

Here, one of the thermochemical processes, namely pyrolysis, has been ana- 
lyzed qualitatively and quantitatively of a given problem. The problem arises when 
the dynamical behavior of inflexion point of thermogravimetric (TG) variation of 
mass curve impedes the predicted stochastic modeling of activation energies to con- 
verge around the neighboring point of it. Many methodologies [1-7] have been 
adopted to fetch the best suitable prediction that can converge up to 90% of the 
measurand value. 

Numerical solution of multireaction modeling of parallel reactions provides a 
better insight of parametric behavior and its effect on the natural phenomenon; 
thereafter, it becomes a lot easier to predict the outcome of solution than that of 
statistical linear models that merely rely on the linear regression of a straight line. 
There are various methods of modeling and evaluating kinetic parameters of pyrol- 
ysis process that can also be categorized as the single reaction model and the multi- 
reaction model [8-12]. 

The demarcation of modeling is based on the end results required by the user; 
therefore, it is necessary to understand the requirement of the given problem. In 
the same context, exploratory and application-based models are of the main con- 
cern [13]. The exploratory models require less computation resources and it is rela- 
tively fast and flexible, but the accuracy is highly compromised. On the other 
hand, the application-based models look for the concreate, robust, mature, and 
precise solution. It requires a little computer-based programming, but solution 
converges without reservation. In the real-time problems, the application-based 
modeling provides an accurate and sophisticated numerical solution which can 
be delineated graphically. 

Hence, in this study the effect of one of such application-based models is con- 
sidered to carry out the simulation process. Thermochemical parameters along with 
distribution pattern are examined through modeling of activation energies by the 
transmuted functions. Moreover, the dynamic behavior of stage decomposition is 
depicted through nonstationary maxima of curve with respect to time. The purpose 
of study is solely based on the application perspective of transmuted function, re- 
ducing the variation gap between Arrhenius parameters, and proposing the new 
concept for modeling activation energies, not to over-fitting of curves. Therefore, a 
new scheme has been proposed where the kinetic parameters are not influenced by 
the inflexion point of two different reaction phases, which is overlooked in the pre- 
vious studies [1, 3, 14]. 
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5.2 The transmuted function of activation 
energy p(E) 


The idea of grouping the parallel reactions to analyze the mechanism of complex re- 
actions is not the new scheme of comprehending the interrelationship between fre- 
quency factor and activation energy. For surface-catalyzed reactions, it had been 
envisaged that the surface has various sites that rely on distribution of activation en- 
ergies among the competitive reactions [15]. The similar concept was adopted to mea- 
sure irreversible changes in an electrical resistance while annealing the evaporated 
metal films [16]. This scheme was highly adopted by the coal pyrolysis community 
and the parallel reactions approach came into light in the 1960s [17, 18]. Later, 
Anthony and Howard (1976) proposed that the activation energy of the parallel reac- 
tions could be expressed as a density function of the normal distribution for coal py- 
rolysis [19]. However, there are some asymmetrical functions tested to validate the 
result for the loose biomass [20-22]. It is not that the derived results from other distri- 
bution functions have higher degree of entropy than the normal distribution. But the 
compensatory effect due to the mathematical mandate in the old theories that get 
drifted off through modification of the prevailing theorem [1]. 

Despite all those efforts for calculating multivariate function [1, 23], still there is a 
scope of assessing the nature of stochastic function as well as selection of appropriate 
thermal programming for pyrolysis reactions. Although it is difficult to interpret TG 
and derivative TG curves exactly, since (i) it is impossible to measure the sample tem- 
perature exactly through thermocouple, (ii) there is a time lag while formation of in- 
termediate stages, (iii) interference of interactive atmosphere or cooling effect, and 
(iv) heat transfer limitation. Experimental curves have constraints and stagnation 
with time; overall, it reflects the gradual variation in property of material, which is 
related to alteration of thermal history. On the other hand, the attribute of mathemat- 
ical model is based on the set protocol that overlooks the miscellaneous factors re- 
lated to thermal decomposition. Transformation time between intermediate stage is 
almost At, which tends to zero for the mathematical model. 

There are some methodologies developed to simulate the TG attribute effec- 
tively such as double-distributed activation energy model (2-DAEM) [5, 14]. The 
scheme manifests to mitigate uncertainties of distribution function to describe 
TG curves by introducing the weight parameter for volatile content *w" [24]. This 
proposition can be valid unless pyrolysis of biomass is carrying out in the unreal- 
istic condition since it is not necessary that the role of volatile content ceases to 
act with the decreasing value of w. Lewellen et al. [25] stated that the residence 
time of volatile content has significant influence on pyrolysis of biomass. It 
decides the extent of secondary pyrolysis. Therefore, the weight factor and time 
factor cannot be correlated, as competitive pathways such as cracking, repolyme- 
rization can be inhibited through formation of some of light volatile gases. The 
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volatile content may initiate the autocatalytic reactions at elevated temperature. It ac- 
celerates the formation of the stable products content which would escape the bio- 
mass matrix to yield volatiles [25]. There are various key issues related to modeling of 
biomass that have been summarized in some of literatures [25-28]. It is complicating 
to find the exact mathematical solution to interpret the TG curves, but the discor- 
dance can be subdued to some extent. There are some other mathematical solutions 
used for condensed phase reactions — isoconversional and model fitting. Both are 
highly recommended for kinetic analysis [13, 29]. 

Regression-based modeling is very promising and time-saving, whether the esti- 
mated results are precise and exhibit the kinetic mechanism of pyrolysis. This must 
be investigated carefully. Isoconversional and model fitting are some of regression- 
based methods used for estimating the kinetic parameters. Isoconversional scheme 
or model-free kinetics does not imbibe the idea of fitting the experimental TG data to 
the phenomenological models, but it does not infer "free from assumption." The 
basic tenets of scheme are that at any given extent of reaction, the same reactions 
occur in the same ratio at different heating rates and the reactions must be indepen- 
dent of temperature. However, this hypothesis can be easily nullified by a system of 
parallel independent reactions whose relative reactivity is proportional to tempera- 
ture and it is varying due to different activation energies [30]. It is also possible that 
the overall reaction pathway of the competitive reactions with different activation en- 
ergies is different at different temperature [31], as Lewellen et al. [25] suggested for 
pyrolysis of cellulose. Unlike isoconversional method, inconsistency for determining 
Arrhenius parameters is observed in the model-fitting method. It is mainly due to 
force-fitting of date to the hypothetical reaction model. The Arrhenius parameters are 
determined by the phenomenological model, which is already assumed. Moreover, 
both temperature and conversion consistently vary with time. This methodology is 
unable to demarcate separately the temperature dependence of rate constant and the 
conversion. Eventually, it causes drastic variation while determining the frequency 
factor and activation energies, and any phenomenological model can be claimed to 
describe the pyrolysis mechanism of biomass. Moreover, the existence of competitive 
reactions with different activation energies affects the overall decomposition rate, as 
these reactions are highly influenced by temperature and extent of conversion; there- 
fore, the activation calculated from this method is a function of T and a [29]. This 
model is highly inaccurate and has been highly discouraged for kinetic analysis [32]. 

However, instrumental noise, response time, inabilities of thermocouples, and 
sensitivity of thermobalance cannot be incorporated in any advance integral model 
or regression schemes. One can filter data to remove suspicious values and reduce 
the gap between A and E variation. One may also adopt statistical testing scheme 
(Fisher's exact test, Dixon Q's test, Bartlett's test, etc.). But if your thermoanalytical 
data size is colossal, these probing schemes become time-consuming to analyze the 
sample size. Distribution in activation energies through conventional logistical dis- 
tribution function [33] is not only the way to introduce a reactivity distribution, but 
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there are several other ways to solve DAEM, such as linear mixing of density func- 
tion p(E) of activation energy: 


p'(E) - (1 1)p(E) -T EI, |r| <1 (5.1) 


where E is the activation energy. 

Unknown activation energy of pyrolysis reaction is said to be transmuted if it 
follows eq. (5.1). However, the same expression holds valid for cumulative distribu- 
tion function of activation energy. Parametrization of well-established distribution 
reduces the effective computation time for deriving more robust and flexible fami- 
lies of distribution. Shaw and Buckle [34] were the leading researchers to explore 
distributional flexibility through introducing a new parameter to an existing distri- 
bution. They figured out that involvement of parameters did not make solution 
cumbersome but enhance the reproducibility of the modeling of natural process. 
The generated or transmuted family imbibes the parent distribution and impart more 
dynamical attribute to simulated results for various types of data than the conven- 
tional logistical approach. Implementation of transmuted family is very vast, and it 
encompasses all the natural phenomenon that demands the high degree of flexibility 
as the entropy of such systems is very high. One of such study is carried out for model- 
ing the metrological data. The transmuted Gumbel distribution is found to be very 
effective for analyzing the snowfall data of a region [35]. The competitiveness of trans- 
muted exponentiated Frechet distribution over exponentiated Frechet distribution for 
modeling wind speed data is performed to gather the statistical properties for a wind 
energy conversion system. The application indicates that the transmuted model's per- 
formance is relatively good to the parent model [36]. Another field of application of 
transmuted function is the mortality or failure of some system with respect to 
time. It is found that the reliable behavior of transmuted distribution functions is 
very useful for modeling reliability data in the manufacturing industries [37]. The 
application as well as modification of the conventional distributions is reported in 
the several literatures [38-40]. 

In the subsequent sections, the approximation techniques are adopted to sim- 
plify the classical DAEM and the conventional density variable from transmuted 
family of distribution, so that the distribution of reactivity can be introduced 
through the proposed scheme. However, interdependence of frequency factor (A) 
and E will still prevail in modeling of DAEM, but it is assumed that frequency fac- 
tors (A) for all the competitive reactions are invariant and assumed to be constant. 
DAEM is represented through double integral equations for first-order reaction 


(eq. (5.2) [1]: 
ceo t 
tzaz [| (- [afa Jore (5.2) 
0 0 
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Density and cumulative functions of activation energy are illustrated in Table 5.1. 


Table 5.1: Distribution functions and their corresponding parameters. 


Distribution Density function p(E) Cumulative function G(E) £»0 
Parameters 
Lindley [41] d (1+ J(E - y))e-9€-» f1- (1+ eE P) e-ae-n\ dAn 
B B 
p-1) (E -[(£x 
Weibull [42] den e ( 7 ) 1-e ( 7 ) B: n, y 
n\n 
-(1+ E-yy-P E-p\-P 
Frechet [43] s6 Po- (E) je) \ P, õ H 
EN 
d 1 = Re 1 E-Eo 
Gaussian [44] et {3 (1+ SEH Eo, 0 
Rayleigh [45] 2x(E- d)exp{ -X(E- ar) 1- exp{ -X(E- oy} dy 


5.2.1 Distributed activation energy modeling 


A complex material such as biomass or fossil fuel comprises several constituents 
whose thermal decomposition varies greatly with time and temperature. The differ- 
ent products are generated at different stages of decomposition. If the parallel reac- 
tions follow the complex pathways, it becomes more complicating to model such 
process, so it is possible that a simple mathematical distribution will not have flexi- 
bility to describe the reaction profile of some reactions. Copula-based modeling [1, 
23] is recently adopted to illustrate the complexity of decomposition reactions and 
to some extent it is found to be very helpful; however, the validation of results for 
different thermal histories is not done yet. Discrete activation energy distributions 
schemes are developed to provide flexibility to the reaction model. One of such 
models was developed for petroleum generation, which was first published by 
Ungerer and Pelet [46]. This scheme is derived from the isoconversional methods, 
which is mainly designed for estimating the unique values of Arrhenius parameters, 
but later it was used to correlate activation energy variation with conversion. The 
detailed study of development of several other models is reported in the literature 
[47]. But it requires special attention for energy spacing between reaction channels, 
which may affect the validity of the model. Moreover, it is purely based on the non- 
linear regression and the numerical integration, which itself is stochastic in nature. 
A regression analysis is based on the relationship between study and explanatory 
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variables through the series of coefficients that decide the characteristics and role 
of independent variables in any event: 


y=f(X, En, Er dy Ay ...,Ax) +E (5.3) 


where € represents the stochastic nature of regression model that can alter the solu- 
tion if activation energy is not discretized evenly. However, this scheme is highly en- 
couraged for estimating the kinetic properties of kerogen pyrolysis rates [13, 48]. A 
major drawback in discrete modeling is the derivation of the structure, that is, of the 
terms in the discrete-activation energy model. When the form of differential equation 
generating data is identified, one can estimate possible terms for the discrete model 
from numerical methods, but not necessarily the associated coefficients. A good 
structure design De, the form of f(Xi, X», ... , Xi, Gi, Go, ...ax)) gets the kinetic 
model-fitted data properly and provide the good predictive results. Redundant terms 
in a model, on the one hand, can lead to over fitting, whereas too few terms make 
model ineffective. There are some practical criteria considered for effective modeling. 
First, the model should be numerically stable. Second, the selected model makes the 
best predictions, which can be performed by testing the given data for different heat- 
ing rates or thermal profile. Asymptotic behavior of model should comply with those 
of the data. Keep the demarcation of discrete and continuous models for abeyance 
and formulate the transmuted function-based scheme. 

To demonstrate the scheme, eq. (5.2) is bifurcated for simplifying the mathemati- 
cal expression of DAEM. It comprises two parts: double exponential term (DExp) and 
the continuous density function of activation energy. Laplace integral method is im- 
plemented to estimate the approximated form of DAEM. DExp terms are an implicit 
function of time and temperature history experienced by the sample. The second part 
is independent of time and depends on the distribution of reactivity. First, the behav- 
ior of implicit function is perused, and then approximations are derived for the non- 
linear profile. 

The ideal thermal histories are seldom implemented for the industrial applica- 
tions; therefore, the computational algorithm is tested for the nonlinear thermal 
profile represented as follows (eq. (5.4)): 


T(K) = as? + To(K) (5.4) 


where EN = am and the value of “a” represents the rate of change of heating 


rate (m) with time. 


Let, f(T) = - & 


Tf 
A 
DExpecexpe: | = | ——— etra) (5.5) 
To 2a E 
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E 


The ratio of Na 


) — co is assumed to be large; thus, the significant contribution 


of integral is at neighborhood of T - T; and f(T;) » 0. This provides the asymptotic 
approximation to the temperature-dependent function: 


" A X -exp(f(T)) 
DExp~exp (; Jatt; T F (Tj) (5.6) 


Correlating the result with the central value and the step size of activation energies, 
eq. (5.6) can be rewritten as 


DExp~expf -exp (= 3 } (5.7) 


w 


Rather than arbitrarily spacing the activation energies, it is better to generate the form 
of function that relies on the explanatory variables of experiments, as the trend of acti- 
vation energy variation with respect to conversion is highly indeterministic [29]. 


Let, g(E) = KE 


Asymptotic approximation results in a rapidly varying function as AT — oo 
0 


and function varies from zero to one as E increases by Ew around E,. To further sim- 


plify the function, KE is to be expanded by the Taylor series about E = E: 


g(E)~g(Es) + (E-E,)g (Es) tL (5.8) 


The initial condition for function g(E) at E, is defined as 


g(E;) - 0 and g (Es) xd ie 
Equating eqs. (5.6) and (5.7), we have 
g(E) = —f(T;) + In > (5.9) 
i É 2f (Ts) a(T; - To) i 


After solving eqs. (5.8) and (5.9), we have 


E; =RT;W(x) 
z RE; (T;) 
“Bs + RT; 
where both E, and Ey are implicit functions of time. W is a Lambert W function and 
the dimensionless group m is denoted by x. 


Note: Lambert W function is defined as one of the real roots of equation of form 
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Approximations of Lambert W function for small and large values of “x” (correlated 
to short and long durations of time) can be represented as 


5.2.2 Characteristic of distribution function 


It is observed from the approximated form of DExp, correlated with activation 
through energy spacing of step size Ey, that the function behaves like a smoothed 
See 
24/ (Tr - To) 


Or (& &) — co and the peak is concentrated at the neighborhood of T= Tr. The 


step function. It rises swiftly from zero to one as the term ( ) increases, 


DExp as implicit function E makes the whole integrand homogenous and it be- 
comes easier to compute the complicated integral problem. But the fact that cannot 
be overshadowed is the relative significance of E,, to the standard deviation of acti- 
vation energies from the average value of activation energy. In eq. (5.2), the DExp is 
multiplied with the density function of activation energy. The mean value as well 
as standard deviation of distribution function is constant; therefore, the shape of 
the whole integrand relies on the study variables E, and E. There must be the lim- 
its that can modulate the form of the model through changing the location of maxi- 
mum of the total integrand. Thus, it provides two different ways of solving DAEM. 
The first case: when the standard deviations of continuous distribution functions 
(listed in Table 5.1) are relatively wide to Ew of DExp, the total integrand is the dis- 
tribution function, but as time increases it is progressively truncated from the left 
side by the sigmoid-like DExp. The location of maximum of the total integrand can 
significantly shifts toward the left side, and the shape of integrand gets skewed. On 
the other hand, the relatively narrow form of the continuous distribution function 
to the width of DExp makes the shape of total integrand like the distribution func- 
tion, but the amplitude is drastically reduced by DExp as time proceeds. The nature 
(symmetrical or asymmetrical) of the distribution function remains unchanged; 
however, the global maxima of the whole integrands get shifted. Conversely, pyrol- 
ysis is not solely decomposition of the biomass, but the simultaneous degradation 
of various constituents (moisture content, hemicellulose, lignin, and cellulose). 
Therefore, it is difficult to describe pyrolysis through a single distribution function. 
It can also be an alternative to combine the two limits together for two distinct re- 
gions of experimental curve. But it is more important to look at the behavior of ther- 
moanalytical data sets rather than applying *hit-and-trial" scheme. The data set is 
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highly unpredictable during devolatization stage [29] and the reason is some fast- 
autocatalytic reactions that trigger the secondary pyrolysis, which is difficult to 
model through the conventional form without any transmutation, or variation of 
limit on o (standard deviation), as the dependency of such reactions relies on the 
residence time of volatiles [25] inside the cellulose matrix; however, there are also 
some other factors, such as grain orientation, size of biomass particle, structural ef- 
fect, and moisture content. Thus, the model must be robust enough to explain the 
transition states that depend on the residence time of volatile. Although the release 
of volatile content changes if heating rate varies from ramping to constant tempera- 
ture cases, so there is no appreciable variation encountered in the amount of re- 
leased volatile until the two parts of the DAEM integral overlap completely [3]. 

DExp as well as distribution together must be varied with time and tempera- 
ture. Analysis is carried out for the relatively narrow distribution function (o — 0). 
The development in the subsequent sections pivoted around the method of the 
moving maximum for the Laplace integral problems (Refer Appendices). The behav- 
ior of the density function as well as DExp is examined at the interior points of acti- 
vation energy domain (p"(E) » 0). The objective of implementing technique is to 
determine the moving maxima of the whole integrand as a function of time so that 
the TG data can be interpreted easily. 


5.2.2.1 Asymptotic scheme 


The postulation for the narrow distribution of the continuous distribution of types 
(Table 5.1) is taking standard deviation o = RT|, where T; denotes temperature at the 
inflexion point and follows the given thermal history. The local maxima y* of corre- 
sponding functions (M(y), J(y), U(y), and H(y)) of distribution functions is sought by 
finding the point where M'(y), J (y), U(y), and H'(y) are zeros. Asymptotes and 
shapes of transmuted function are decided by following the two limits [49]. 


Proposition 5.1. The asymptotic of eq. (5.1) if p(E) — O are 


(i) p (E) (1 t)p(E) 
(ii) G'(E)— (14 1)G(E) 


Proposition 5.2. The asymptotic of eq. (5.1) as E — ce are 


(i) p' (E) -(1- 1)p(E) 
(ii) G (E) -(1- 1)G(E) 


Modeling and analyzing material properties in the field of engineering are very com- 
plicating and cumbersome. Sometimes, the data is too big to predict precisely through 
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the conventional distribution function, especially the life prediction of machine 
parts. The density functions are widely implemented for studying fatigue and endur- 
ance life of engineering device and materials. Despite the wide range of application of 
distribution function, sometimes it is unable to explain different types of data that has 
nonmonotonic attribute; for instance, conversion of material at different thermal his- 
tories. Therefore, it is required to develop a flexible parametric model that can appro- 
priately be embedded into a large model by adding an additional shape parameter. A 
significant contribution has been carried out for transmutation of the conventional 
model, thus enhancing flexibility in modeling of various types of life prediction data 
of the engineering components. The new transmutation models are developed and 
comprehensively perused in some of literatures [35, 36, 38, 49]. Therefore, some of the 
distribution types are modified to predict the thermoanalytical data efficiently through 
activation energy modeling. 


5.2.2.1.1 Case A: transmuted Lindley (refer Appendix A) 
Invoking eq. (5.1) and replacing the distribution function with transmuted function 
p'(yi) by proper transformation in eq. (5.2), we have 


1-a ~ | e| -exp (=>) ae (5.10) 
0 


Putting the value of p‘(y,), we get 
Ysı -yı 
1-a ~y | exp{ — exp (sch ) HDD +(1-k) “on (5.11) 
W. 
0 


Equation (5.11) can be written as a convex combination of functions fı (y1) and Eil: 


ZEE e (2 UE ljexp{ - w (yı - 1)] 


+ EEN 


(5.12) 


where k is a mixing fraction. 
Rewrite eq. (5.12) into two different parts and solve it for the same value of y: 


1-a-[((1-ag) * (1-a1)] (5.13) 


Solve eq. (5.13) separately for 1- a. 
Let, 


s i 
1 a= [55 e| exp(” a d ljexp{ - wı(yı - 1)}dyı (5.14) 
0 
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Suppose that 


Vii) = exp (125) (yi - 1) 


The function Vin) is at its extremum if V’(y,) = 0. 
Differentiating eq. (5.15) w.r.t y;, we get 


Nader um 
Vamos -u 
Ywı 
Let the location of maximum is denoted by y,°, then 
yit —ys-Ywi In(@1Yw1) 


The second derivative of eq. (5.15) w.r.t. y; is given by 


Replacing Vi) in eq. (5.15) by its Taylor expansion, we have 


VQ) Vy*) + HD , CT Wii) eee 


After simplification, we get 


(5.15) 


(1- a9) ~1.25(yi° -1)k (vo ran) D D oša 4) exp{-@i(Ywit+ yif —1)} 


Now consider the other term (1 — a4) 


(ET ZE | a4 (1 — ver) —exp =>) Joen -wi(y-1)}dy 
w1 
0 


Similarly, the required expression for (1— a) is, 


(5.16) 


(5.17) 


(1-24) ~1.250(1-k) (eet! D D Gi + j exp{ - ai(ywi +y: — 1)} 


(5.18) 


Substituting eqs. (5.16) and (5.17) into eq. (5.13), we have the resultant solution for 


transmuted Lindley function at different limits 
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p(E)—0 wi 


lim (1-a )-125 Da sei (nt D Gg 4) 
exp{ - 01 (Yui + y? - 1) (0 wi(yi® - 1)} 


Jim (1-a)-125 (1-7) (1- ovii (eere E + 1) 


exp{ - wi (Ym +y? -1) {0 + wy - 0] 


5.2.2.1.2 Case B: transmuted Weibull (refer Appendix B) 
From egs. (5.1) and (5.2), we get 


1-a- | wen) -exp (5 2) Jena (5.19) 
J Yw2 


Substituting the density (p(y2) from Table 5.1 into eq. (5.1) for the Weibull distribu- 
tion, we have 


1-a-(1* nos | (ya - jte -exp =>) o (0 z uf } dy (5.20) 
0 


The final expression after transmutation is expressed through eq. (5.20). 
Assume that M(y2) is a function whose maxima is at y; = y2°, then 


M(y2) = - exp =>) i (w0 - uf (5.21) 


Find the first- and the second-order derivatives of eq. (5.21) and substitute its Taylor 
series expansion into eq. (5.20): 


Equating M’(y2) = 0, we get 


y =14 (B-1)ywW ( g- a (zz) "e sis =) }) 


Likewise, 


exp (s "2 
w. p-2 
yu: B(B - 1) (y; -1) (5.22) 


MO) = 
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e Ir uer 7. @ 
M(y2)~ M(y2°) + tege M Defi. (5.23) 
After substituting M (y2) from eq. (5.22) into eq. (5.23) and replacing M(y;) by its 
expansion in eq. (5.20), we get the asymptotic expression for the Weibull distribu- 
tion for different limits: 


mw B(y2° - 1) 


nud ~(1+T) ve IS CH (erf(y2°) +1) 


Yw2 


exp{ - wj (y;* - 1)? (yf + Byw- 


B e 
jim (1- a) —(1- T) ve Gë e Yr (erf(y2°) 
—oo b Mila 

{( Lj 


exp{ wën uf" (yf + Bys2 - 1 


5.2.2.1.3 Case C: Transmuted Gaussian (refer Appendix C) 


Replacing the Gaussian distribution from eq. (5.2) with transmuted Gaussian func- 
tion, we have 


1-a- | Eoexpf -exp (” = d bound (5.24) 
WwW 
0 


After putting the value of the Gaussian distribution from Table 5.1 into eq. (5.1) and 
substituting into eq. (5.24), the transmuted form will be represented as 


i-a ien | exp{ exp (1523) w3(V3 y] dy; (5.25) 
0 


Suppose that 


Ilys) = ep (52) (ys - 1? 


The first and the second derivatives of J(y3) w.r.t “y3” are given by eqs. (5.26) and 
(5.27), respectively 
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U” O)| = 205 K ` a 1) (5.27) 


Ww 


The location of maximum of function J (y3) is obtained by J’ (y3) = 0. Hence 


1 Ys — 2) 
7 — Lt Ma WI ex 
5 CN (zia d Yw 


Substituting the values of 
we have 


J” (y3°)| in the Taylor series expansion of J(y3) about y3°, 


a (y -yst 


BI rn, — 


J(y3)- Jy? 


After simplification and rearrangement, the required solution of transmuted 
Gaussian function is expressed as 


exp{ - W3(y3° — 1)(2yw + y5^ — 1) } 


lim (1-a)~(1+T)(erf(y3°) +1) for p(E) ^ 0 
p(E)0 (m 1 
yw T ) 
m (1- a)~(1- r)(erf(ys^) + 1) Sut «dn Og" Uta Up) for E — oe 
m cm 


5.2.2.1.4 Case D: Transmuted Frechet (refer Appendix D) 
The expression of transmuted Frechet is obtained after replacing density function 
from transmuted density function. 

Equation (5.28) is the required expression which is subjected to investigation: 


1-a~ | H exp - exp (==) yuan (5.28) 
4 Yw3 


Replacing transmuted expression through substituting the value of density function 
from Table 5.1, we have 


1-a~ (14 ||) | exp Il - exp (==) Jeu "2 -4) 6*9 e- aal 2 T len, 
0 
(5.29) 


To solve eq. (5.29), suppose that 
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U(ya) = exp (12) {wa Du -1)) ^ (5.30) 
1-a- (1er) Lenin, (5.31) 


0 


Obtain the first- and the second-order derivatives of U(y,) as it is performed in the 
other cases: 


ys -ya* 
; exp » = 
ue, "SE ) - |o; Reg - 1) -! 
W. 


The location of maximum value is derived by U’(y,) = 0. 


1 
1 1 lp - 1| ys-1 
e =1+YyYw3 W ( ) exp( ) 
d CS = Ywslplw4 P! B Yw3 


After simplifying the second-order derivative of eq. (5.30), we get 


" , e 1 lp -1l 
"volo ene nue (927 ve 


Expand U(y,) with the Taylor series expression about y4 = y4*, 


2 n 


(y — af U (DI dues (5.32) 


U(ya)~U(ya") + 2 


After replacing U(y,) from its Taylor expansion given by eq. (5.32), we get the solu- 
tion for the transmuted Frechet as 


e 4l, lol 
lim (1-a)~(1+1)0.707 Ava D 0h 


us e 4y 0-1 
pta (Cram 


Yw3 


(erf(ys*) +1) 


ep) wën ET Ywslp| -o«-15]] 


lpin(ys* -1)P 
ej lp -1| 
(E52 ien 


exo — n P (ya? - 1)" (yuslo| - (Va? —1)) 3 


Jim (1- a) —(1- 1)0.707(erf(y?) +1) 


Note: Valid for the negative shape parameter (p « 0). 
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5.2.2.1.5 Case E: Transmuted Rayleigh distribution (refer Appendix E) 
Two-parameter Rayleigh distribution is used for demonstrating the scheme. The 
transmuted distribution function is obtained by replacing the density function from 
its expression given in Table 5.1. Equation (5.33) is the transmuted expression for 
first-order chemical reactions: 


1-a ~ eel el ECH posay (5.33) 
9 w4 


After obtaining p‘(ys) from eq. (5.1), we get the expression for the expected solution 


1-a~ | exp{ -exp (8) pus (y 1)exp{ vis (Dis 1) bays (5.34) 


Assume that 


H(ys) = exp) viii, - 1) 


w4 


Then eq. (5.34) can be represented as 


co 


1-a- | 2wsdly-1)exp(H(ys) dy (5.35) 
0 


The first derivative of H(ys) w.r.t. “ys” is 


For deriving the location of maximum, equate H'(y;) = O 


1 ys, — 1 
$214 w( ex ( ) ) 
d dii 20sywa? P Y wa 


“yse” is the location where the function H(ys) attains its maximum. 
Similarly, after simplification, the second derivation is obtained as 


24 5 E 1 
H’(ys°) = B; (5 — k 1) | (5.36) 
Yw4 
The Taylor series expansion of function H(ys) about ys* is given by 


c s 2 "74 
Hos) Host) « 9-397. se). (5.37) 
After putting the value of H” (ys?) from eq. (5.36) into eq. (5.37) and replacing H(y;) 
in eq. (5.35) with its Taylor series expansion (eq. (5.37)), we have the required ex- 


pression for the transmuted Rayleigh distribution: 
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lim 1-a~(1+T) 
p(E)—0 


lim 1- a—(1- 7) 


E—oo 


5.3 Application and computation analysis of TG data 


Elemental analysis of sample subjected to thermal degradation is conducted through 
CHNS-O analyzer. The CHNS-O analyzer (vario MACRO) is used to evaluate the chemi- 
cal composition of pine needles. Before using the CHNS-O analyzer, the furnace is 
heated up to 1,473 K for 30 min. Once it gets heated, the capsule forms of sample are 
fed into a slot of rotating disk. The flow of oxygen maintains the catalytic combustion, 
whereas helium gas acts as a carrier gas stream so that the volatile gases with water 
(CO>, SO;, NO,, and H;O) can be absorbed at different reduction columns. These tubes 
are in-between combustion chamber and detection unit (computer-based software), 
where the gas mixture is separated into its components through purge or trap chroma- 
tography. Thereafter, each component is separately detected by thermal conductivity 
detector (TCD). The gases CO, HO, and SO, are absorbed in a sequence, whereas N; 
passes through all the three columns and it is finally detected by TCD. TCD generates 
the electrical signal that is proportional to the concentration of elementary compo- 
nents of biomass. The reference material is considered to be the powder of birch leaf. 
The bomb calorimeter at constant volume is used to measure the calorific value of 
pine needles. Thermal evaluation of pine waste is based on TG differential thermal 
analyzer (Diamond TG/DTA, Perkin Elmer, USA). The sample of 5.68 mg is placed in 
an electrobalance from 303 to 923 K. To avoid the measurement error due to buoyant 
effect at the high temperature, the horizontal differential-type balance is considered 
for experimentation. Thermocouple-type “R” is referred to measure the temperature 
inside the furnace. The volumetric flow of noninteractive gas is set at 200 mL min”. 
Indium and tin are used as reference materials for differential thermograms. 

Table 5.2 illustrates the elemental characteristics and higher heating value of 
pine waste on the dry basis. TG data of pine sample is used to implement the trans- 
mutation of DAEM, so that qualitative assessment of different distribution functions 
can be performed. Computational as well as algorithm designing is carried out on 
MATLAB 2017b. The nonlinear thermal history is tested for narrow distribution na- 
ture of density functions. The different assigned values of parameters are repeatedly 
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Table 5.2: Elemental composition of pine needles. 


C% H% N% 096 S% Ash content (%) Higher heating value (MJ/kg) 
53.7 6.48 0.56 36.13 0.10 2.92 22.23 


examined in every iterative loop until it does not satisfy the given condition. If any 
parametric value does not validate the constraint, system searches for different val- 
ues. Once the solution is obtained, the program terminates itself. 


5.4 Results and discussion 


The numerical solution of transmuted modeling of DAEM is performed by using the 
Laplace integral method. The proposed scheme is exhaustively examined for differ- 
ent distribution functions. The significant role of moving maximum with time on 
the energy distribution pattern of released volatile is illustrated graphically. There 
are three different limits imposed on each distribution case and the results are de- 
rived separately after transmutation. 


5.4.1 Effect of transmutation on the Lindley distribution 


The effect of transmuted Lindley distribution is shown in Figure 5.1(a)-(f). The sensi- 
tively of the Lindley distribution with respect to time and temperature is relatively 
higher than that of other types of distribution functions analyzed for the energy model- 
ing purpose. There is a steep increase in the conversion rate of biomass with variation 
of kinetic parameters in Figure 5.1(a)- (f). Elevation of frequency factor causes the re- 
maining mass proportion curves to shift upward with relatively high (Table A.1) resid- 
ual mass of biomass. However, the variation of residual mass of biomass is controlled 
to some extent through transmutation, but relative change in the residual mass of bio- 
mass is negligibly small. Transmutation of the Lindley distribution infers that the 
higher activation energy of 103 kJ/mol can be easily modeled through Lindley distribu- 
tion. With increase in activation energy, the remaining mass fraction gets shifted 
down; however, there is no significant change in the attribute of (1-a) curves with the 
variation in the thermal history of sample. The maximum upper limit of activation 
should not be more than E. « 5.72 x 10? kJ/mol, while the average activation energy 
should not be more than 103. Comparative sketch of the transmuted Lindley distribu- 
tion with the classical DAEM is shown in Figure 5.2. The fact concluded through the 
Lindley distribution is that Lindley distribution is not able to describe the pyrolysis 
mechanism more effectively than the other distribution types, despite having the 
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Figure 5.1: The effect of transmutation on the numerical solution of Lindley-based DAEM: (a)- (d) 
kinetic parameters; e and (f) characteristic of nonlinear thermal profile. 
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Figure 5.2: Comparative sketch of t-DAEM (Lindley) with the conventional DAEM form. 


highest value of coefficient of regression. It may be possible that the Lindley distribu- 
tion is not able to extrapolate pyrolysis effectively, and behavior is like the phenome- 


nological models. 
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5.4.2 Effect of transmutation on the Weibull distribution 


Unlike the Gaussian distribution, the Weibull distribution is abounded with different 
parameters, which is also one of the reasons of higher redundancy than the Gaussian. 
The effect of transmuted Weibull distribution is illustrated in Figure 5.3(a)- (f). The ther- 
mal decomposition of biomass takes some time to reach the critical temperature for ini- 
tiating the process. The Weibull distribution is unable to distinguish the response time 
at different activation energies; therefore, all the curves are highly overlapped without 
any significant variation at the end states of pyrolysis (Figure 5.3(a)). The effect of dif- 
ferent limit on the transmuted (Table B.1) Weibull distribution is depicted in Figure 5.3 
(b)-(f). As the frequency factor increases, the remaining mass fraction curves get 
shifted upward, which results in shifting of critical decomposition temperature. Unlike 
the transmuted Gaussian distribution, the curves are converging for the lower 
range of activation energies, or the relative width between activation energies of 
parallel reaction must be narrow for the better simulation of pyrolysis through 
Weibull distribution. At the elevated heating rate with time, the devolatilization 
stage gets shifted to the left. The similar trend is obtained for the initial experimental 
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Figure 5.3: The effect of transmutation on the numerical solution of the Weibull-based DAEM: (a)- 
(d) kinetic parameters; (e) and (f) characteristic of nonlinear thermal profile. 
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temperature. Although the Weibull distribution is not able to interpolate the char for- 
mation stage effectively, despite the devolatilization stage effectively described the 
transmuted Weibull distribution. But the objective must not be overfitting of simu- 
lated solution to the thermoanalytical data, as unfamiliarity with the mechanism of 
autocatalytic reactions and the residence time of volatiles may provide erroneous so- 
lution. The maximum upper limit of activation energy for the transmuted Weibull dis- 
tribution must not be more than 130 kJ/mol and the average activation energy (kJ/ 
mol) should be within the domain of (69, 89] for the better convergent rate. 
Comparative demonstration of two different schemes is shown in Figure 5.4. The 
transmuted form of Weibull describes the devolatization stage of decomposition of 
biomass efficiently, but the residence time of volatile content increases with tempera- 
ture and the whole biomass is transformed into volatile gases (Figure 5.4). The char 
yield became asymptotic with temperature. 


Experimental TG data 
— Transmuted Weibull distribution 
Weibull distribution 


300 400 500 600 700 800 900 
Temperature (K) 


Figure 5.4: Comparative sketch of t-DAEM (Weibull) with the conventional DAEM form. 


5.4.3 Effect of transmutation on the Gaussian distribution 


The comprehensive effect of transmutation on the Gaussian DAEM is depicted 
graphically in Figure 5.5(a)- (f). At the initial stage of pyrolysis, the remaining mass 
fraction (1-a) of biomass is always near to 1, unless there is interactive atmosphere 
or adsorption, or surface oxidation processes occur during thermal degradation of 
biomass. But the graphical representation of the classical DAEM (Figure 5.5(a)) is ex- 
actly 1, which is quite idealistic for pragmatic pyrolysis condition. As the frequency 
factor increases, there is an appreciable displacement of the upper segments of re- 
maining mass fraction (1-a) than that of lower halves at the constant rate of change 
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Figure 5.5: The effect of transmutation on the numerical solution of the Gaussian-based DAEM: (a)- 
(d) kinetic parameters; (e) and (f) characteristic of nonlinear thermal profile. 


of heating rate with time. Despite drastic variation of distribution parameters, con- 
verging ability of the numerical solution is improved through transmutation. Thus, 
the homogeneity of shape is maintained throughout the process. The most important 
facet of transmutation is illustrated in Figure 5.5(b) and (c). The qualitative aspect of 
transmutation of distribution of activation energy is known through the characteristic 
of explanatory variable rather than the coefficient of regression analysis. The overlap- 
ping of curves happens if the density function of activation energy approaches zero, 
whereas the nonconvergent solution converges to solution, if activation energies 
spectrum is very broad for the Gaussian distribution; therefore, the Gaussian distribu- 
tion can be the best fitting model for those reactions (Table C.1) where the activation 
energy variation among reactions is colossal. The effect of initial temperature gets de- 
volatization stage shifted to the right without any significant variation in “drying” as 
well as "char formation" stages. The residual mass obtained after simulating the pyrol- 
ysis process is approximately the same with + 1.896 difference. The effect of activation 
energy variation is depicted in Figure 5.5(d). The increase in the activation energy 
causes the intermediate stages to shift left. According to the analysis, the maximum 
upper bound of activation energy should not be more than E~ x 185 kJ/mol and the 
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maximum permissible limit for the average value of activation energy (kJ/mol) must lie 
in the domain of [120, 150). Inclusion of additional parameters supports the conver- 
gence rate of the numerical solution obtained by the Gaussian distribution. The effect 
of transmutation is mainly on the intermediate states of decomposition, which plays a 
pivotal role in the kinetic mechanism of pyrolysis. However, the influence of transmu- 
tation on the end stages of pyrolysis is difficult to interpolate, yet the variation gap and 
derivation of the classical DAEM from practical aspect of pyrolysis are overshadowed 
by transmutation and the form becomes more homogeneous than the classical form of 
DAEM. The nature of Arrhenius parameters and its effect on the thermal decomposition 
of biomass can be easily distinguished through transmutation. Comparative analysis of 
transmuted DAEM with classical DAEM is depicted in Figure 5.6. It clearly demarcated 
the controlling ability of transmutation to describe pyrolysis process far more effec- 
tively than the classical DAEM. 
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— Gaussian distribution 
—- Transmuted Gaussian distribution 
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Figure 5.6: Comparative sketch of t-DAEM (Gaussian) with the conventional DAEM form. 


5.4.4 Effect of transmutation on the Frechet distribution 


Graphical representation of the transmuted Frechet distribution is illustrated in 
Figure 5.7 (a)- (f). The increase in the frequency factor increases the conversion rate, 
which results in residual mass approach to the experimental condition at the reduced 
timescale. However, the process is accelerated through introduction of linear mixing, 
and the remaining mass curves cause to shift upward as E — 0. Therefore, the trans- 
muted Frechet holds validation for set of reactions having relatively low variation of 
activation energies. The upper limit of activation energy is 190 kJ/mol, while the aver- 
age activation energy (E kJ/mol) must not be more than E, x 120. Activation energy 
modeling through transmuted Frechet distribution provides good result only at lower 
initial temperature, and the variation of heating rate must not be very high. The 


92 —— Alok Dhaundiyal, Suraj B. Singh 


1, = 0.0008 (at Aj) 
Ty = 0.2 (at Aj) 
— t, = 0.008 (at Aj) 
t, = 0.02 (at A5) 
- —T, = 0.008 (at A3) 
“1 = 0.2 (at Aj) 


400 500 600 700 
Temperature (K) Temperature (K) 


—T = 0.05 (at Aj) 0.8 


T, = 0.04 (at Aj) — Ej4 7 176 kJ/mol 


* 0.6 — t, = 0.05 (at A;) — 0.6 E2 = 190 kJ/mol 
d — t = 0.04 (at A) i E43 = 233.45 kJ/mol 
= 04 2 E = 04 
T, = 0.05 (at A3) 
0.2 — 1 = 0.04 (at A3) 0.2 
0 m. i 
300 400 500 600 700 800 900 300 400 500 600 700 800 900 
Temperature (K) Temperature (K) 
() 1 O t 
0.8 04 = 0.0175 °C-min? Kä 
` — a, = 0.0155 *C-min? 0.8 \ 
— 0.6 a} = 0.0125 "C-min? ` 06 
F SH 
= 0.4 x SV 04. Ty, = 400K S 
S —T95 = 303 K > 
0.2 02. Toz = 270K SS 
0 0 — 
300 400 500 600 700 800 900 300 400 500 600 700 800 900 
Temperature (K) Temperature (K) 


Figure 5.7: The effect of transmutation on the numerical solution of Frechet-based DAEM: (a)—(d) 
kinetic parameters; (e) and (f) characteristic of nonlinear thermal profile. 


predicted transmuted DAEM solution with respect to the classical DAEM is shown in 
Figure 5.8. The transmuted Frechet distribution does not provide much better conver- 
gence than that of the Gaussian, the Weibull, and the Rayleigh distributions. 
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Figure 5.8: Comparative sketch of t-DAEM (Frechet) with the conventional DAEM form. 
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Comparative analysis of Laplace temperature with other numerical techniques is 
based on the attribute of numerical solution. However, the Copula distribution [1, 23] 
and 2-DAEM [5] are alternative forms of linear mixing of two different functions, but 
the emphasis is laid upon fitting of data rather than behavior of kinetic parameters 
with respect to time and the thermal history. Dynamic nature of inflexion points in the 
wide distribution case [21], and neglecting the time factor in the kinetic modeling, 
may provide a better fitting but the measured values would be quite far away from the 
pragmatic solution. Another way of solving the double integral of DAEM is through 
Simpson methods [21]. But without knowing the error component while spacing the 
activation energies leads to the fallacious inference. The concept of analyzing the ther- 
mal data is based on the stability of local maxima of curves and removing the inflex- 
ion point H”(ys°)#0 out of the numerical solution. The derived solutions for the 
transmuted distributions rely upon the attribute of y* with time. In the beginning, the 
moving maxima remains constant with temperature, but later it increases linearly as 
time proceeds. The narrow distribution is not able to interpolate the initial decomposi- 
tion phase, as the Lambert W function is almost dormant in the beginning stage of 
simulation, while the Lindley distribution is related to the nature of logarithms. 
Therefore, spontaneity of reactions modeled through Lindley reactions is quite differ- 
ent from other transmuted forms. This problem can be tackled if the similar experi- 
ment is conducted at the complex thermal history. The obtained values of kinetic 
parameter for pine needles are homogenous with the isoconventional methods with- 
out any compensatory effect [29], which is a positive benchmark of the Laplace inte- 
gral. However, the applicability of scheme is limited to only those functions that have 
the nonzero second derivatives and exhibit nonmonotonicity with time or tempera- 
ture. The effect of proposed scheme is based on the significance of the whole inte- 
grand about its central value, which is overlooked in the previous studies [2, 21]. 


5.4.5 Effect of transmutation on the Rayleigh distribution 


The effect of the transmutation on the activation energy modeling through the 
Rayleigh distribution is depicted in Figure 5.9(a)-(f). The resolution capacity of 
Rayleigh distribution with respect to variation of parameters is far more than that of 
the Gaussian and the Weibull distributions. One can observe a fine demarcation of 
behavioral variation of (1-a) with the frequency factor as illustrated in Figure 5.5(a). 
The increasing value of frequency increases the dehydration temperature and the 
(1-a) curves are deviated to the right side without any variation in the char formation 
temperature. However, the deviation of the remaining mass fraction is controlled 
through transmutation of Rayleigh distribution (Figure 5.9(b)-(f). Likewise, the 
Gaussian distribution and the Rayleigh distribution support the reactions that occur at 
higher activation energies. There is no significant variation in the intermediate stages 
of decomposition except the time of initiation of reactions that increases with 
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Figure 5.9: The effect of transmutation on the numerical solution of Rayleigh-based DAEM: (a)- (d) 
kinetic parameters; (e) and (f) characteristic of nonlinear thermal profile. 


increasing activation energies, which is obviously much practical. The upper bound of 
activation energy variation for the effective representation of pyrolysis reaction is 
220 kJ/mol, whereas the feasible (Table EU boundary limits of variation in average 


activation energy (kJ mol’) is within the range of (96, 127]. Increasing the initial ex- 


perimental temperature causes the remaining mass fraction curves (1-a) to shift up- 
ward with higher residual mass than that of lower initial temperature. Similarly, the 
decomposition temperature at the initial stage of pyrolysis is elevated with drastic vari- 
ation of heating rate with time (Figure 5.9(e)). However, perturbance of Rayleigh distri- 
bution due to ingression of linear mixing does not affect the form of distribution 


drastically (Figure 5.10). 
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— Transmuted Rayleigh distribution 
Experimental TG data 
— Rayleigh distribution 


0 1 1 — 
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Temperature (K) 


Figure 5.10: Comparative sketch of t-DAEM (Rayleigh) with the conventional DAEM form. 


5.5 Conclusion 


The narrow distribution of transmuted function is adopted to derive the numerical 
solution. Transmutation of the Gaussian distribution, the Weibull distribution, and 
the Rayleigh distribution provides the better results than that of transmutation 
form of Lindley and Frechet distributions. Every distribution function has its own 
unique attribute for representing the variation of activation energies. The Gaussian 
distribution and the Rayleigh distribution are found to be good for developing the 
model for higher activation energy reactions of maximum of 202 kJ/mol, while the 
Weibull distribution holds good relationship for the activation energy range of 80- 
130 kJ/mol. The modeling of activation energy through Frechet distribution is valid 
for the negative value of shape factor, whereas Lindley has anomality with the mea- 
sured value of pine needles from the isoconversional method. However, there is a 
likelihood of existence of those reactions whose activation energy is of 103 kJ/mol. 
The maximum value of activation energy for the transmuted Frechet distribution is 
found to be 190 kJ/mol. The frequency factor for all the cases is estimated to lie in 
domain of 103-1,013 min '. The variation of activation energy with time, tempera- 
ture, energy spacing, and thermal profile and the significant contribution of the 
whole integrand are investigated through knowing the attribute of moving maxi- 
mum of the whole integrand. However, the computational effort of scheme is rela- 
tively high to other schemes [1, 23, 50], as algorithm evaluates the location of 
maximum with respect to time. 
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Appendices 


To demonstrate the scheme, the Laplace method of integration is chosen to simplify 
the DAEM problem of iterative computation of double integrals. This method is use- 
ful to solve the integral of the form 


I- | e" BO s(x) dx 


a 


where a and b can be finite or infinite. Here Q(x) and s(x) are replaced by their 
local expansion when x - xo. The same concept is used for approximating the inte- 
gral form as mentioned earlier. 

The general solution of Laplace method for interior points is given by 


1 


I1-17724e- Poo) s(x,) Bro 
0 


Note: s(xo)z0 
However, c(ye) < 1, so its contribution to integral is negligibly small for inter- 
val of 0 « ye « 1. 


A. The Lindley (three-parameter) distribution 


Table A.1: Kinetic parameters for the Lindley distribution functions. 


Average activation energy Frequency factor R (p(E) — 0) (E — ee) 
(kJ/mol) (min?) After transmutation After transmutation 
T R? T R? 
1.0543 x 10^ A 2x10'^^ 0.95 0.0008 0.93 0.05 0.94 
0.02 0.94 0.04 0.95 
1.4906 x 10^ A; 2x10? 0.93 0.0008 0.94 0.05 0.94 
0.02 0.94 0.04 0.93 
2.7733 x 10? As; 2x10 0.96 0.0008 0.97 0.05 0.97 
0.02 0.97 0.04 0.98 

oi 

p(E) (1+ y(E- y)e TEL E> y>0,0+ >0 (A.1) 


SEET 
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This can be solved through the convex combination of two different functions, but 
provided that vfi (E), f(E)eC and CC R 
Equation (A.1) can be rewritten as 


P(E) = kk (E) + (1- KE (A.2) 
where k is the mixing fraction and o is a scale parameter 


0 


- _1 E 
CET. Fc and ET 


The mean and the variance of the Lindley distribution are Eo; and ol, respectively 


gaz CWO) quc. Le ye) EE 


CES 8^(6. 6) 


Energy is rescaled as 


y- 


MV | En 


> Ya = ER. and Yw = Ew 
y y 

The expressions for fi(E) and f; 
fi(E) = Op(E - y)exp(- €(E - y) (A.3) 

f(E) = €'exp(- (E- y)) (A.4) 


— 


E) functions are 


The simplified form obtained after incorporating the energy rescaling 
fiy) =w- 1)exp{- wi (yi - 0) 
hı) -€ exp(- ani - 1) 
The resultant form of density function in terms of energy rescaling factor is 


p) 7 kii) +0- kbh) 


Note: The whole integral transforms into convex combination of Gamma functions 


(FQ) EIS un) if (35) — o 


w 
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B. Weibull distribution (three parameters) 


Table B.1: Kinetic parameters for the Weibull distribution functions. 


Average activation energy Frequency factor R (p(E) — 0) (E — eo) 
(kJ/mol) (min?) After transmutation After transmutation 
T R? T R? 
89 A 10? 0.920 0.0008 0.94 0.05 0.90 
0.02 0.96 0.04 0.91 
99 A» 10? 0.910 0.0008 0.94 0.05 0.89 
0.02 0.95 0.04 0.90 
69 A3 107 0.900 0.0008 0.93 0.05 0.88 
0.02 0.95 0.04 0.89 
B 
(B-1) E-y 
TUNE 
p(E) == | — e" (B.1) 
n n 


This is the expression for the Weibull distribution function for three parameters. 
Energy rescaling factor is given as follows: 


E E E 
Y2=—, ya 7 — and ym = — (B.2) 
y y y 


After simplifying (B.1.), the required form of density function p(y2) is 
P(Y2) = Pot ya -1fle (202-9) 


where 


The mean and the variance of the Weibull distribution are Eo; and 0°, respectively: 


nl reell 
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C. Gaussian distribution 


Table C.1: Kinetic parameters for the Gaussian distribution functions. 


Average activation energy Frequency factor R (p(E) — 0) (E — eo) 
(kJ/mol) (min?) After transmutation After transmutation 
E A 
w T R? T R? 
120 A, 2x101? 0.9315 0.0008 0.93 0.05 0.97 
0.02 0.92 0.04 0.94 
119.6 A2 2 x 107 0.9308 0.0008 0.93 0.05 0.94 
0.02 0.92 0.04 0.94 
150 A; 2x10? 0.9307 0.0008 — 0.93 0.05 0.96 
0.02 0.92 0.04 0.95 


The required expression for Gaussian distribution is given as follows: 


Ge) 
p(E)- ut 26 (C.1) 


The energy correction factor for the Gaussian distribution is expressed as 


E E E 
=g. er — and yw = — (C.2) 
0 y y 


Simplifying eq. (C.1) through incorporating energy factors, we have 


| oil Ey V 2 
2 
Putting w3 = EN in eq. (C.3), the required form is 


pos) exp w3(y3 1) 
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D. Frechet distribution 


Table D.1: Kinetic parameters for the Frechet distribution function. 


Average activation energy Frequency factor R (p(E) — 0) (p(E) — 0) 
(kJ/mol) (min?) After transmutation After transmutation 
r R? T R? 
119.6 Ay 2x10°° 0.94 0.008 0.93 0.05 0.96 
0.02 0.92 0.04 0.95 
140 A; 2x10? 0.80 0.008 0.80 0.05 0.79 
0.02 0.81 0.04 0.80 
171 A, 2x10? 0.15 0.008 0.16 0.05 0.13 
0.02 0.97 0.04 0.14 

p (E-y -(1+p) (race 


Equation (D.1) represents the density function of the Frechet distribution. 
The energy factor for distribution is given by 
E Es Ew 
Ines Yas = — and ya = — (D.2) 
H" p > y 


Introducing the energy rescaling factors in eq. (D.1), we have 


5 -(1+p) NE 
ES RI cu (D.3) 


Putting w4 = E in eq. (D.3), the desired expression will be 


pya) = Bwn -P+D (y, - 1) - (+0) o- {wa (vs -1)) ^ 
The mean and the variance of the Frechet distribution are Eo, and ol, respectively, 


1 
Es -uenr(1- 7)» 1) 


céiert 


104 — Alok Dhaundiyal, Suraj B. Singh 


E. Rayleigh distribution 


Table E.1: Kinetic parameters for the Rayleigh distribution functions. 


Average activation Frequency factor R (p(E) — 0) (E — ee) 
energy (kJ/mol) (min?) After transmutation After transmutation 
T R? T R? 
127 A 2 x 108 0.94 0.0008 0.940 0.05 0.95 
0.02 0.935 0.04 0.96 
79 A; 2 x 10° 0.93 0.0008 0.941 0.05 0.94 
0.02 0.937 0.04 0.95 
234 A3 2x10?! 0.93 0.0008 0.930 0.05 0.95 
0.02 0.92 0.04 0.94 
The density function of two parameters (2-P) is expressed as 
p(E) - 2E - p)exp{-x(E-$)"} (E.1) 
The energy factor for the Rayleigh distribution is expressed as 
E E, Ew 
ys7 4 Yos=— and yw - —— 
p >| $ "A 
Rearranging eq. (E.1), we will have the required form of density function: 
Ws 2 
p(ys) =2 $ (Ys 1)exp{ ws(ys - 1) \ (E.2) 


where 


p 


"5" 1-A) 


The mean and the variance of the Rayleigh distribution are Eos and gie, 


respectively, 


Jn 1 Di 
E S andio e se” 
pat pe gs X X 


where dh kJ/mol is the location parameter 


The physical significance of location parameter is the minimum energy re- 


quired to initiate the chemical reaction. 


Shruti Tikhe, Kanchan Khare, S. N. Londhe 
6 Air quality forecasting using soft 
computing techniques 


Abstract: Air pollution is a fundamental problem among the critical challenges of 
modern society. Acute as well as long-term exposure may pose serious health prob- 
lems. Real-time air quality predictions are necessary to forecast the pollutant con- 
centrations. This information can be used to issue early air quality warnings that 
allow government and people to take preventive measures. 

The authors offer application of soft computing techniques of artificial neural 
networks (ANN) and genetic programming (GP) for air quality parameter forecast- 
ing, a few time steps in advance. Seasonal forecasting models are also useful for 
proactive measures. Availability of continuous field data sometimes is a challenge. 
It is time-consuming and expensive. Modeling options are suggested for situations 
when there are constraints about data availability and the urgent need to under- 
stand and broadcast the air quality. 

A case study of Pune, a metropolitan city in India, demonstrating the strengths 
of GP is presented. We have evaluated the models by a correlation coefficient, 
d-statistics, root mean squared error, and mean bias error. 


Keywords: Air quality, soft computing, ANN, GP 


6.1 Introduction 


Air is a fundamental need for survival and development of all lives on the Earth. It 
has major impact on health and also influences the economic development of the 
country [18]. Air quality describes the chemical state of the atmosphere at a particu- 
lar time and place. Rapid industrialization, development of cities, and concerned 
human activities are the major causes for an alarming increase in air pollution lev- 
els and degradation of the air quality especially in almost all the metro cities of the 
world [23]. Short-term exposure to high levels of pollution may pose severe health 
impacts such as eye irritation, difficulty in breathing, and cardiovascular and pul- 
monary health effects. Health concerns associated with long-term exposure can be 
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listed as the threat of cancer, risk of premature death, and also damage to the 
body’s immune, respiratory, neurological, and reproductive systems. Senior citi- 
zens, kids, and individuals with preexisting heart and lung maladies along with di- 
abetics are at greater risk for air pollution-related health effects [19]. The World 
Health Organization has issued guidelines related to permissible levels of pollu- 
tants with an objective to protect human health and the environment. In spite of 
significant progress in understanding the process of emission and fates of the air 
pollutants as well as in decreasing their ambient levels in urban localities for the 
last five decades, air pollution is estimated to destroy the lives of about 4 million 
people every year worldwide [34]. It is therefore worthwhile to model and forecast 
the air quality. 


6.2 Issues and challenges 


Air quality is described by the concentrations of major air pollutants such as sul- 
fur dioxide, nitrogen oxides, and respirable suspended particulate matter (SO), 
NO,, and RSPM) that reflect the status of air pollution in the environment [36]. Air 
quality assessment studies are essential to determine the trends and effects of air 
pollutants on human lives, animals, soil, and environment. Source contribution 
to pollution can also be identified and accordingly controlling strategies can be 
decided [11]. 

Air pollution comprises a mixture of gases and particles in harmful amounts 
that are released into the atmosphere due to either natural or human activities. The 
leading sources of air pollution are natural as well as anthropogenic. Natural phe- 
nomena such as volcanic eruptions and forest fires release harmful substances 
such as SO;, CO;, NO;, and CO. 

Anthropogenic sources include the pollutants entering from a fixed point, mo- 
bile sources, and area sources in addition to the burning of fossil fuel. Pollutants 
emitted from human-made sources include hydrogen, oxygen, nitrogen, sulfur, 
metal compounds, and particulate matter. 

They affect badly on human health and the environment. The global impact of 
air pollution is on climate change and particularly on global warming. 

During the last five decades, scientific research has been carried out to model 
as well as forecast the air quality. Various methods and techniques have been 
evolved for modeling of air pollutants. For effective air quality management, moni- 
toring and record-keeping of the pollutant concentration data become a vital pre- 
requisite. Many countries have their real-time air quality monitoring as well as 
forecasting systems. Developed countries have their established air quality monitor- 
ing system, but the issue is with the developing countries where availability of the 
continuous data is a significant challenge. 
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6.3 Air quality modeling 


It has been observed that air quality assessment and modeling has been carried out 
using a variety of approaches. The complexity of the problem and the availability of 
the resources decide the choice of method [6]. Traditionally, hard computing ap- 
proaches have been used to carry out air quality modeling and forecasting. They 
include either physical, deterministic, stochastic, statistical, or numerical models 
[30]. Deterministic models require emission inventory and meteorological variables 
to calculate the pollutant concentrations. These models have been used tradition- 
ally to calculate air pollutants in urban areas [20]. 

Physical models are costly, data-intensive, and challenging to construct and 
operate; as a result, their use has been restricted only to understand the complex 
urban dispersion phenomenon [1]. Majority of the dispersion models developed to 
evaluate and predict the pollutants in urban areas are “causal” in nature and 
therefore fail to predict the “extreme” concentrations [9]. In air quality manage- 
ment, forecasting and broadcasting of the peak concentrations (extremes) are 
very important to give sufficient warning to the citizens. Statistical techniques 
represent only a specific monitoring station; as a result, they cannot be extended 
to other regions with different meteorological conditions [14]. The interrelation- 
ship between human, climate, and air pollution is extremely complex; as a result, 
it is very difficult to be represented in deterministic models without developing a 
separate statistical model [6]. As meteorology affects air quality and this relation- 
ship is complex and nonlinear, data-driven approaches are suitable; however, 
they are also site-specific. They provide tools to facilitate the conversion of the 
data sets into a convenient form and convey a better understanding of the pro- 
cesses hidden in these data sets by expressing the mathematical statement of the 
relations between various parameters in the data [2]. Soft computing tools have 
been developed as an alternative to deterministic and physical models. 


6.4 Forecasting of air quality parameters 


Air quality parameters are stochastic in nature; thus, short-term predictions are 
possible based on past meteorological data [22]. Soft computing tools are modern 
forecasting approaches. Resource availability and accuracy requirements decide 
the choice of the tool. Chosen tool should be robust enough to take into account the 
data fluctuations without compromising the accuracy. Extensive experimentation is 
required to identify such a tool [12]. 
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6.4.1 Artificial neural networks 


There are many data-driven approaches known in computer science today. One of 
the most popular approach is artificial neural networks (ANNs). ANNs are biologi- 
cally inspired, parallel computational tools. They consist of simple highly intercon- 
nected processing elements that process the input similar to the human brain [28]. 
ANN is broadly classified as a feed-forward neural network and a recurrent neural 
network. Figure 6.1 shows a feed-forward network [15], which typically has three 
layers, namely input, output and hidden layer. Input layer neurons receive the 
input variables from the problem, and output layer neuron/neurons yield the net- 
work output (desired result). ANN works similar to the human brain on two fronts. 
First one is acquiring knowledge from the input data via the network of neurons 
through a learning process and interneuron connection strengths are utilized to 
store the acquired knowledge. The weights [between the input layer and hidden 
layer as well as hidden layer and output layer] represent the relative strength of 
specific input parameters in the computation of output. An artificial neuron weighs 
and sums up the inputs and compares the results with a predefined threshold. 
Iterations are carried out until the stopping criteria are reached. The network- 
generated output heavily depends on the type of activation functions (also called 
transfer or squashing function) used. Three types of transfer functions are generally 
used, namely, step, ramp, and sigmoid function. The sigmoid function is the most 


Input Hidden Output 
layer layer layer 


SO, (t+1) 


SO, (t-2) 


Figure 6.1: The feedforward network. 
Source: The ASCE Task Committee, 2000 [15]. 
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commonly used transfer function for air quality modeling [13]. The recurrent net- 
work has a feedback path which makes it sequential rather than combinatorial. 
Once cyclic connections are included, the network becomes a nonlinear dynamic 
system. Backpropagation is one of the most widely used algorithms in the feed- 
forward type of networks, which is based on the steepest descent method. 

The performance of ANN models entirely depends on the quality and the quantity 
of the data. ANNs were considered as black box earlier but some researchers have 
shown that the network incorporates all of the physics underlying the process since 
they are embedded in the data and ANN is not merely playing with numbers [25]. 
Similarly, the works of Londhe and Shah [10] and Londhe and Panchang [26] throw a 
light on the work of knowledge extraction from the trained network. Deciding ANN 
architecture is a complex process that demands the skill of the modeler. There is no 
fixed rule to determine the data division and the number of neurons in the hidden 
layer. It has been reported that ANN models may suffer from overfitting due to the 
addition of too many hidden neurons or executing a large number of iterations. Some 
researchers have found that ANNs have a capability of prediction of extreme events 
when combined with other techniques such as wavelets [4]. 


6.4.2 Genetic programming 


Genetic programming (GP) is an evolutionary data-driven modeling approach that 
is similar to genetic algorithms (GA) and follows the Darwinian principle of the sur- 
vival of the fittest and obtains the solution to the problem through the process of 
crossover, mutation, and reproduction. GP can be conveniently used as a regression 
tool as the solution of GP is a computer program or an equation as against a set of 
numbers in GA [24]. 

The process of GP starts with a population of randomly generated computer 
program initially constructed from the data set on which genetic operations are per- 
formed using function and terminal set. Function set comprises the operator to be 
used such as addition, subtraction, logarithm, and square root. Terminal set con- 
sists of values such as inputs, constants, and temporary variables on which the 
function set operates. From the population so generated, four programs are ran- 
domly selected, and a tournament is conducted. 

GP measures the performance of each program, and the two programs are se- 
lected based on performance. GP algorithm copies the two winner programs and 
transforms these copies into two new programs via crossover and mutation as per 
fitness [17]. Creation of offsprings continues iteratively till a specified number of off- 
springs are created in the generation, and until a specified number of generations are 
produced. The resulting offspring at the end of the process is either a computer pro- 
gram or an equation. It is the solution to the problem. The generated equation or pro- 
gram can be directly applied to unseen data to obtain required predictions [25]. 
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Figure 6.2 shows a typical GP flowchart. GP can evolve a formula or an equation 
by automatically selecting the significant inputs contributing to the model devel- 
opment and simultaneously ignoring the variables with the least contribution to 
the model. Three forms of GP are available, namely standard or tree-based GP, 
linear GP (LGP), and automatic induction of machine (AIM) code GP. Tree-type 
GP operates on parse trees, which is constituted by a function set and a terminal 
set. Crossover operation is carried out by randomly swapping the subtrees be- 
tween the selected individuals [2]. The typical GP process is continued until the 
stopping criterion is met. LGP evolves sequences of instructions using an “imper- 
ative programming language” or a machine [8]. In LGP, the expressions of a func- 
tional programming language (like LISP) are substituted by programs of 
imperative language (like C/C++) [3]. Noneffective codes (introns) are identified 
and eliminated. It allows the corresponding useful instructions to be extracted 
from a program, and the process can be accelerated [7]. AIMGP is a particular 
form of LGP. For the given dataset, AIMGP writes a computer program and results 
with the “Best Program” as well as “Best Team.” It creates a total of 30 programs 
in a project and produces five teams of 1,3,5,7,9 program combinations [5]. 


6.5 Case study in India 


In the present case study, one day and two days ahead air pollutant forecasting 
models have been developed for Pune city (18.5218.5204°N, 73.8567°E), which is 
one of the highly polluted metropolitan cities located in the Maharashtra state of 
India. It is the premier industrial center and the fastest developing city of the coun- 
try [37]. The city is located on the western side of the Maharashtra state on the 
Deccan plateau at an elevation of 560 m above mean sea level at Karachi. Rapid 
industrialization along with boom in construction sector has resulted in exponen- 
tial vehicular population growth, which has resulted in increased pollutant concen- 
trations and worsening of the air quality of Pune for the last few years [32]. Central 
Pollution Control Board (CPCB) of India has declared GO. NO,, and RSPM as the 
criteria air pollutants which decide the air quality [35]. The national ambient air 
quality standards of India sets the maximum permissible limit for daily average 
concentrations of SO, NO,, and RSPM as 80, 80, 100 pg/cum, respectively, for resi- 
dential urban areas. For the last decade, it has been observed that the upper limits 
have been crossed for SO,, NO,, and RSPM with a maximum value recorded as high 
as 87,158 and 680 ug/cum, respectively, for the study area [33]. 

Daily values of the pollutant concentrations recorded by Maharashtra pollution 
control board from January 2005 to July 2015 were used for the development of an- 
nual and seasonal models. In order to obtain the complete annual picture about the 
air quality of Pune, the annual models are developed by arranging the yearly data 
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Define parameters 
SO, (t), SO, (t-1), SO, (t-2), SO, (t1) 


Define fitness function 


Create first generation of 
chromosomes 


termination 
criteria 


Figure 6.2: Typical GP Flowchart. 
Source: Yadav & Sharma, 2019. 
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serially [16]. Similarly, pollutant concentrations vary with the environmental condi- 
tions as observed in different seasons. In India, four seasons such as winter 
(Jan-Feb), summer (Mar-May), monsoon (Jun-Sept), and post monsoon (Oct-Dec) 
are observed [24]. Seasonal models are also developed by dividing the data into dif- 
ferent seasons. 


6.5.1 Experimentation using ANN and GP 


Air quality is time-dependent phenomenon. In this chapter, using soft computing 
approaches of ANN and GP, annual and seasonal models have been developed with 
the time series data of the criteria air pollutant concentrations. Uniform seasonal 
patterns and linearity between the inputs were assumed while developing all the 
forecasting models [26]. The data pertaining to the pollutant concentrations (SO;, 
NO,, and RSPM) recorded from January 2005 to July 2015 have been considered for 
model development. The input optimization has been carried out using correlation 
analysis, considering the annual temporal data spanning over the week. Inputs 
were added one by one, and the multiple linear regression (MLR) models were de- 
veloped and evaluated. The correlation coefficient of the MLR models suggested 
that three, four, and three predecessor values shall be used as model inputs for 
SO;, NO,, and RSPM models, respectively. The same number of inputs was used for 
respective seasonal models as well. 
Temporal annual 1 day ahead forecasting models can be written as [29] 


SO»(t +1) =f (SO2(t), SO2(t 1), SO2(t — 2)) (6.1) 
NO,(t +1) = f(NO,(t), NO,(t — 1), NO(t - 2), NO,(t-3)) (6.2) 
RSPM (t 1) = f(RSPM (t), RSPM (t — 1), RSPM (t 2)) (6.3) 


Temporal annual 2 days ahead forecasting models can be written as 


S02(t + 2) =f(SO2(t), SOx(t - 1), SO2(t - 2)) (6.4) 
NO, (t2) -f(NO, (t), NO, (t-1), NO, (t- 2), NO, (t= 3)) (6.5) 
RSPM (t + 2) = f(RSPM (t), RSPM (t- 1), RSPM (t- 2)) (6.6) 


and temporal seasonal 1 day ahead forecasting models can be written as 


SO;(t + 1) =f(SO2(t), SO2(t — 1), SO2(t - 2)) (6) 
7 
for winter/summer/monsoon/post monsoon 


NO, (t+1) 2 f(NO, (t), NO, (t-1), NO, (t- 2), NO, (t- 3)) (6.8) 
8 
for winter/summer/ monsoon/postmonsoon 
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RSPM (t + 1) - f(RSPM (t), RSPM (t — 1), RSPM (t - 2)) 


(6.9) 


for winter/summer/monsoon/postmonsoon 


Similarly, temporal seasonal 2 days ahead forecasting models can be written as 


S02(t + 2) = f(SOx(t), SOx(t — 1), SO2(t - 2)) 


for winter /summer/monsoon/postmonsoon (6.10) 

NO, (t+2)=f(NO, (t), NO, (t— 1), NO, (t- 2), NO, (t-3)) 
for winter/summer/monsoon/postmonsoon (521) 
RSPM (t + 2) =f (RSPM (t), RSPM (t - 1), RSPM (t - 2)) ef 


for winter /summer/monsoon/postmonsoon 


One day and 2 days ahead, pollutants are predicted using ANN and GP with control 
parameters listed in Tables 6.1 and 6.2. 


Table 6.1: Control parameters used in ANN. 


S. no. 


1 


Items 


Network architecture 


Criteria used in this chapter 


Input neurons = 3, 4, and 3 for S02, NO,, and RSPM 
Output neurons = 1 variable for each model 

Hidden neurons = smallest number of neuron which yields 
minimum prediction error on the validation data set 


2 Neuron activation function Input neuron = identity function 
Output neuron = identity function 
Hidden neuron = hyperbolic tangent function logsig and 
purelin for all the models 

3 Learning parameters The learning parameters converge to the network 
configuration and give best performance on validation data 
set with least epochs 

4 Criteria for initialization of Network weights are uniformly distributed in the range of 

the initial weights -1to1 

5 Training algorithm Levenberg Marquardt 

6 Stopping criteria Performance goal/epochs 

7 Performance indicator R, RMSE, d, MBE 
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Table 6.2: Control parameters used in GP. 


S. no. Parameters used Value 

1 Maximum initial tree size 15-20 

2 Maximum tree size 45-50 

3 Population size 50-1,000 
4 No. of children produced 100-3,000 
5 Mutation 0.05 

6 Crossover rate 0.4-1.00 
7 Objective type R?, RMSE 


6.6 Results and discussion 


We have demonstrated the application of GP, which is one of the evolutionary algo- 
rithms, a subset of machine learning to forecast the air quality parameters a few time 
steps in advance. One day and two days ahead, temporal annual and temporal sea- 
sonal models are developed and evaluated using ANN and GP. Owing to the numer- 
ous advantages, ANN has been commonly used for air quality forecasting [30]. It has 
the capability of adaptive learning self-organization and handling of the nonlinear 
systems [20]. We have used ANN as a benchmarking tool in this work. GP is relatively 
a new approach that has been used for air quality forecasting of Pune city. ANN and 
GP models have been developed and tested for unseen inputs, and the qualitative 
and quantitative performance is assessed by the correlation coefficient (r), root mean 
squared error (RMSE), d-statistics (d), and mean bias error (MBE). 


6.6.1 Annual models 

Temporal one day and two days ahead forecasting of air quality parameters has 
been carried out using ANN and GP. 

6.6.1.1 Forecasting of SO; 

Sulfur dioxide is an air pollutant that is responsible for acid rain, haze, and many 


health-related problems. It is produced when the coal is burnt to generate electric- 
ity. As per CPCB report, 2010, the chief source of SO; in Pune city is industry 
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(73.53%) followed by mobile sources (12.70%) and brick kilns (6%). Industrial area 
in Pune has a very high fuel consumption. It is the largest forging industry in Asia. 
Observed concentrations of SO; are slightly higher than the permissible limits. Both 
the soft computing approaches, namely, ANN and GP worked reasonably well as 
evident from the performance indicators, but the performance of GP is slightly bet- 
ter than ANN for the present case study (refer Table 6.3). In the process of predic- 
tions using ANN, the structure is defined initially, and the learning algorithm finds 
the weights. In GP, functions are defined, which always results in an optimal solu- 
tion. One day ahead prediction results are better than two days ahead as it can be 
observed from Fig. 6.3 and Fig. 6.4 as it can be seen with increased r. 

Correlation coefficient of SO; (t+1) for GP is 0.65 and for SO; (t+ 2) is 0.6 and 
d-statistic of SO, (t + 1) for GP 0.8 and for SO, (t + 2) for GP is 0.75. The peak observed 
in testing for one day ahead prediction was 50 pg/cum and predicted by the ANN 
model is 44.35 pg/cum. The peak predicted by GP is 45.60 pg/cum. The observed 
peak for two days ahead forecasting was 46 pg/cum. Peak prediction by ANN is 
33.81 ug/cum, and by GP it was 38.14 ug/cum. 


6.6.1.2 Forecasting of NO, 


NO, is produced in the air from the reaction of nitrogen and oxygen gases, espe- 
cially at high temperatures, during combustion. The amount of nitrogen oxides 
emitted into the atmosphere as air pollutant can be significant, especially in the 
areas of high motor vehicle traffic in Pune city. Increased levels of nitrogen oxides 
can cause damage to the human respiratory tract and increase the risk of respira- 
tory infections and asthma. Chronic lung diseases can be a result of long-term expo- 
sure to high levels of nitrogen dioxides. According to the CPCB report 2010, the 
significant contribution of 9596 of NO, in Pune is from vehicles. The industrial con- 
tribution is limited to 296 as industrial area is confined within the city limits. 
Remaining 396 of NO, is produced from domestic and commercial fuel burning for 
cooking. Modeling and forecasting of NO, can help to warn the residents, especially 
sensitive groups. It was observed that the concentrations of NO, exceeded the an- 
nual standards during study period. 

One day and two days ahead NO, forecasting models have been developed 
and evaluated using ANN and GP. Alike SO,, it can be seen that GP results are 
better than ANN and one day ahead predictions are superior than two days ahead 
(Tables 6.3 and 6.4 and also from Fig. 6.5 and Fig. 6.6). It can be seen that NO, 
forecasts are more accurate than SO; forecasts. As far as errors are concerned, 
RMSE values are higher than as compared to SO; models, and negative MBE is ob- 
served for NO, models indicating that the models are under predicting. 

The peak observed in testing for one day ahead prediction was 115 pg/cum 
and predicted by the ANN model is 53.64 pug/cum. The peak predicted by GP is 
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Table 6.3: One day ahead temporal annual model results. 


Performance evaluator Models Pollutant NO, RSPM 
SO; 

r ANN 0.6673 0.8110 0.8374 
GP 0.6776 0.8258 0.8434 

RMSE ANN 6.1876 11.6846 32.4299 
GP 6.0417 12.2770 31.6747 

d ANN 0.7989 0.8863 0.8956 
GP 0.7995 0.8943 0.9088 

MBE ANN 0.5077 0.2663 -2.1348 
GP 0.2041 -0.0391 -0.5727 

Table 6.4: Two days ahead temporal annual model results. 

Performance evaluator Models Pollutant NO, RSPM 

so, 

r ANN 0.6523 0.8109 0.7006 
GP 0.6487 0.7752 0.7864 

RMSE ANN 6.2901 11.6856 42.3596 
GP 6.2298 12.6171 36.4318 

d ANN 0.7880 0.8862 0.8173 
GP 0.7585 0.8591 0.8713 

MBE ANN 0.4669 -0.2707 -2.1363 
GP 0.3714 -0.1105 -1.2505 


123.45 ug/cum. Similarly, the observed peak for two days ahead forecasting was 
155 pug/cum, and peak predicted by ANN and GP is 80.26 and 83.83 pg/cum, re- 
spectively, which indicates that GP predicts the peak closer to the observed 


values. 
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6.6.1.3 Forecasting of RSPM 


Particulate matters with an aerodynamic diameter less than or equal to 10 um are 
classified as respirable suspended particulates (RSP). Particulates are produced 
from combustion processes, vehicles, and industrial sources. They are microscopi- 
cally small and can enter the lungs and penetrate deeply into the respiratory sys- 
tem and cause difficulty in breathing as well as respiratory and cardiovascular 
effects such as bronchitis and asthma. Sensitive groups, including pregnant 
woman, children, and older people, are at high risk. A warmer climate increases 
the risk of illnesses and death from extreme heat and poor air quality. Pune city 
has a hot and semiarid climate, which further increases the risk [32]. Vehicular 
emission is a major source of particulate matter in Pune and the highest recorded 
concentration of RSPM is 680 pg/cum during the study period. RSPM is the pollut- 
ant responsible for deciding the air quality. Modeling and forecasting of RSPM are 
essential to calculate the air quality index, which can be easily understood by the 
community. 

From the observation table 6.3 and 6.4 and also from Fig. 6.7 and Fig. 6.8, it 
can be seen that similar to SO, and NO, results, one day ahead forecasts of RSPM 
are better than two days ahead, and the GP models performed better than ANN. The 
fact that air pollution is a time-dependent phenomenon is verified from the results. 
There is a strong correlation between previous concentrations of RSPM with the 
forecasted values, which is evident from r values indicating a stronger correlation. 
“r” js a linear measure, and the air quality model evaluations based on only “r” 
mostly fail due to the presence of lag between source emission quantity and the am- 
bient pollutant concentration. Temporal models fail to take into account the lag ef- 
fect as the meteorological parameters are absent as inputs. The lag is observed due 
to inversion condition, which indicates the accumulation of pollutants in the ambi- 
ent environment during odd hours of the day when there are no source emissions 
recorded [13]. In such conditions for air quality modeling, “d” statistics is the rele- 
vant evaluation criteria. *d" values of 0.90 and 0.87 are recorded for RSPM one day 
and two days ahead forecasting, which indicates that about 9096 and 8796 model 
results are error free. Negative MBE demonstrates that both the tools, ANN and GP, 
underpredict the phenomenon. Lower RMSE values (31.67 ug/cum for one day and 
36.43 ug/cum for two days ahead forecasting) in the data values ranging between 17 
and 680 ug/cum indicate the satisfactory performance of GP. 

Alike SO, and NO, models, peak predicted by GP for RSPM one day ahead 
model is closer to the actual value (observed peak is 291 ug/cum and the peak pre- 
dicted by ANN is 242.55 pg/cum and by GP is 244.12 ug/cum), whereas the observed 
testing peak for two days ahead prediction was 680, but the peak predicted by ANN 
is 145.36 ug/cum and by GP is 153.28 ug/cum. This could be due to the fact that 
RSPM concertation of 680 pg/cum was observed only once during Diwali 
(festival) day. 
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6.6.2 Seasonal models 


Seasonality is one of the essential factors affecting air quality. Variations in the pol- 
lutant concentrations are observed in different seasons due to the interaction of the 
various meteorological parameters with pollutant concentrations [27]. Seasonal 
models are formulated and evaluated using ANN and GP to get the complete picture 
of the air quality of Pune. 

In Pune city, winter is observed during January and February. Air quality gets 
worsened during this time due to inversion condition. In winter time, a temperature 
inversion takes place when a layer of warmer air traps cold air near the ground. 
During inversion, the air becomes stagnant and pollutants are caught close to the 
land, which may result in the formation of smog. 

Summer starts at the beginning of March and ends at the end of May. As com- 
pared to winter, summer heat prevents the inversion, which improves the air qual- 
ity slightly. Summer always seems to drive out the dense clouds of pollution that 
suffocate many Indian citizens. Air pollutant concentrations are observed to be 
higher during the summer time because the heat and sunlight primarily react with 
the air along with all the chemical compounds in it. When these compounds react 
with the nitrogen oxide emissions present in the air, it results in smog. According to 
medical studies, air pollution in summer increases the risk of heart attack [31]. 

June to September is declared as monsoon months in Pune [26]. Rainfall is a 
cleansing mechanism that results in washing away of the pollutants during mon- 
soon [26]. Rain makes the ground wet and mixing of the dust into the air again be- 
comes difficult resulting in lesser pollutant concentration. Ozone formation rate is 
also lower in monsoon as the rain clouds obstruct the correct amount of tempera- 
ture and sunlight required for ozone formation. 

October and November are the two months of post monsoon (or retreating mon- 
soon) in which the wind becomes weaker. Skies are clear, but humidity in the atmo- 
sphere increases. Diwali festival happens to be in this season. The air quality in this 
season decreases due to firecrackers. This condition increases the chance of pollu- 
tion episodes. 

All four seasons have different characteristics and from the observation tables, 
it can be seen that both the tools (ANN and GP) worked well, but GP results are 
better than ANN. 

From the observation Tables 6.5—6.7, it can be seen that one day ahead predic- 
tions are superior to two days ahead forecasts. Results indicate that SO; models have 
shown better performance in summer and post-monsoon season. NO, models have 
performed well in monsoon and post-monsoon and all RSPM models performed very 
well in all the seasons except for winter. The probable reason could be in winter 
there are slight movements of the air masses resulting in stagnation of the pollutants. 

It was observed that maximum concentrations of SO, NO,, and RSPM were re- 
corded in winter, whereas the minimum concentrations have been recorded in 
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Table 6.5: Seasonal forecasting of S03, 


Pollutant Time of prediction Season Tool Performance evaluator 
r RMSE d MBE 
SO; 1 day Winter ANN 0.4734 8.4711 0.6688 0.7368 
ahead 
GP 0.5577 7.5759 0.6963 0.3866 
Summer ANN 0.6544 4.8672 0.7408 —-0.0325 
GP 0.6546 4.7015 0.8267 0.3179 
Monsoon ANN 0.5692 5.1757 0.7408 —-0.0325 
GP 0.7067 4.3710 0.8267 0.8267 
Post ANN 0.6338 7.8940 0.7865 0.6965 
monsoon 
GP 0.7066 6.8280 0.8219 0.3703 
2 days ahead Winter ANN 0.3024 9.0665 0.5142 1.0553 
GP 0.3588 8.6694 0.5503 —-0.0517 
Summer ANN 0.6022 5.4609 0.6044 2.0816 
GP 0.5914 5.0841 0.7040 0.9720 
Monsoon ANN 0.5556 5.3268 0.7079 0.0609 
GP 0.5734 5.0645 0.7919 0.1971 
Post ANN 0.4560 8.7618 0.6392 0.4464 
monsoon 
GP 0.4582 8.6370 0.6161 0.4400 


monsoon as well as winter. Heavy fog during early mornings of winter supports the 
highest concentrations during winter. Rainfall was observed during study time in 
winter for few days resulting in the minimum pollutant concentrations. Rainfall acts 
as a cleansing agent and results in washing away of the pollutants in monsoon. From 
the meteorological records, it can be seen that the quantity of rainfall was sufficient 
to reduce the concentrations of the contaminants in monsoon. 


6.7 Modeling options for situations when there are 
constraints on the data availability 


Quality and consistency of the data decides the model prediction accuracy. 
Availability of continuous and accurate meteorological and pollutant concentration 
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Table 6.6: Seasonal forecasting of NO,. 


Pollutant Time of prediction Season Tool Performance evaluator 
r RMSE d MBE 
NO, 1 day Winter ANN 0.6495 17.2756 0.7647 -1.1737 
ahead 


GP 0.6846 16.4143 0.7849 -1.4931 


Summer ANN 0.6758 19.3564 0.8681 -0.4552 


GP 0.6920 18.7333 0.8810 -0.1921 


Monsoon ANN 0.7221 8.1978 0.8163 0.5864 


GP 0.7784 7.4318 0.8687 0.3330 


Post ANN 0.7248 13.1202 0.8074 -1.7946 
monsoon 


GP 0.8080 11.0645 0.8848 0.0178 


2 days ahead Winter ANN 0.5625 18.6040 0.6885 -1.5606 


GP 0.6511 17.0322 0.7722 -0.6801 


Summer ANN 0.5860 8.9691 0.7328 1.1892 


GP 0.5946 8.7939 0.7301 0.7584 


Monsoon ANN 0.6352 9.4085 0.7155 2.0252 


GP 0.6775 8.7317 0.7903 0.5660 


Post ANN  Á 0.6245 14.9824 0.7135 -1.4086 
monsoon 


GP 0.63843 14.4177 0.7526 -0.0676 


data is a challenging task in developing countries. Adverse meteorological condi- 
tions, instrumental errors, and so on are also some of the reasons for nonavailabil- 
ity of continuous and accurate data of all the meteorological parameters and 
pollutant concentrations. Atmospheric inversion is a vital condition deciding the 
air quality, which must be taken into account while modeling the air quality pa- 
rameters. Data of meteorological parameters such as wind speed, rainfall, solar 
radiation, wind direction, humidity, mixing height, visibility, and color cover 
should be included as model inputs. Processing time for data collection and dis- 
semination is quite high; as a result, it becomes difficult to include meteorology 
as a model input as in the case of this work. In such situations, the models devel- 
oped should be robust enough so that they can be useful for short-term prediction 
of criteria air pollutants. An attempt has been made to suggest various air quality 
modeling alternatives. 
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Table 6.7: Seasonal forecasting of RSPM. 


Pollutant Time of prediction Season Tool Performance evaluator 
r RMSE d MBE 
RSPM 1 day Winter ANN 0.5244 30.8937 0.6563 9.4319 
ahead 
GP 0.5553 28.8443 0.6839 3.1098 
Summer ANN 0.7935 19.3564 0.8681 —0.4552 
GP 0.8850 18.7333 0.8810 -0.1921 
Monsoon ANN 0.7806 19.3564 0.8681 -0.4552 
GP 0.7964 18.7333 0.8810 -0.1921 
Post ANN 0.6675 34.4305 0.7947 -0.5861 
monsoon 
GP 0.6826 33.6118 0.7914 -0.3259 
2 days ahead Winter ANN 0.3941 34.4683 0.8810 11.6724 
GP 0.3801 32.5242 0.8824 6.3909 
Summer ANN 0.7729 52.3487 0.7881 -9.3464 
GP 0.7926 47.1344 0.8596 -3.0583 
Monsoon ANN 0.6942 22.2994 0.8072 -0.5772 
GP 0.7104 21.7756 0.8122 -0.4177 
Post ANN 0.4932 81.7014 0.4789 -31.2984 
monsoon 
GP 0.5527 14.4177 0.6815 -0.02559 


When meteorological parameters are available, the correlation of inputs (meteo- 
rological parameters) with the output (pollutant concentrations) can be calculated 
and data can be arranged in the descending order of their correlation with the output. 
Only three inputs (from the available inputs) with high correlation with the output 
shall be considered and called as top three causes. Those having low correlation with 
the output are called as bottom three causes. The modeling options can be presented 


as [29] 


Model type 1 - model with top three inputs 


Modeltype2 - model with bottom three inputs 


(6.13) 
(6.14) 


Model type 3 - model with bottom three and toponeinput (6.15) 


Modeltype4 - temporal models 


(6.16) 
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Of the various meteorological parameters listed above, the temperature difference 
is the crucial factor for atmospheric inversion and, at the same time, it is effortless 
to record and disseminate. Temperature model can be developed in which previous 
values of the temperature difference can be added on the basis of correlation along 
with previous pollutant concentrations in case of the temporal model, and such 
modified temporal model can be presented as 


Model type 5 - temperaturemodels (6.17) 


With these modeling options, criteria air pollutants can be predicted on a continu- 
ous basis. 


6.8 Conclusion and future directions 


We have presented the application of soft computing techniques to forecast air 
quality parameters one and two days in advance. Models are developed with GP, 
and the results are compared with established forecasting tool of ANN. Temporal 
annual and seasonal models are developed for all the three criteria air pollutants 
(SO;, NO,, and RSPM) and found that model results are matching with the recorded 
values. One day ahead forecasting results are better than two days ahead. GP 
proved to be better than ANN in terms of accuracy of forecast and found equally 
competent for peak predictions for this case study. 

In air quality modeling, researcher always prefers a short-term and continuous 
model to broadcast the current status of the air quality. This information is useful 
to plan the activities and safeguard the health of sensitive groups. In developing 
countries like India, there is a problem of availability of continuous data. In such 
situations, various modeling options are proposed. 

The major limitation of the study is the absence of meteorological parame- 
ters for model development. Model inputs were selected using a linear method 
of correlation for the air quality forecasting problem, which is actually a non- 
linear phenomenon. There were certain assumptions made while developing 
forecasting models such as seasonal patterns were assumed to be uniform, 
meteorological processes were simplified, and linearity between the inputs was 
assumed which is seldom true. The studies in future could be directed toward 
the better understanding of the tool GP and combining it with other air quality 
forecasting approaches. Comprehensive models can be prepared by considering 
the vertical and spatial distribution of the pollutants along with meteorological 
parameters. 
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Palash Dutta, Bornali Saikia 

7 Arithmetic operations on generalized 
semielliptic intuitionistic fuzzy numbers 
and their application in multicriteria 
decision making 


Abstract: This chapter attempts to initiate a novel generalized semielliptic intuition- 
istic fuzzy number (GSEIFN) along with basic arithmetic operations on GSEIFNs. 
Furthermore, an advanced ranking method for GSEIFN is proposed. In the end, a 
multicriteria decision-making problem has been carried out by using the proposed 
GSEIFN and the ranking approach to exhibit the legality and applicability. 


7.1 Introduction 


Multicriteria decision making (MCDM) is a common activity in our daily life, such as 
choosing a car, selecting a candidate in an interview for the post of a teacher, and 
choosing or buying the best book or the other things from online shopping apps. It is 
an action or a procedure of choosing the best option from the accessible options 
through multiple criteria. In making decision, due to the presence of conflicting criteria 
or vague information, it is impossible to make proper decision by using a traditional 
set theory. To deal with such uncertain or vague data, Zadeh [1] introduced the fuzzy 
set theory (FST), which deals with those elements that partially belonged to the set. 
Later, Atanassov [2] extended FST, which is known as intuitionistic fuzzy set (IFS), 
which deals with both elements that are partially belonged to the set and not belonged 
to the set. IFS deals with such situations in an enhanced manner in comparison to FST. 
Due to growing complexities in making decisions, several intuitionistic fuzzy numbers 
(IFNs) have been developed [3-5]. More often, for simplicity, triangular or trapezoidal 
IFNs (TrIFNs) are employed to solve decision-making problems. However, in some real- 
world situations, the interactions between available system factors are found to be 
more complicated and in such situations triangular or trapezoidal types of IFNs fail to 
provide proper result. To make a proper decision, it is more essential to reconstruct an 
appropriate IFN. To overcome such difficulties, generalized semielliptic (GSE) type of 
IFN may be developed by making more suitable to handle such complex environments. 

In MCDM problems, ranking of IFNs plays an efficient role. In literature, a few 
numbers of ranking approaches are encountered based on value and ambiguity [6], 
centroid point [7], centroid [3], magnitude [8], area [9], score and accuracy function 
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[10], possibility mean [11], mean value [12], and improved possibility degree [13]. 
The approaches [3, 4, 6, 14-24] focused on the value as well as the ambiguity-based 
ranking method of IFNs and applied in MCDM problems. The approaches in [7] and 
[3] dealt with the centroid-based ranking method, the approach in [21] presented a 
ranking technique for IFNs by using the distance of each IFN from the fuzzy origin, 
and the approach in [9] developed a ranking method based on (a, 8)-cuts and area of 
IFNs. Notwithstanding having some primacy of ranking approaches, these approaches 
have some deficiencies and do not provide a proper result, which produces predica- 
ment for judgment architects in any decision-making problem. To overcome such lim- 
itations, a new ranking approach has been developed based on the consummated 
exponential mean. Furthermore, a comparative analysis of the present ranking ap- 
proach and other available techniques has been presented to exhibit the validity of 
the same. In this chapter, an effort has been made to introduce generalized semiellip- 
tic IFN (GSEIFN) and perform its fundamental operations of GSEIFNs. Finally, an 
MCDM problem is being solved to demonstrate the existence and significance of pro- 
posed IFNs and proposed ranking approach. 


7.2 Preliminaries 


Definition 7.1 ([2). An IFS B defined on the universal set X is denoted by 
B= (G6 Up(X), va(x)),x € X}. ug:X — [0,1] and vg:X — [0, 1] are respective membership 
function (MF) and nonmembership function (NMF) such that 0 < ug(x) + vg(x) x 1. The 
amount ztg(x) = 1— Hell — vg(x) is called the degree of hesitancy. 


Definition 7.2 (elt, An a-cut of B is defined as 

"B= (xjug(x) >a} 
For IFN a-cut produces a closed interval, say ^B = [Lp(a), Rg(a)]. 
Definition 7.3 (elt, An B-cut of B is defined as 

B= {x\va(x) <B} 
For an IFN the a-cut produces a closed interval, say PB = [Lg(B), Reif, 
Definition 7.4 (elt An (a, )-cut of B is defined as 

(B= (xlug(x) >a, va(x) « B) 


such that O<a<1,0<fB<1, andO<a+f<1. 
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Definition 7.5 ([18]). The values of IFN B of MF and NMF are denoted by V,,(B) 
and V,(B), respectively, as 


VB) = [aca + Ry(a)]ada 


1 
V.(B)= | itat) + muta - ap 
n 
Definition 7.6 ([18]). The ambiguities of IFN B of MF and NMF are denoted by 
A,,(B) and A,(B), respectively, as 


Definition 7.7. The mean of IFN B are denoted M,,(B) and M,(B), respectively, as 


M,(B) = d [Ls(a) + Re(a)]da 


1 1 
wus = 5| (Lo(B) + Ro(B)dB 
n 
Definition 7.8 ([5]). The general expression of an ellipse centered at (p,q) is 


(x-pY | (y-ay 


fee s 


Assuming that g=0 and k - 1 to construct a normal semielliptic IFN (NSEIFN) P, we 
have 


Then the needed MF is 


2 
up(x) 2 4/1- 5 P) »p-h<x<pth 


In the same way, the NMF can also be constructed as 


(x- p) 
hi 


vp(x) =1-4/1 p-hzxszp«h 
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where h< h; and p is the core of the NSEIFN P, while h and h; indicate support/do- 
main MF and NMF, respectively. 


Definition 7.9. The GSEIFN P = Sg(p, hy, hz;w1, H) on universal set of discourse X is 
defined by the MF and NMF, respectively, as 


2 
up(x) zm DH sp-hs<xspth 
1 
NOS 
vp(x) z1- (1-1) 4/1 "P „p-ħ<x<p+h 
2 


Definition 7.10. Consider P = Se(p, hy, ho; w1, N4) be a GSEIFN. Then (a, f)-cut of the 


GSEIEN is 
v on {E} enie] 
fp- p-r 1- {es Lac: EN 


7.3 Elementary operations on GSEIFNs 


Here h; x ho. 


= 


Here, we use (a, B)-cut method to perform some arithmetic operations of GSEIFNs 
and also discuss some numerical examples. Suppose that P= Sg(p, hy, baam, N4) and 
Q= Sg(q, ki, ko; Wo, ];) are two GSEIFNs whose MF and NMFs are defined as follows: 


2 
x- 
Up(X) = Wi 1- í JP sp-hs<x<pt+th 
1 
2 
vp(x) -1- (1-7) 4/1 oP sp-hn<x<pt+h 
2 
Here hı x hy 
(x- ail 


Holzlzmwaull- B ,q- hxxxq*k 
1 
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vox) =1-(1-1)4/1 E q-k<x<qtk 
2 


Here k x k. 
Now we find out that the (a, B)-cuts of P and Q are 


rene Diech 

ve pn fe En ET] 
e e H7 | 
ta- bn en nmn] | 


Here we also use w; = min(w;, w2) and 7, = max(n,, n) to perform arithmetic opera- 
tions. 


N 


7.3.1 Addition of GSEIFNs 


For the evaluation of sum of GSEIFNs P and Q, it is required to sum the a-cuts and 
B-cuts of P and Q via interval operations and here we assume w; = min(wi, w2), 


n, = DANI, na): 


2 2 
reer0-|ipoa)-roknft-{2} wrom {E | (7.1) 
2 2 
teata loso med HÄ risch BH (3) 


Now for getting proper MF yp, 9(x) and NMF vp, (x), it is important to equate ini- 
tial component and final component of eqs. (7.1) and (7.2) with x which give 


ug tg (er ta fi [EP] 
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-p10 turk) (SY pras tere ft (108 


Next, it is needed to express a in terms of x by taking a2 0, a < wn and ax wi, az 0 
in eq. (7.1) we have 


We get the domain of x, 


€ [(p-h) + (q- kı), (p* a). x € [(p +q), (p h) + (a ki)] 


In the same way, f can be expressed in terms of x by taking B<7,, 820 and £20, 
B xn, in eq. (7.2) we have 


Si 2 x 2 
B=1-(1-n,)4/1 eo ,B=1-(1-,)\/1 eet 


We get the domain of x, 


€ [(p- h2) + (a- k2), (P+ a). x € [(p* q), (p + h2) + (q + bi 
Thus, the MF pp, g(x) and NMF vp. q(x) of P+ Q are 


EE oce Iph) + (a-k), (pfi) (a) 


veso=1-0-m)y1- e ipi) «(ak (P+) + (a+) 


7.3.2 Subtraction of GSEIFNs 


For evaluation of subtraction of GSEIFNs P and Q, here, it is also required to sub- 
tract the corresponding a-cuts and f-cuts of P and Q via interval operations and 
considering wn = min(w;, w2), 7, = max(n;; n2), we have 


"P-*Q- (v-a q) - (hi +kı)4j1- (p-q) * (hi a)i]1- n (7.3) 
2 2 
hp-ig- — o-a) doe [T | (r4) 
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Now, for getting proper MF up. g(x) and NMF vp. och of P-Q, it is essential to 
equate both initial and final components with x eqs. (7.3) and (7.4) which give 


x-(p- 4) - (hi +ki) i (E x= (p-a) -atk ien) 


Next, it is further essential to express a in terms of x by considering a2 0, axwi 
and a x w, a > O in eq. (7.3) we have 


e EE 
Aa=wi4/1 { Ee E Wull { ake 


We get the domain of x, 


x € (p- hı) - (q+ kı), (p- a).x € [(p - a). (p 3) - (a- ki)] 


In the same fashion, f can be expressed in terms of x by taking f xr, B>0 and 
B 20, B xn, in eq. (7.4) we have 


B=1-(1-m)4/1 (re OV pat (1-m)4/1 pea 


We get the domain of x, 


x € [(p - h2) - (q+ k2), (p - q)];x € [(p - q), (p + hz) - (q - k2)] 
Thus, the MF p. (x) and NMF vp. o(x) of P - Q are 


ne arm xD eech (a+ ky) (ph) - (7) 


= _ 2 
vyo-q=1-(1-m)yf1- PGP EON xe ipt) - el o) - (ak) 


7.3.3 Product of GSEIFNs 


For evaluating the product of GSEIFNs P and Q, it is required to product the a-cuts 
and f-cuts of P and Q via interval operations: 
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hat Etat 


DEEN (7.5) 
cl BECH D HEN i 


(pn GEF) (eE) 


To find MFyupo(x) and NMF vpo(x), here too it is needed to equate first and last com- 
ponents to x in eqs. (7.5) and (7.6), which gives 


rH e9 ET) 
exce (on) 
"eH EID T 
REH 


Putting w; = min(w;, w2) and expressing a in terms of x 


bb (- EN - (pla + qhi) | 1- CH +pq-x=0 
hk, (- ER + (pla + qh) 1- ny +pq-x=0 


2 
Next, solving the expression by taking 1— (a) - y^, we have, 
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- (- (pla + qh))) + y {- (pki +qħ)¥ - 4hikı(pq- x) 
2hik, 


2 


(pl; + qhi) - "TS + qh) - 4hika(pq - x) 
&=wı,]|1 
2hikı 


Similarly, 


2 


(pk; + qhi) - "TS + qh) - 4hiks(pq - x) 
&=wı,]|1 
2hik, 


Setting a>0, a < w; and a< w, a> 0 in eq. (7.5), we have 


x € [p- hi)(a - ki). pq].x € [pa (p + hi) (q+ ki)] 


Putting 7, =max(7,,7,) and expressing f in terms of x 


mifi- EM - (pk +qħ2)4/1- [EE ag 


mifi- EM + (pk; * qh;)4]1— [EBV an 


2 
Again solving the expression by taking 1— EA =q’, we have 


- (- (pk qho)) + y/{- (blo + aha) Y - 4lle(pq - x) 
2hjk, 


q= 


2 


B=1-(1-m),/1 E ke + ho) - y eme 
2K2 


Similarly 


2 


1, f fed -y (pk + gha)? -4hklpa -x) 
B Ei ( i) 2hjko 
Setting £ «1, B2 0 and f 2 0, B <n; in eq. (7.6) we have 


x € [(p- hi)(a - ka), pq]; x € [pa. (p + h)(a + k2)] 
Thus, the MF jp; (x) and NMF vpo(x) of PQ are 
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2 


Hpo(X) =w 


(pk, + qhi) - "TS * gh)’ - 4hyk (pq - x) 
2hik, 


(p-h)(q-k)sxs(p*hi))(q-*k) 


2 


Vpo(X) 21- (1—11;) 


(pk; + qh;) - V Oks + qhi) - 4hik;(pq - x) 
2h;k; 


(p- h)(q-Io) <x< (p h»)(q + Io) 


7.3.4 Division of GSEIFNs 


For evaluation division of GSEIFNs P by Q, it is essential to divide the a-cuts by £- 
cuts of P and Q via interval operations: 


y y 
p-h 1-{ {a} pth, 1-{ {a} 


S (7.7) 
Eu. {a} q-k fay’ 
2 y 
af SEHE us 


Bp p-m 
BQ” 232 
ge kd n- {EEV o-bdi- (52) 
Now, for getting exact MF ps (X) and NMF va(x) of P/Q, it is required to equate first 
and last components to x in eqs. (7.7) and (7.8), we have 


hyy/1- bal, rrhdi: fal’ 
EU {a} -fa 


N 
" 
DE 
—— 
N 


KSC cuf (E 


Next, express a in terms of x by considering w; = min Dun, w2), «420, axw; and 
axw;,az20 in eq. (7.7) we have 
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2 2 
_ |, fp-x | A f p-xq 
ROME 1 pa > MA ux 


We get the domain of x 


e 4 » E CA 
qtk’ ol q q-k 


In the same fashion, f can be expressed in terms of x by taking 7, =max(n,,75), 
<n, ß z0 and f 20, ß <n in eq. (7.8) we have 


B-i-ü nay JE) GE TEEN 


We get the domain of x 
x KS x ES 
q+k q q q-k 
Thus, the MF p(x) and NMF vp(x) of 5 are 
q Q 


2 
8 p-xq p-h pth 
npe) -ws ; Ud ? qi q-k 


2 
= p-xq p-h pth 
ui Gay erac] P wks qe 


It should be noted that fundamental operations on GSEIFNs produces GSEIFN. 


7.3.5 Illustrative examples 


Here, numerical examples of GSEIFNs are being discussed with the help of above 
arithmetic operations. 

Suppose that P= Sg(10, 4, 5; 0.2, 0.5) and Q=Sz(8, 2, 3; 0.5, 0.5) are GSEIFNs and 
their respective MFs up(x), Holz) and NMFs vp(x), vo(x) are as follows: 


x-10 


Hp(x) =0.2 1-( ) x € [6,14] 
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U(x) = 0.54/1- ey x € [6,10] 


2 


2) , X € [5,11] 


3 


vo(x) 21- (1- 0.5/1 ( 


The (a, B)-cuts of P and Q are 


7.3.5.1 Sum of GSEIFNs 


Sum of GSEIFNs P and Q is given by P+ Q whose MF yp, 5 (x) and NMF vp, g(x) are 


2 
Hp, g(x) =0.4/1- K =) , x € [12,24] 


2 
Vp. o(X) 21- 0.5 e) , x € [10,26] 


which shows this is again a GSEIFN. 


7.3.5.2 Subtraction of GSEIFNs 


Subtraction of GSEIFNs P and Q is given by P- Q where MF Ur o(X) and NMF 
Vp- Q(X) are 
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2 
Hip. o (X) =0.24/1 (52) , X € [-4,8] 


2 
Vp-o(X) -1- 0.54]1 (53). x [- 6,10] 


which shows that subtraction of two GSEIFNs P and Q is again a GSEIFN. 


7.3.5.3 Product of GSEIFNs 


Product of GSEIFNs P and Q is given by PQ where MF ji; (x) and NMF vpo(x) are 


52 - 4/ (52)? - 32(80 - x) 


20.,1|1 
Hpo(X) 16 > 


36<x<140 


2 


70 - 4/ (70)? - 60(80 - x) 


Vpo(X) —-1-0.,1!1 30 » 


25<x<165 


which shows that multiplication of two GSEIFNs P and Q is again a GSEIFN. 


7.3.5.4 Division of GSEIFNs 


Division of GSEIFNs P and Q is given by 6 whose MF u p (x) and NMF vp (x) are 


which shows that division of two GSEIFNs P and Q is again a GSEIFN. 
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7.4 Mean and value of GSEIFN 


Here, an attempt has been made to study mean and value of the newly defined 
GSEIFN by taking (a, B)-cuts. Suppose that P = Sein, hy, hz;w1, lisa GSEIFN and its 
MF as well as NMF are as presented in Section 7.2. The mean of GSEIFN is evaluated 
from Definition 7.7: 


=pwı 


Similarly for NMF M, (P) =p(1-n,). 
The value of GSEIFN P is calculated as 
w1 


V,(P) = j [Lp(a) + Rp(a)]ada 


= pw; 


Similarly for NMF V,(P) = p(1-7,)’. 


7.5 Comparing the proposed IFN with the other 
existing IFNs 


In this section, we compare the introduced IFN with the other existing IFNs by 
using some numerical examples and existing ranking approaches. Mainly we com- 
pare GSEIFN with the normal trapezoidal (NTrIFN), normal triangular IFN (NTIFN), 
and NSEIFN with the help Rezvani [17] and Z. Xu’s [25] ranking method as discussed 
below. 


7.5.1 Ranking approach of Rezvani 
Rezvani [17] proposed the ranking approach for TrIFNs based on the value. We can 


easily evaluate the rank of two different TrIFNs by using this method. By this ap- 
proach, values for the TrIFN P= < qj, 2,05, 04:41 05, A3, d'4 > are evaluated as 
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_ a +2a2 + 203 t Q4 

d 6 

For two different TrIFNs A and B, rank can be evaluated by using the following 
procedure: 


WA) 


If V (A) > V,(B) then A» B 
If V,(A) < V,(B) then A <B 
This approach gives the counterintuitive output in some situations. For example, 


for the NTrIFNs A= < (4, 5.5, 6, 8);1,0> and B= < (3.5,5,7,7.5);1,0> values of A and 
B denoted by V,A and V,,B are evaluated as 


2x5.54+2 
Viet "7 X6*8 033 
542 2 ; 
vje a X7+73 83 


This gives the same output for non identical NTrIFNs A and B. It is problematic for 
the decision makers to identify the best one. To get rid of this situation, NTrIFN can 
be reconstructed in the form of GSEIFN as given below. 

Suppose that A= Sz(6, 2, 2;0.5,0.5) and B = Sg(5.5, 2, 2; 0.1, 0.7) are two GSEIFNs 
whose MFs and NMFs are 


and 


Then a-cuts and f-cuts of A and B are 


2 2 
hea |G so hee P uut 
a- le-ai- (e ni (8 
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BR ke air Don BEA 


Through Definition 7.5, the values for A and B can be measured as 


V,(A) = IR z Dr IO) en SÉ ny ada - 15 
V,(B) = i (55-2 j1- fay £55421/1- EN oda 0.055 


which shows that V,,(A) > V,(B) and this indicates that A > B. It is witness that draw- 
backs of earlier method can be overcome. 


7.5.2 Ranking approach of Xu 


Using score and accuracy function, Xu [25] developed this ranking approach [25]. 
Suppose that A= < a1, a2, a3;W1, Uu and B= «bi bj Daa, Hz are two triangular 
IFNs. Then by Zus ranking method, score functions denoted by S(A), SCH) and ac- 
curacy functions denoted by H(A), H(B) are as follows: 


S(A) = wn -N S(B) = W2- i, H(A) = wi n H(B) = w2 +N, 


Then the rank of A and B can be determined by using the following procedure: 

If S(A) < S(B) then A < B 

If S(A) = S(B) then we check the value of accuracy function 

If H(A) > H(B) then A» B 

If H(A) < H(B) then A <B 

For the two NTIFNs A= <(4,5.5,6);1,0> and B= < (3.5,5,7);1,0>, using Xu's 
approach the score and accuracy values are obtained as 


S(A) =1, S(B) = 1, H(A) -LH(B) -1 


which provides the identical result. 
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Similarly, NSEIFNs A = < (4,2, 3);1,0> and B= <(8,3,5);1,0> whose MFs and 
NMFs are 


and 


For these IFNs, value score and accuracy functions are 


S(A) =1, S(B) = 1, H(A) =1,H(B) = 


These IFNs also provide identical results. So it is difficult for the experts to distin- 
guish the best alternative by using NSEIFN. So we reconstruct the IFN in the form of 
GSE type. Suppose that A = Sz(4, 2, 3; 0.1, 0.05) and B= Sz(8, 3, 5; 0.2, 0.03) whose MF 
and NMFs are 


and 


=1-0.974/1- € [3,13] 


Then the values of score functions and accuracy functions are 
S(A) = 0.05, S(B) = 0.17, H(A) = 0.15, H(B) = 0.23 


which shows that A < B. From this comparative analysis, it can be opined that the pres- 
ent approach has the ability to overcome those drawbacks of NTIFNs and NSEIFNs. 
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7.6 Ranking approach of IFN based on consummated 
exponential mean 
Here, a novel ranking approach of IFNs based on the consummated exponential 


mean has been presented. 
Let B be an IFN whose (a, B)-cuts are considered as follows: 


"B= [Lg(a), Ra(a)), ^B- Leif, Re(B)] 


Now, evaluate the mean of B as 


for normal w=1 and 7 =0. 
Now the consummated exponential mean of B is given by M(B) and is defined as 


M(B) = eM (B) + gMy(B) 


Suppose that P and Q are two IFNs, then the ranking procedure is as follows: 
If M(P) » M(Q) then P» Q 
If M(P) < M(Q) then P<Q 


Consummated exponential mean of generalized trapezoidal IFN (GTrIFN) 
P= < [py p2, P3, Pa; wi], Ip, P2; P3» P43] > is defined by 


M(P) - eU mz +Pa)} 4 esa «va «ps p) 


Consummated exponential mean of generalized triangular IFN (GTIFN) P< 
[Di Pa Dai Wi], Ip, P2; Pit! > is defined by 


M(P) ln: «p4)) 4 ef y sen) 
Consummated exponential mean of GSEIFN P (p, hy, hw; 2) is defined by 


M(P) = e?" + e-m) 


7.6.1 Numerical illustration 


Here, a discussion has been made on some numerical examples with the help of 
proposed ranking approach. Suppose that P= < (1,2,3, 4; 0.2), (0.5, 2, 3, 4.5; 0.4) > 
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are two GTrIFNs. Then the consummated exponential mean of P and Q are evalu- 
ated as follows: 


M(P) eld 24344] y e{ 0524344) 


= eo + el^ 


= 6.130411 


Similarly, M(Q) = 3.799609, which shows that M(P) > MIO), that is the rank of order 
P»Q. 

Suppose that P= «(1,2,3;0.2),(0.5,2,3;0.4) > are two GTIFNs. Then the con- 
summated exponential mean of P and Q are evaluated as follows: 


M(P)- eof 220+2+3)} 4 oly (05+2+3)} 


W 


=e A + e 


= 4.572042 


Similarly, M(Q) = 2.645911, 

which shows that M(P) » M(Q), that is, the rank of order P» Q. 

Suppose that P= Sz(5, 2,3; 0.2, 0.3) and P= Sg(4,1, 2; 0.1,0.3) are two GSEIFNs, 
then the consummated exponential mean of P and Q are 


M(P) - 4.137349, M(Q) - 3.717366 


which shows that M(P) > M(Q), that is, the rank of order P > Q. 


7.6.2 Comparative study of the present approach with the 
available approaches 


Here, we compare the present ranking approach with the other existing methods in 
literature as shown in Table 7.1. 

From the above table it is observed that the ranking approaches developed in 
[6, 10, 15, 17, 26, 27] as well as in [9] produce the same result for nonidentical IFNs. To 
overcome these drawbacks and to get a proper result, we should develop a proper 
ranking approach. In this chapter, we introduced a consummated exponential mean- 
based ranking approach and from the comparative Table 7.1, it is encountered that the 
present ranking approach has the capability to overcome the limitations of the avail- 
able methods. 
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Table 7.1: Comparative table. 


IFNs Available Present 
approaches approach 


P= <(4,5.5,6);1,0>, Q= < (3.5,5,7);1,0> [25], P- Q P>Q 


P= < (0.45, 0.63, 0.83);0.73,0.2>, Q= < (0.48, 0.64, 0.82); 0.72, 0.2» [10], P=Q P<Q 


P= ((0, 0.25, 0.3);1,0), Q= ((0.1, 0.2, 0.4);1, 0) [15], P- Q P«Q 

< (0, 0.25, 0.3);1,0>, Q= < (0.1, 0.2, 0.4);1, 0> [6], P- Q P«Q 
P= < (0.20, 0.35, 0.40, 0.55);1,0>, Q= < (0.10, 0.40, 0.45, 0.45);1, 0 > [26],P=Q  P»Q 
P= <(4,5.5,6,8);1,0>, Q= <(3.5,5,7,7.5);1,0> [17], P=Q P>Q 
P= < (2.57, 2.73, 2.83);0.73, 0.2» , Q= < (2.58, 2.74, 2.82);0.72, 0.2 > [27],P=Q P>Q 
P= <(3.5,5,7,7.5);0.5,0.5>, Q= <(4,5.5, 6, 8);0.7, 0.3 > [7], P>Q P<Q 
P= < (0.20, 0.35, 0.40, 0.55);1,0>, Q= < (0.10, 0.40, 0.45, 0.45);1, 0 > [28], P<Q P>Q 
P= < (4,4, 4);1,0>, Q= < (2,2, 2);1,0 > [9], P=Q P>Q 


7.7 MCDM via GSEIFN using the present ranking 
approach 


Here, a methodology of MCDM problem based on the proposed ranking method and 
proposed IFN has been discussed. 

Suppose that P = (P,, P5, ..., P4) be a set of m alternatives and N = (Ni, N>, ..., Nm} 
be a set of n criteria for each alternative P;. Let wo = (Wa; We, Wem) be the weight 
of criteria where We, € [0,1] and we, + Wc, + --- * wc, =1 to choose the best alterna- 
tive among Pys for i-1,2,..,n based on the consummated exponential mean 
method. The methodology is discussed as follows: 


Step 1: Form a table of fuzzy rating criteria by choosing fuzzy criteria in terms of 
GSEIFNs ` po = (pi; Wy, Ni), (i21,2,3,.,m;j-1,2,3,.,n) and py eR and 
weight of criteria We. 

Step 2: A decision matrix B= (Pij)mxn has been formed by using the above choos- 
ing table of fuzzy rating criteria: 


N, hi P; 2 Pr 

N> (Pu; Wii» Ny) (py W12, Ny) wee (Pins ; Win; Nin) 
B= | (Pa; Wa, Na) (Pas Wna) = (Dani Wans Non) 

Nn 


(Pmi Wm Nm) (Bis Wma Nma) -~ (Binns Wns Nmn) 
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Step 3: If the given data are inadequate, then normalize the choosing decision ma- 


trix B= (by), by using the equation 
. WaqDij 
dij z 100 (7.9) 


Then the normalized decision matrix will be of the form 


" P, P, Cé D 
N> (aui wio ni) (ous Wi» Mo) = (ami Wim Min) 
R= (an; wa, ia (4225 Wo, Mol eM (Gn; Won» Non) 
Now Lier Wan Mal Ier Wha» Hl  — (Qian; Wa Mn) 


Step 4: After getting the normalized decision matrix R = (aj) 
using the following equation: 


evaluate R(P;) by 


mxn? 


R(P;) = Y MO. i=1,2,...,n (7.10) 
k=1 


where M(P;) Ng 1S the consummated exponential mean-based rank of GPIFNs. 


Step 5: The decision makers choose the alternative of maximum R(P;). 


7.7.1 Numerical example of MCDM 


In this section, an MCDM problem has been discussed to exhibit the validity and 
applicability of the present ranking method and proposed IFN. Suppose that a stu- 
dent searching for a good college among available colleges, namely, D. P2, P3 for 
him located in a city on the basis of four criteria represented by GSEIFNs: 

(1) N;- Transportation 

(2) N5- Yearly cost 

(3) N3= Teaching quality 

(4) N,= Yearly performance 


Let the weight of criteria w, = (0.1,0.2,0.3, 0.4) and the fuzzy rating of criteria of 
each alternative as given in Table 7.2. 


152 —— Palash Dutta, Bornali Saikia 


The decision matrix B = (bj) is formed by using Table 7.2: 


Table 7.2: Rating of criteria. 


P, P; Ps 
N: (40, 20, 30; 0.2, 0.3) (90, 40, 50; 0.4, 0.5) (50, 30, 40; 0.3, 0.4) 
N2 (30,10, 20; 0.1, 0.2) (40, 20, 30; 0.2, 0.3) (60, 40, 50; 0.2, 0.5) 
N3 (50, 30, 40; 0.3, 0.4) (80, 30, 40; 0.4, 0.5) (90, 40, 50; 0.4, 0.5) 
N, (60, 40, 50; 0.2, 0.5) (70, 20, 50; 0.3, 0.5) (90, 40, 50; 0.4, 0.5) 


Table 7.3: Table of ranking result by the proposed method. 


Alternatives of 
decision-making problem 


Results by using 
proposed method 


Py 8.5352206 


P2 8.7712329 


P3 8.9213067 


(40, 20, 30; 0.2,0.3) (90, 40, 50; 0.4, 0.5) (50, 30, 40; 0.3, 0.4) 
(30,10, 20; 0.1,0.2) (40, 20, 30; 0.2,0.3) (60, 40, 50; 0.2, 0.5) 
~ | (50,30, 40; 0.3,0.4) (80,30, 40; 0.4, 0.5) (90, 40, 50; 0.4, 0.5) 
(60, 40, 50; 0.2,0.5) (70, 20, 50; 0.3, 0.5) (90, 40, 50; 0.4, 0.5) 


Now normalized the above matrix by using eq. (7.9) as follows: 


0.04, 0.02, 0.03; 0.0002, 0.0003) 
0.06, 0.02, 0.04; 0.0002, 0.0004) 
(0.15, 0.09, 0.12; 0.0009, 0.0012) 
0.24, 0.16, 0.2; 0.0008, 0.002) 


0.09, 0.04, 0.05; 0.0004, 0.0005) 
0.08, 0.04, 0.06; 0.0004, 0.0006) 
0.24, 0.09, 0.12; 0.0012, 0.0015) 
0.28, 0.08, 0.2; 0.0012, 0.002) 


R= 


Da Mc 


0.05, 0.03, 0.04; 0.0003, 0.0004) 
0.12, 0.08, 0.1; 0.004, 0.001) 
0.27, 0.12, 0.15; 0.0012, 0.0015) 
0.36, 0.16, 0.2; 0.0016, 0.002) 


By using eq. (7.10), result of the decision-making problem is obtained. 
By observing Table 7.3 we have P; > P; > P3; therefore, P; college will be chosen 
by the student. 
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7.7.2 Conclusion 


GSEIFN has been introduced in this chapter. Here we performed some arithmetic oper- 
ations via (a,f)-cut technique and it is exhibited that fundamental operations on 
GSEIFNs also give GSEIFN. There are several types of IFNs that have been developed 
by many researchers. In most of the MCDM problems, researchers used the data or 
information that are in terms of trapezoidal and triangular IFNs. But in some cases, 
these IFNs give incorrect solutions by using NTrIFN, NTIFN, and NSEIFN. In this chap- 
ter, we have shown that the proposed IFN overcomes these limitations in NTrIFN, 
NTIFN, and NSEIFN, although a new ranking method has been developed by using 
the exponential mean of MF and NMF and comparing this method with the other ex- 
isting methods with the help of some existing ranking approaches and numerical ex- 
amples. Finally, an MCDM problem has been done by using the proposed IFN and the 
present novel ranking method to showcase the validity and applicability. 
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Laxminarayan Sahoo 
8 Method for solving intuitionistic fuzzy 
assignment problem 


Abstract: The method to find an answer for assignment problem (AP) under intui- 
tionistic fuzzy domain is proposed in this chapter. Due to the irregular rising and 
falling of the present market economy, here we have assumed that the assignment 
costs are not always fixed. Therefore, the assignment costs are imprecise in nature. 
In the existing literature, different approaches have been used, which are interval, 
fuzzy, stochastic, and fuzzy-stochastic approaches to represent the impreciseness. 
In this chapter, we have represented impreciseness taking intuitionistic fuzzy num- 
bers (IFN). The proposed method is hinged on ranking of IFN and use of well- 
known Hungarian method. Here, we have used a newly proposed centroid concept 
ranking method for IFNs. In this chapter, we have solved AP where costs for assign- 
ment are taken as triangular IFNs. A numerical example has been considered to de- 
rive the optimal result and also to adorn the applicability of the suggested method. 
In the end, concluding remarks and future research of the proposed approach have 
been presented. 


Keywords: assignment problem, intuitionistic fuzzy numbers, Hungarian method, 
impreciseness, triangular intuitionistic fuzzy numbers 


8.1 Introduction 


Assignment problem (AP) is effectively used in solving realistic problems for finding 
the best alternative solutions. AP plays an essential role in managerial decision- 
making. It is related to the problem of production planning, scheduling, and engi- 
neering design problem. The basic task of AP is to allocate many tasks/jobs in a more 
efficient way such that ideal assignment can be obtained in an acceptable limit. In an 
AP, n workers and n jobs and effectiveness of each worker for each job, the problem 
involves in assigning each worker to only one job so that the total amount of effec- 
tiveness is optimized. In reality, different procedures and algorithms, namely, linear 
programming method [1-3], Hungarian algorithm [4], neural network method [5], ge- 
netic algorithm [6], and so on were used to solve the AP. In traditional AP, it is as- 
sumed that all the parameters related to AP are fixed valued. Notwithstanding in 
reality, these parameters are not fixed due to uncertainty and lack of proper informa- 
tion. The uncertainty may occur resulting from some factors: (i) decision maker has 
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no sense about the assignment parameters when a job is to be scheduled at the be- 
ginning time, as a result some uncertainty may occur with regard to assignment pa- 
rameters; (ii) nowadays present market economy is always flooded, that is, always 
unstable due to competing market and, as a consequence, the assignment cost is to- 
tally changeable. To succeed in dealing with these types of condition, the problem 
can be framed using the notion of uncertainty and parameters are served as impre- 
cise/uncertain in nature. In such possible situations, fuzzy set (FS) theory plays an 
important role in picking up such situation. The theory of FS was originated by 
Zadeh [7] in 1965 and it dealt with imprecision and vagueness in real-world situa- 
tions. In the year 1970, Belmann and Zadeh [8] introduced the notion of decision- 
making problems assuming uncertainty. Here, we have served imprecise parameters 
considering fuzzy numbers. Therefore, the fuzzy AP (FAP) yields an efficient frame- 
work that solves real-life problems with uncertain information. In the last few years, 
several efforts have been made in the existing relevant works for solving FAP. A num- 
ber of researchers [9-15] have formulated AP in different sectors/areas such as firm 
management, vehicle routing, and network management considering fuzzy parame- 
ters and/or interval parameters. For solving these problems, they have also used 
several mathematical techniques/methods like linear programming technique and 
genetic algorithm. Chen and Chen [16] described about fuzzy optimal solution of 
the AP based on the ranking of generalized fuzzy numbers. Fuzzy risk analysis 
based on ranking fuzzy numbers using a-cuts was proposed by Chen and Wang. 
Rommelfanger [17] discussed a ranking of fuzzy cost present in the FAP, which takes 
more advantages over the existing fuzzy ranking methods. Kar et al. [18] unfolded the 
fuzzy generalized AP with restriction on available cost. Solution of FAP with ranking 
of generalized trapezoidal fuzzy numbers was resolved by Thangavelu et al. [19]. 
Andal et al. [20] discussed a technique to minimize cost flow of FAP considering 
fuzzy membership. Pramanik and Biswas [21] introduced multiobjective AP with gen- 
eralized trapezoidal fuzzy numbers. Mukherjee and Basu [22] suggested the applica- 
tion of fuzzy ranking method for solving AP with fuzzy cost. Very recently, Sahoo 
and Ghosh [23] solved AP with linguistic costs. 

In fuzzy optimization technique, the degrees of receiving objective functions and 
constraints are considered. FS theory has also been evolved in several fields and its 
different modifications/upgradations have come to light. The most important general- 
izations of FS theory is an intuitionistic FS (IFS) which was originated by Atanassov 
[24] in 1986. The concept of IFS is an alternative approach to define FS in the case 
where information is not available sufficiently for defining impreciseness by means of 
ordinary FS. An IFS is characterized by the degree of acceptance and degree of rejec- 
tion. The IFS theory has been applied in different fields including decision-making. 
For more details, one may refer to the works of Jana and Roy [25], Roy et al. [26], and 
Kumar and Hussian [27]. 

In this chapter, we have considered intuitionistic FAP (IFAP). To solve the IFAP, 
we have used a new ranking method based on centroid concept introduced by Arun 
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Prakash et al. [28]. Then the AP has been remodeled into an AP whose parameters 
are fixed/crisp valued and solved by linear programming method and/or Hungarian 
method [4]. An example has been studied to obtain the optimal result and also to 
illustrate the applicability of the proposed method. The rest of the chapter is pre- 
sented as follows: 

In Section 8.2, brief notions about IFSs and fuzzy numbers are described. 
Mathematical formulation of AP is given in Section 8.3. Section 8.4 gives the solu- 
tion methodology of the IFAP. In Section 8.5, a numerical example is given for illus- 
tration purpose. Section 8.6 concludes the chapter with further scope of research. 


8.2 Preliminaries 


In this section, a brief idea about IFSs and fuzzy numbers used in this chapter is 
presented. 


Definition 8.1 ([23]). Let X be a given nonempty set. An IFS A in X is defined by 
A= (Gua 0). ya(0):x € X} where ju; (x): X — [0,1] and y; (x): X — [0,1] define, re- 
spectively, the degree of acceptance and degree of rejection of the element x € X to 
the set A which is a subset of X and every x € X, 0 < u; (x) + y; (x) <1. Again each IFS 
A in X, the function 7z;(x) -1- u; (x) - y;(x) is called the degree of hesitation. If 
7t; (x) =0 for all x € X then IFS becomes FS. 


Definition 8.2. An IFS A is called normal if there exists xo € X such that ju; (xo) 2 1 
and y; (xo) - 0. 


Definition 8.3. An IFS A is called convex if for all xix; € X and O<A< m 
(Axı + (17 2)x2) 2 Min {pug (1), y3 (X2) } and y; (oa + (1 Ae) > Min y 63), ya 09) J- 


Definition 8.4. Support of IFS A with universal set X is denoted by Support (A) 
and is defined by Support (A) = {x:u;(x) > 0 and yj(x) «1, x e X}. 


Definition 8.5. An IFS A= { (x, Ha (x), yq(x)):x € X}of the real line is called an IFN if 
(i) Ais intuitionistic fuzzy (IF) normal; 

(ii) A is IF convex; 

(iii) ua) is upper semicontinuous and y; (x) is lower semicontinuous; and 


(iv) Support(A) is bounded. 
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Definition 8.6. An IFN A = (a), a2.a3)(ā, a2, 3) where à; < a4 < a < a3 < ā; is said to 
be triangular IFN if its membership and nonmembership functions, respectively, 
are defined as follows: 


x-a D 
SEN ifa <x<@ 
E - d3-X i 
Hy (x) au, ee 
0 otherwise 
Q5)-X “cs 
and TEN ifa,<x<a@ 
» Ms x-a d W 
ya(X) = EET ifa; xxxds 
1 otherwise 


Definition 8.7 ([27]). The centroid point (x(A), y(A)) of triangular IFN A = (a1, a2.a3) 
(44, a5, a3) can be determined by 


dı + A2 03 
“A= — 
201 — A2 + Q3 
a fmt 
1 
ENER 
y,(4) 3 
and j 
Ai 
y,(A) 3 


Definition 8.8 ([27]). The ranking function of the triangular IFN A is defined by 
1 2 1 2\2 
DEENEN 


Theorem 8.1. Let A and B be two triangular IFNs then p(A+C)>p(B+C) > 
A+C>B+C. 


m 
+ 

8 
V 

f 
+ 

D 

t 

box 
ma 
+ 


Now if, A> B then o 


Theorem 8.2. Let A and B be two triangular IFNs, then, 
A<B if p(A) «p(B) 


A» B if p(A) > p(B) 
A=B if p(A) = p(B) 
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8.3 Mathematical formulation of assignment 
problem 


Suppose there are n jobs to be performed and n workers/machines are available for 

doing these jobs. Here we have assumed that each worker/machine can do each job 

at a time, though with swing degree of efficiency. Let cj be the cost or time if the ith 

person/machine is assigned the jth job, the problem is to find an assignment so 

that the entire cost of payment or time for performing all jobs are minimum. 
Mathematically, AP can be modeled as follows: 


n n 
Minimize z= 3 3" cyxy, (5j21,2,.., n) 
i=1j=1 


1 if ith worker/machine is assigned jth job (8.1) 


Subject to x= 
^ l 0 if ith worker /machine is not assigned jth job 


and 


where xj; denotes that jth job is to be assigned to the ith worker/machine. 

In most of the cases, it is observed that the cost of payment/time of assignment 
is not precise due to the fluctuation of market economy and, hence, we have as- 
sumed that the assignment costs/times are imprecise and have also considered that 
they are fuzzy numbers Cj(i-1(1)nand j=1(1)n). Then AP (8.1) becomes FAP, 
which can be written as follows: 


n n 
Minimize Z- 35^ Gage, (Lj -1,2,... n) 
fent 


1 if ith worker/machine is assigned jth job (8.2) 


Subject to xj = 
i l 0 if ith worker/machine is not assigned jth job 


and 


where xj; denotes that jth job is to be assigned to the ith person/machine. 
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Theorem s 3. If xj-Xj minimizes z= xx vé Do Over all xj - 0 or 1 such that 
i-1jz 

Y EAR Y Am =1 then xj = Xj also minimizes z= Y vé e mn where c č ij = Či - Ui- Vj 

j=1 j=1 ni^ 1j=1 


Proof. Let, xi = Xj minimizes the total cost z= 3 3 i-i yxy over all xj - 0 or 1 and 
isi 


Za zw 


Now it is to be shown that x; = X;; also minimizes a new cost Z- » 3 (Cy — ŭi - Vj) 
i=1j=1 
Xij where u; and v; are imprecise cost constants subtracted from cost matrix (Cj) ER 


3 


n n 


So, here 2’ is independent of vue, and hence Z is minimized when Z is mini- 
mized and vice versa. 


Theorem 8.4. In an AP with cost matrix (Cij) if all čij 2 0 then a feasible solution 


ENEE nxn? 
xj is optimal when 3" Y; Go - 0. 
i-1j-1 


~ n n ~ 
Proof. Since all cj > 0 and all xj > 0, therefore (-- 2 3 Gm | 2 0. Hence, the min- 
i-1jz1 ~u 
imum possible value of z- 2 2 Ca can attain is always 0. Then any feasible so- 
ait 1j- 


lution x; is optimal iff 5 Sg Cjxj - 0. 
erte 


8.3.3 The Hungarian method 


The following method applies Theorem 8.4 to a given nx ncost matrix to obtain an 

optimal assignment [4]. 

Step1: For each row i,i=1,2,...,n, let ci, r € (5,2,.., n] be the smallest entry and 
calculate cj - Cir, foreach iandr € {1,2,...,n}. 
Obviously Min{ cj - Cir} =0, foreachiandr € {1,2,...,n}. 

Step 2: For each column j,j=1,2,....n, let cy,k € (,2,.., n] be the smallest entry 
and calculate cj — cj, foreachjandk € {1,2,...,n}. 
Obviously Min{ cj - cj} =0, foreachjandk € {1,2,...,n}. 

Step 3: Find minimum number of lines l (say) in such a way that l 1 ines covered 
all the zero entries of the cost matrix. 

Step 4: Now, either /=n or l<n. If [- n, then an optimal assignment is arrived. (ii) 
If 1<n then go to Step 5. 

Step 5: Find the smallest entry not covered by any line, subtract this from each 
uncovered row and column, and then add it to junction entry and go to 
Step 3. 
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8.4 Solution methodology 


Suggested approach of solution procedure of linguistic AP is as given below: 


Step 1: 
Step 2: 


Step 3: 


Step 4: 


Find p(¢j) of each Ci; using ranking function defined in Definition 8.8. 
Check whether the given AP is balanced or not. 

(i) Ifthe AP is a balanced then go to Step 4. 

(ii) If the AP is unbalanced then go to Step 3. 

Make dummy rows and/or columns with zero cost so as to form a balanced 
AP. 

Solved by method A: (existing linear programming problem method) and/ 
or by method B: (the Hungarian method described in Section 8.3.3) to ob- 
tain the appropriate assignments. 


n n 
n İ=1j=1 
subject to 3 sl, j=1,2,...,n 
i=1 


n 
RE: ENEE 
j=l 


Xij € {0,1} 
Method B: Find an optimal assignment using Hungarian algorithm. 


8.5 Numerical example 


In this section, we have considered an IFAP with three resources Wi, W2, W3 and 
three jobs Jı, J2, J3 with IF assignment cost. The cost matrix (Cil) s is given below 
whose elements are triangular IFNs. 


Activities 


Resource 


A h h 
W, (4,7,14),7,24) (12,14,16)(11,14,28) (8,15,25)(7,15,30) 
W; (6,18,25)(5,18,28) (13,19,28)(12,19,36) (20,25,33)(19,25,38) 
Ws; (20,25,38)(19,25,42) (2,11,14)(1,11,19) (6,12,15)(5,12,20) 


Solution: Here, it is seen that the IFAP is a balanced one. Now, we have obtained 
R(cy) of each cjusing ranking function mentioned in Definition 8.8. Now the re- 
duced problem is as follows: 
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Activities 
A h h 
Resource Wi 12.02 15.57 17.41 
W; 15.67 22.49 27.38 
W3 29.58 8.33 11.35 


Then, the AP can be modeled in the following linear programming problem as 
follows: 


Minimize z = 12.02xy, + 15.67x2; + 29.58x3, + 15.571? + 22.49X2 + 8.33X32 (8.3) 
8.3 
+ 17.4143 + 27.38X23 + 11.35x33 


Subject to 
Xu + Xa + X31, 71; Xi + Xz t X32 = 1; Xi3 + X73 t X33 = 1; 
Xu +X + Xy3 = 1; Xa +X2 + X93 = 1; X31 + X32 + X33 = 1; 
xj € {0,1}, i=1,2,3, 4andj=1,2,3 


Solving problem (8.3) we get the optimal assignments W; — J, W; — Jj, W — Jo. 
Also proceeding by the Hungarian method, the optimal assignments are Wi — Jo, 
W, — Jı, W3 — J3. Now, for optimal assignments W; — J3, W2 — J4,W; — Jz the mini- 
mum assignment cost in terms of IF is (8,15,25)(7,15,30) + (6,18,25)(5,18,28) + (2,11,14) 
(1,11,19), that is, (16,44,64)(13, 44, 77) and for optimal assignments Wi > Jo, 
W,— Jı, W3 — Jz the minimum assignment cost in terms of IF is (12,14,16) 
(11,14,28) + (6,18,25)(5,18,28) + (6,12,15)(5,12,20), that is, (24, 44, 56)(21, 44, 76). 
More effectively, in reality the assignment cost cannot be taken as fixed/precise one. 
So, in this chapter we have solved an AP considering each cost of the assignment is 
IF. For this purpose, to solve the AP with IF cost we have used the ranking function 
method defined in Definition 8.8 For the IF optimal costs (16, 44, 64)(13, 44, 77), we 
have obtained p((16, 44, 64) (13, 44, 77)) 2 31.95 and p((24, 44, 56) (21, 44, 76)) = 33.59. 
Using, Theorem 8.2, we have seen that p((16,44,64)(13, 44, 77)) « p((24, 44, 56) 
(21, 44, 76)); hence, (16,44,64)(13, 44, 77) is the minimum assignment cost. 


8.6 Conclusion 


In this chapter, AP with IF assignment cost/time is considered and its solution pro- 
cedure has been established. AP is one of the most essential problems in decision- 
making. In everyday situations, parameters of the AP are imprecise. The AP with IF 
cost is more reasonable than the AP with precise parameters because most of the 
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real-life events are from uncertain domain. The proposed procedure presented here 
is very effortless and easy for implementation. This solution procedure is based on 
IF representation of assignment parameters and ranking measure of IF parameters. 
Here, a newly developed ranking measure based on centroid concept has been 
used. Then, the AP has been transformed to an AP whose parameters are crisp/pre- 
cise valued and solved by the existing Hungarian method/linear programming 
method. Lastly, an IFAP has been solved, and evaluated results have been con- 
ferred and measured. It may be asserted that the entire solution procedure set out 
in this chapter can be put in to solve other realistic decision-making problems in 
the near future demanding IF parameters. 
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9 Optimization of EDM process through 
evolutionary computing and fuzzy MCDM 
techniques 


Abstract: This chapter intends to explore the best possible set of process parameters, 
namely, current, voltage, and pulse on time during electric discharge machining 
(EDM) process in order to examine the operational changes pertaining to the rate of 
material removed and the rate with which electrode wears. Here the workpiece 
selected is high carbon high chromium die steel, and the titanium nitride—coated 
copper electrode is used. Initially, the data set is trained and validated and finally 
tested by employing artificial neural network. Then multiobjective genetic algorithm 
(M-GA) is applied in order to find better parametric combinations for the response 
values. Finally, multicriteria decision making has been used in order to obtain the 
noblest parametric combination from the data of M-GA, where two contradictory re- 
sponses like high material removal rate and low electrode wear rate can be achieved. 


Keywords: HCHCr tool, MRR, EWR, ANN, GA, fuzzy, MCDM 


9.1 Introduction 


As there is advancement in the usage of day-to-day materials, the secondary 
manufacturing processes have become very essential in order to sustain in that 
challenging manufacturing industries. Electrical discharge machining (EDM) is ex- 
tensively popular in the nonconventional machining process, while manufacturing 
various industrial critical components related to plastic molding, automobile, aero- 
space, and so on. The process utilizes thermal energy in the form of controlled 
sparks to remove materials from the conductive workpiece. Due to this distinctive 
feature of the process, it provides an edge in machining complex contours on high- 
temperature strength-resistant materials over their contemporary traditional coun- 
terparts [1]. In the EDM process, series of sparks are generated from the conductive 
electrode due to the necessary potential difference created between the narrow tool 
work gaps. The complete machining zone is flooded with die electric fluid, which 
simultaneously cools the system as well as helps in flushing the debris produced 
during machining. The temperature generated during this spark emission across 
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the tool work gap ranges between 8,000 and 12,000 °C, which is capable enough in 
melting all types of conductive electrodes that are difficult to machine materials. 
The parameters that influence the machining process most are duty cycle, pulse off 
time, pulse on time, current, voltage, spark gap type of electrode, and so on. 

Few past research works on EDM have been elaborated here in order to select the 
proper design of experiment and control parameter. Kanthababu [2] has introduced 
the basic ideas of single and multiobjective optimization, and also the evolutionary 
algorithms, as well as gives a brief outline of their relevant characteristics. Some of 
the metaheuristic algorithm normally applied while solving critical multiobjective 
functions are shown explicitly. Palanikumar et al. [3] have discussed the optimization 
of performance characteristics by using the Taguchi method along with an evolution- 
ary gray-fuzzy system with respect to several responses like material removal rate 
(MRR), surface roughness (Ra), and specific cutting pressure (Pr). Their results show 
that the optimization technique is very significant while performing machining oper- 
ations of Glass Fiber Reinforced Polymer (GFRP) composites to achieve good surface 
quality with minimum low tool wear concurrently. Pandian Vasant [4] during solving 
of fuzzy programming problems in industrial production systems has operated with 
three metaheuristic optimization techniques, where a decision maker takes the final 
call. Here in order to unravel the comprehensive challenges of genetic algorithms 
(GA) practically, mesh adaptive direct search methods are taken into consideration 
along with the pattern search. Meyer et al. [5] have investigated problems through 
distinctive market research with respect to small data sets through statistical techni- 
ques. With these small data sets, they have added three different neural network 
models to this investigation. Deo et al. [6] have predicted the wind speed of 15 differ- 
ent places in Queensland, Australia, by employing artificial neural network (ANN) 
along with GA. Sahu and Mahapatra [7] have undergone EDM on titanium alloy work- 
piece. Three different tool materials, namely AlSiMg, copper, and graphite, are used 
while altering different control parameters. Finally, the surface finish, namely, Ra, 
Rt, and Rz is determined. Bose et al. [8] have performed optimization of EDM process 
using artificial bee colony algorithm for multiresponse optimization and whale opti- 
mization algorithm in case of single response optimization. Bose and Pain [9] have 
focused on the optimization of electrochemical grinding process, while machining 
Alumina-aluminium (Al,03—Al), an interpenetrating phase composite (IPC) by 
using bioinspired metaheuristics. Prabhu and Kumar [10] have investigated the 
electrochemical deburring method using nickel-coated tool on AlgOg, material 
successfully for deburring complex geometrical cross-holes. The input parameters 
such as current, electrolyte concentration, current density, and time varied in 
order to evaluate the performance of the tool as well as the analysis of machining 
characteristics are done. 

From the brief literature review, the primary goal of this chapter is to explore the 
best possible set of process parameters, namely, current (J), voltage (V), and pulse on 
time (ton) in order to determine the changes in MRR and electrode wear rate (EWR) 
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during EDM of high carbon high chromium die steel (HCHCr) workpiece by applying 
the titanium nitride-coated copper electrode. Applying the regression analysis, the 
initial investigation is commenced. Then the data set is trained and validated and fi- 
nally tested applying ANN. Then multiobjective GA (M-GA) is applied to those data in 
order to find a better input-output parametric combination. Finally, to find out the 
noble parametric combination to satisfy the two contradictory responses like high 
MRR and low EWR, multicriteria decision making (MCDM) is used. 


9.2 Experimental design 


The experimental investigation is done using die sinking type EDM (ACTSPARK SP1, 
China make). A positive polarity is maintained while using titanium nitride-coated 
copper tool electrode. The dielectric used is EDM oil 30 having 0.80 specific gravity 
(at temp. 25 °C) and 3.11 cSt. viscosity (at 38 °C). The flushing pressure of 0.2 kgf/cm? 
is maintained in between the copper tool (12 x 12 mm?) and the workpiece during ma- 
chining. The tool work gap is made constant by using a servo control system. 

The experimental setup is shown in Figure 9.1. 


Figure 9.1: Machining setup of EDM process. 


168 —— Goutam Kumar Bose, Pritam Pain 


Keeping in view of the earlier research works done in this setup of previous experi- 
mentation, three control parameters are chosen. These parameters having three lev- 
els of variation are considered for 17 experimental runs. The input parameters with 
their levels of variation are illustrated in Table 9.1. Some of the insignificant factors 
retained are fixed. 


Table 9.1: Process controller with different levels of variation. 


Process variables Notation Units Levels of variation 

Ly L; Ls 
Current I amp 6 8 10 
Voltage V volt 45 50 55 


Pulse on time ton US 50 100 200 


The MRR and EWR are calculated as the weight of the workpiece material removed 
over time and are expressed as g/min, respectively. 


9.3 Regression analysis 


Generally, in case of one dependable variable, let us say “y” which is dependent on 
an independent variable “k.” The correlation among these two types of variables is 
characterized as a regression model. A strong correspondence exists relating to the 
design of experiments and regression analysis. A regression model is developed by 
using a subset of the existing regression which involves contradictory objectives as 
follows: 
— The regression model should incorporate maximum regressors for composing 
the resultant equation beneficial for the analysis. 
— There should be minimum regressors in the regression model in order to lower 
the cost involved in gathering and sustaining the information comprehensively. 


MRR = - ( - 0.016 + 0.0264*I — 0.00021*V + 0.000074*ton — 0.000164*I*V 
+ 0.000007*I*ton — 0.000004*V*ton) (9.1) 


EWR = - 0.0564 + 0.00675*I + 0.000867*V + 0.000041*ton — 0.000075*T*V 
— 0.000005*T*ton + 0.000001* V*ton (9.2) 
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9.4 Artificial neural network 


An ANN is a bioinspired metaheuristic. A neural network consists of nodes that are 
interconnected thoroughly. Each neuron sends an input signal to the adjacent neu- 
ron. It then subsequently exudes the output signals (O;, i=1 to j) to the subsequent 
neuron. Every input P; multiplies by the corresponding weight W; in synapses and 
then transfer to the hidden layer as follows: 


n 
= $ WP (9.3) 
i=1 


Output from every node will generate only when the input reaches the threshold 
level of that node. This output is a sigmoidal function in general as follows 


0; =f (Ii) =1/(1+e ul (9.4) 


Figure 9.2 shows the architecture of a node of ANN. 


P1 
Net input 

P2 i 0, 

= È W, P; 

t=1 
Net output 0j 
=0 

Pha 

Pa 


Figure 9.2: Illustration of a synapse of artificial neural network. 

The commonly used neural network architectures are as follows: 
— Self-organizing networks 
— Feedforward networks 


— Backpropagation networks 


The diagrammatic representation of an ANN has been illustrated in Figure 9.3. 


9.4.1 Result analysis applying ANN 


The design of experiment is conducted by using soft computing technique through 
Matlab R2015a. Several algorithms are considered to facilitate the calculation such 
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Input layer 

Output layer 
Input layer 
Output layer 


Hidden layer Hidden layer 


Feedforward Feedback 


Figure 9.3: Schematic design of a neural network. 


as divider and function for dividing the data set, for training Levenberg-Marquardt, 
and mean squared error (MSE) for performance study. In Figure 9.4, a simulink 
model is illustrated having the four parameters under consideration in order to op- 
timize the different responses. 


Neural Network 
Hidden Output 


Algorithms 
Data Division: Random (dividerand) 
Training: Levenberg-Marquardt (trainlm) 


Performance: Mean Squared Error (mse) 
Calculations: MATLAB 


Figure 9.4: Simulink diagram. 


The gradient as well as the Jacobian is determined by using backpropagation algo- 
rithm techniques. 

By applying the chain rule of calculus in the reverse order, the backpropagation 
computation is achieved in the overall network. Here in the ANN model MSE is for- 
mulated between the output of the network and the corresponding target value. 
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An epoch represents the set of initial value and/or target value (training vectors) in 
an ANN network which is multiplied by its corresponding weights. Initially 1,000 
epochs, which is the maximum number possible, are set to stop the training for all. 
The graphical representation of the regression displays the correlation among the tar- 
gets and network outputs. The regression plots illustrate the training, validation, test- 
ing, and subsequent overall data. The target, which is the difference between the 
perfect result and outputs, is illustrated by dashed lines. Whereas the best fit line is 
illustrated as a continuous line, which signifies the linear regression line relating to 
the targets and outputs. When a direct correlation exists between the targets and the 
outputs ten R becomes one, else it becomes zero. 


9.4.1.1 MRR analysis 


Now while computing using ANN the following data distribution is made. For train- 
ing, validation, and testing 60%, 15%, and 25% data are taken into consideration, 
respectively. Here MSE is achieved after three number of iterations. The training will 
terminate when the value of the gradient is less than 1.00 x 10" and when the valida- 
tion check number is 6. During computation, the gradient observed is 1.5 x 10 ? and 
the validation checks reach 4 after successive iterations. 

Figure 9.5 illustrates the execution of neural network training for MRR. 


Progress 

Performance: 0.00835 |- 3.8 SE | 0. 00 
Gradient: 0.0185 31. 00e-07 
Mu: mc 1.00e+10 


Validation Checks: 


Figure 9.5: Neural network training for MRR. 


It is observed that the best validation result is 0.00068473 at epoch 1. The valida- 
tion performance graph for MRR is shown in Figure 9.6. 

It is shown in Figure 9.6 that at 1, the validation experiment is minimum while 
training continued till 4 before terminating. The plot shows a significant relation- 
ship between testing and validation. The test curves and the validation curves indi- 
cate close resemblance among themselves. 

Figure 9.7 demonstrates the following results of R: training = 0.96358, valida- 
tion = 1, testing = 0.98834, and overall = 0.86907. Hence, it can be inferred that the 
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Best validation performance is 0.00068473 at epoch 1 


Ee —a Train 


— Validation 


— Test 


107? Best 


10720 


Mean squared error (mse) 


1075 


0 0.5 1 1.5 2 2.5 3 3.5 4 
4 Epochs 


Figure 9.6: Validation performance graph for output MRR. 


training data fits good here and both the results of validation and testing are more 
than 0.9. 


9.4.1.2 EWR analysis 


At this point while computing using ANN, the following data distribution is made. 
For training, validation, and testing the data under consideration are 60%, 20%, 
and 20%, respectively. Here the MSE is achieved after the first iteration. The train- 
ing will terminate when the value of the gradient is less than 1.00 x 107 and when 
the validation check number reaches 6. During computation gradient observed is 
2.68 x 10 ? and the validation checks reaches 1 after successive iterations. 

Training for EWR has been illustrated in Figure 9.8. 

It is observed that at the fourth epoch, the finest validation result is attained, 
that is, 1.1425 x 10 ?. Figure 9.9 shows the validation performance for output EWR. 

Figure 9.9 shows that at 4 the validation experiment is minimum while training 
continued till 5 before terminating. The plot highlights a considerable correlation 
between testing and validation. The test curves and the validation curves show 
close similarity between them. 
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Figure 9.7: Regression graph of output MRR. 


Figure 9.10 exhibits the following results of R: training = 0.99986, validation = 1, 
testing = 1, and overall = 0.97029. Hence, it can be inferred that the training data fits 
good here and both the results of validation and testing are more than 0.9. 


9.5 Genetic algorithm 


Depending on the Darwin’s theory of survival of the fittest, a metaheuristics approach 
like GA has been established [11]. It is a nature-inspired metaheuristic that helps 
in optimizing multifaceted design problems, which incorporates a combination of 
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Progress 

Epoch: 0 1000 
Time: 0:00:00 

Performance: 0.000108 0.00 
Gradient: 0.000289 1.00e-07 
Mu: 0.00100 1.00e+10 
Validation Checks: of 6 


Figure 9.8: Neural network training for the output EWR. 


Best validation performance is 1.1425e-05 at epoch 4 
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Figure 9.9: Performance graph for the output EWR. 


constant discrete variables, irregular and lying in nonconvex design space. It is well 
competent in locating the optimum solution globally with maximum probability. 

The first changes on the original bounded function are done by means of pen- 
alty function and transformed it into unbounded function as follows: 


P - variable 


Minimize f(P)+R y D(gi(P)) 
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Figure 9.10: Regression plot of EWR. 
Subjected to p, « p, « py”, x21,2,.., n (9.7) 


where ® denotes the penalty function as follows: 
o-(X)- (xy (9.8) 


where 


eI if oa 35 
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R denotes the penalty constant for the parameters. 
Boosting the fitness function, F(X), the subsequent transformation function is 
minimalized as follows: 


F(P) = Fmax g (re) BEN = Finax -f'(P) (9.10) 
i=1 
where 
Fmax >f (P) (9.11) 


f(P) is the fitness function. 


This fitness function usually performs either YES or NO function. This arrange- 
ment can be looked upon as a genetic chromosome. 


9.5.1 Result analysis applying GA 


Here dual optimization is performed with the help of the regression equations (9.1) 
and (9.2), respectively. The response functions, that is, maximization of MRR and 
minimization of EWR are ambiguous in nature. But they are to be simultaneously 
optimized satisfactorily. Now the input parameters with their boundary condition 
are taken into consideration from Table 9.1. 

The computation of the algorithm is performed with the following considera- 
tions. The population type is a double vector having a population size of 50. The 
functions are selected in two tournaments. During the analysis, 0.8 is the crossover 
fraction with 1 as the intermediate crossover function. The forward migration with 
0.2 fraction and interval of 20 has been considered. Here, 0.35 fraction is set for the 
Pareto front population. The stopping criterion by considering 100 stall generations 
is 100 x numbers of variables, where the stall limit is infinity. The tolerance of the 
function is 0.0001 with tolerance for the constraint is 0.001. 

Therefore, the total number of iterations completed is 158. There are 17 number 
of parametric combinations along with their objective functions. The final result is 
shown in Table 9.2. When the average change in the dispersion of the Pareto solu- 
tion is achieved toward its limiting condition, the program is terminated. The con- 
tradictory outputs are optimized simultaneously by this GA, which is represented 
by various graphs in Figure 9.11. From this figure it is imperative that the control 
parameters fluctuate within a range. Here MRR value ranges between 0.085 and 
0.163 g/min, respectively, while EWR varies between 0.006 and 0.018 g/min, re- 
spectively. Therefore, for attaining at an optimal parametric combination where 
this dual objective gets satisfied, fuzzy gray relational analysis (FGRA) is done. 
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Table 9.2: Variation of parameters with their responses. 


S. no. Current Voltage Pulse on time MRR EWR 
1 6 45 50 0.085 0.006 
2 10 45 52 0.163 0.018 
3 9 45 51 0.136 0.014 
4 8 45 50 0.116 0.011 
5 8 45 51 0.126 0.012 
6 7 45 50 0.096 0.007 
7 9 45 50 0.152 0.017 
8 6 45 51 0.092 0.007 
9 9 45 51 0.138 0.014 
10 9 45 51 0.151 0.016 
11 8 45 52 0.131 0.013 
12 8 45 50 0.119 0.011 
13 10 45 52 0.158 0.017 
14 7 45 51 0.099 0.008 
15 10 45 51 0.159 0.018 
16 9 45 52 0.142 0.015 


17 6 45 51 0.087 0.006 


9.6 Multiobjective gray relation analysis 


Gray relation analysis (GRA) is a multiobjective optimization technique developed 
by Deng in 1989. Ever since its inception, it has found its presence in experimental 
design where the model can be unspecified or the information can be insufficient, 
thereby reflecting a mixture of familiar and unfamiliar. Thus, GRA bridges the cor- 
relation among the expected result (finest/perfect) with the actual experimental 
findings. Gray grade is the average of the gray coefficient. The calculated gray 
grade denotes that the relation is robust if the value is toward 1. If the value is to- 
ward 0, the stability of the model is very poor. 

In order to obtain the maximized performance value for some parameter like 
MRR, then the equation is normalized as follows: 
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Yong - Min[Ypq,p = 1, 2,... n] 


12 
Max|Yy p = 1, 2,....n] - Min[Ypg,p = 1,2, ..., n] Si 


Xpq = 


In order to obtain the minimized performance value for some parameter like EWR, 
then the equation is normalized as follows: 


Max|[Yj,, i 2 1,2, ..., n] — Ypq 
Max|Yy, p = 1,2, ... n] - Min[Ypg,p = 1,2, .... n] 


Ka (9.13) 


The correlation among definite normalized data along with reference can be ex- 
pressed by the gray relation coefficient as follows: 


Vmin + GV max 
Vg + GV max 


Y (Xoq» Xp) = [D- 1,2, .., n and q-1,2,..., m] (9.14) 


where ` Va- a — Kal, Vmin = Min[Vp p - 12, ..,n and q=1,2,..., m] and Vmax = Max 
[Vog P 21,2, .., n and q- 1,2, .., m]. 

¢ denotes the prominent coefficient (vary within 0-1). 

In general, a suitable fraction can agree with the prominent coefficient as 
follows: 


1 m 
T (Xo, Xp E» Y( (Xog X pa) (9.15) 
where m denotes the quantity of process control parameter. 

A single arrangement of process parameter can be achieved by the help of 
GRA. Here the parametric combination having rank 1 satisfies the multiresponse 
optimization. 


9.7 Concept of fuzzy set 


Whenever the problem of uncertainty in decision making arises, fuzzy set concept 
is applied. Instead of the numerical values, the linguistic calculation of weights is 
done here. The transformation of the decision matrix to the fuzzy decision matrix 
can be influenced by the decision makers' fuzzy ratings (Keufmann & Gupta, 1985). 
Finally, it creates the normalized fuzzy decision matrix. 

An element X, which is a part of a fuzzy set p can be expressed by its member- 
ship function p(x), where “x” denotes the number of the membership and d maps 
all elements in X to the actual number in the limit of 0-1. At this point, triangular 
fuzzy number (TFN), p, can be expressed as a triplet (pi, p2, p3) and the membership 
function is defined [12] as follows: 
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0, X€p1 
X-p1 
<p> «6 EE 
ux) - is (9.16) 
Eeer P2SX<p3 
0, X> p3 


The defuzzification is generation of a crisp value that relates to the formation of a 
nonfuzzy number from a fuzzy number 

Here best nonfuzzy performance (BNP) number is calculated by using “centroid 
of area” technique: 


BNP = Ip: — pl S (ps - pi)] +pı (9.17) 


Selecting the linguistic terms for criteria, 
linguistics terms are utilized in order to give the weight of various criteria. By 
applying TFNs, the linguistic terms can be easily expressed as given in Table 9.3. 


Table 9.3: Linguistic variables in case of criteria. 


Linguistic variables Fuzzy number 


Very high (VH) (0.9,1.0,1.0) 
High (H) (0.7,0.9,1.0) 
Moderately high (MH) (0.5,0.7,0.9) 
Moderate (M) (0.3,0.5,0.7) 
Moderately low (ML) (0.1,0.3,0.5) 
Low (L) (0.0,0.1,0.3) 
Very low (VL) (0.0,0.0, 0.1) 


9.7.1 Fuzzy gray relational analysis 


In this chapter, the two objective functions considered are MRR and EWR. They re- 
flect different consequences while machining. During experimentation, more im- 
portance is imparted on EWR as compared to MRR [13]. While applying fuzzy logic, 
the weights of different criteria are shown in Table 9.4. 


Table 9.4: Weights selected for the outputs. 


Responses Linguistic variables Fuzzy number BNP 


MRR MH 0.5,0.7,0.9 0.436 
EWR H 0.7,0.9,1 0.564 
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The highest gray relation grade value reflects that the result has reached the opti- 
mal condition. Therefore, the optimized machining condition is reached when the 
gray relation grade value is maximum. From the output of GA 17 number of results 
are obtained, which are basically the various experimental combinations. Now 
while applying MCDM, the conflicting nature of the objectives, that is, high MRR 
and low EWR, gets simultaneously satisfied. 


9.7.2 Process control parameter optimization 
Considering the respective weightage for MRR to be 43.6% and for EWR as 56.4%, 


respectively, in the fuzzy set, the gray relation coefficient is determined. The corre- 
sponding gray grade calculated is illustrated in Table 9.5. 


Table 9.5: Gray relation coefficient and grades with corresponding ranks. 


Experiment sequence Objectives Gray grade Position 
MRR EWR 
1 0.085 0.006 0.65 2 
2 0.163 0.018 0.68 1 
3 0.136 0.014 0.51 13 
4 0.116 0.011 0.51 14 
5 0.126 0.012 0.50 17 
6 0.096 0.007 0.57 8 
7 0.152 0.017 0.58 7 
8 0.092 0.007 0.59 6 
9 0.138 0.014 0.51 12 
10 0.151 0.016 0.57 9 
11 0.131 0.013 0.50 16 
12 0.119 0.011 0.50 15 
13 0.158 0.017 0.63 5 
14 0.099 0.008 0.56 10 
15 0.159 0.018 0.64 3 
16 0.142 0.015 0.53 11 


17 0.087 0.006 0.63 4 
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It is imperative from the table that the second experiment gives the maximum gray 
grade. So, it can be inferred that this experimental sequence satisfies multicriteria 
decision-making condition. Here the parametric combination of current (10 amp), 
voltage (45 volts), and pulse on time (52 us) is best suited for obtaining high MRR 
and low EWR. Finally, validation experiment is conducted, while applying the 
above-mentioned experimental combination through gray relational grade the MRR 
is 0.163 and 0.0180 g/min, respectively. The comparison between predicted results 
and experimental results is shown in Table 9.6. 


Table 9.6: Results of predicted performance experimental findings. 


Criteria Estimated result Analytical result 
MRR 0.163 0.162 
EWR 0.0180 0.0182 
Gray relation grade 0.68 0.60 


Enhancement of the gray relation grade by 0.08 


9.8 Conclusion 


Nowadays, several nature-inspired metaheuristics help in solving complex optimiza- 
tion problems effectively. The multiobjective yet contradictory natured objective func- 
tions are easily optimized with these algorithms. The fundamental concepts of these 
algorithms are taken from nature. ANN can be effectively used to train as well as vali- 
date the experimental results. GA provides global optimal solutions to the multiobjec- 
tive optimization problems. Finally, from the results obtained from the GA, a single 
parametric combination is derived using fuzzy multicriteria decision-making tech- 
nique where contradictory objectives are satisfied simultaneously. Therefore, it is evi- 
dent that EDM process can effectively machine HCHCr steel while varying current, 
voltage, and pulse on time, respectively, by using titanium nitride-coated copper 
tool. Through ANN, it is evident that the objective to maximize MRR and minimize 
EWR has been successfully accomplished since the value for R in both the cases is 1. 
From the GA algorithm, the values of MRR range between 0.085 and 0.163 g/min, and 
the values of EWR vary between 0.006 and 0.018 g/min, respectively. Considering the 
range of values of MRR and EWR, finally FGRA is applied. Lastly from the analysis, the 
parametric combination having current (10 amp), voltage (45 volts), and pulse on time 
(52 us) provides the maximum MRR and minimum EWR simultaneously. 

Now EDM is a very complex process having numerous control parameters. Today it 
has emerged as a key machining system in the manufacturing domain. Commercially, it 
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is also proving to be very economical and productive. The ranges of these parameters 
are quite large. But through this experimental study, the exact value of the control pa- 
rameters through which the effective machining can be done has been established. This 
will be useful to the operators who are working on it commercially. 
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10 Fuzzy reliability of system using different 
types of bifuzzy failure rates of 
components 


Abstract: The present study proposes to determine the fuzzy reliability of systems 
using different types of conflicting bifuzzy numbers Till now, to evaluate the fuzzy 
reliability used the failure rates of components of systems follow the same type of 
fuzzy set or conflicting bifuzzy set. However, in real-life problems, such type of situ- 
ation rarely occurs. Therefore, in this study, a new method has been introduced to 
determine the fuzzy reliability of a system having components following different 
types of conflicting bifuzzy failure rates. Using the introduced method, membership 
functions and nonmembership functions of fuzzy reliability of a series, parallel, 
parallel-series, and series—parallel systems are evaluated. Numerical problems are 
also taken to describe the proposed study. 


Keywords: fuzzy system reliability, conflicting bifuzzy number, failure rate, system 


10.1 Introduction 


The traditional reliability of an item is the probability that item will perform a re- 
quired function without failure under stated conditions for a stated period of time. 
More generally, reliability is the capacity of parts, components, equipments, prod- 
ucts, and systems to perform their required functions for desired periods of time 
without failure, in specified environments and with a desired confidence. However, 
in real-life problems, the collected data or information are not always accurate 
hence the deduction of appropriate value of probability becomes very difficult in 
many cases. To handle this problem, fuzzy set theory has been used to analyze the 
reliability of a system. The fuzzy set theory introduced by Zadeh [1]. 

To understand the concept of fuzzy set, let us consider a system having three 
independent components. Evidently, the system is in working stage when all three 
components are in working condition simultaneously. However, when one or two 
components are failed, the system will operate in an offended state. In this stage, 
the system is neither fully working nor fully failed, but is in some intermediate 
state. This implies that the binary state assumption for delineate system failure and 
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success may no longer be suitable. To handle this problem, we fuzzify the system 
failure. In real-life problems a person may assume that a commodity belongs to a 
set to a few degrees, but he is not sure that a commodity belongs to a set. In other 
words, there may be a hesitance or uncertainty about the membership function of 
object in A. An eventually solution is to use intuitionistic fuzzy sets (IFS), defined 
by Atanassov (1986). 

The concept of IFS is the generalizations of the fuzzy set theory of Zadeh [1]. 
The membership function is considered only in fuzzy set but membership degree 
a(x) and nonmembership function v4(x) are considered in IFS and the sum of 
membership and nonmembership degree is less than or equal to 1. Suppose if the 
sum of the value of membership and nonmembership function is greater than 1 
then how do we fulfill the condition, for example, in a case where the membership 
function H, (x) = 0.8 and the nonmembership function v4(x) equals to 0.6? This situ- 
ation gives us a chance to redefine a solution lastly to reform the concept of IFS. 
Based on this reason Zamali et al. [2] gave a new concept of conflicting bifuzzy set 
(CBFS). In bifuzzy set, the sum of membership and nonmembership degree is less 
than 2. Hence, CBFS theory is an extension of IFS theory. Imran et al. [3] presented 
the definition and graphical representation of conflicting bifuzzy sets based on in- 
tuitionistic evaluation. 

Fuzzy reliability theory was introduced and developed by Wen et al. [4, 5], 
Chen and Mon [6], Onisawa and Kacprzyk [7], and Cai [8]. The recent compilation of 
papers by Onisawa and Kacprzyk [7] gave several different approaches for fuzzy 
reliability. Kumar and Yadav [9] proposed a new method for analyzing fuzzy reli- 
ability of systems using different types of intuitionistic fuzzy failure rates of compo- 
nents. Aliev et al. [10] analyzed fuzzy reliability of systems using time dependent 
fuzzy set of components. Mon and Cheng [11] developed a new method for fuzzy reli- 
ability of systems for components with different membership functions. Kumar et al. 
[12] developed a new method for evaluating fuzzy reliability using intuitionistic fuzzy 
number. Mon et al. [11] analyzed the fuzzy reliability of different systems with differ- 
ent membership functions. 

In the above-mentioned studies, it is supposed that the failure rates of components 
of systems follow the conflicting bifuzzy numbers (CBFNs) to evaluate the fuzzy reli- 
ability of systems. In this study, we have constructed the membership degree and non- 
membership degree of the fuzzy system reliability using different types of CBFNs. 
Using introduced method, membership degree and non-membership degree of fuzzy 
system reliability of series, parallel, parallel-series, and series-parallel systems are 
evaluated. Numerical problems are also taken to exemplify the proposed study. 
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10.2 Preliminaries 


Definition 10.1 ([1]). If a set Y be fixed, then a fuzzy set Z is given by 
Z- (o. uz): € Y] 


where uz(x) € [0, 1] is the membership function of the element y € Y. 


Definition 10.2 (Atanassov 1986). If a set Z be fixed then an intuitionistic fuzzy set 
(IFS) A is given by 


A= {(z, uy (2), vq(z):z e Z)} 


where j;:Z — [0, 1] is the membership function of positive z (zc Z) and 
v4:Z — [0, 1] is the non- membership function of negative z, (z € Z) w.r.t. A and 
Oxgu;-tvixl. 

Zamali [2] in his study proposed the new notion called a bifuzzy set which is 
extension of IFS. He defines CBFS. 


Definition 10.3 ([2]). If a set A is fixed then a CBFS Z of A is an object having the 
given form as follows 


Z={(a,uz(a),vz(a):a € A)} 
where uz: A — [0, 1] is the membership function of positive a, (a € A) with respect 


to Z and v;:A — [0, 1] is the nonmembership function of negative a, (a € A) with 
respect to Z and 0 < u;(a) + vz(a) <2. 


10.2.1 a-Level set or a-Cut set of a CBFS 


a-Cut of a CBFS À is defined as 


Aq = (x e xiu) 2a, };0<a 


10.2.2 B-Level set or B-Cut set of a CBFS 


B-Cut of a CBFS A is defined as 


Ap = {x e X:vq(x) «Bh B<1 
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10.2.3 Triangular CBFS 


A triangular CBFS A is defined as 
A=(m-Lm,m+n;m-l',m,m+n’) 
where m c R is the center, |» O and n» 0 are the left and right dispersion of the 


membership function of A, respectively. |» 0 and n' » 0 are the left and right 
spreads of the nonmembership function of A, respectively. 


10.2.4 Trapezoidal CBFS 
A= (p-La.ns+np-l,qnstn) 


where p,q,r,s € R is the center, 1» 0 and n» 0 are the left and right dispersion of 
the 11; (x), respectively. IT >0 and n' » 0 are the left and right dispersion of the va(x), 
respectively. 


10.3 Bifuzzy reliability function or survival function 


The reliability function R (t) of systems in terms of failure rate A (t) is given by fol- 
lowing formula 


t 


R(t) = exp l- poa, t>0 (10.1) 
0 


Let us consider the failure rate is constant, that is, (t) =A, then the reliability is as 
follows: 


R(t) = exp( - At) (10.2) 


If failure rate is not precise, then failure rate is considered as CBFS. 
Let failure rate is considered as CBFS A defined on Z given below: 


A={ < Z, py (Z),Vi(Z) >:z € Z} 


Two types of cuts for CBFS are defined as follows: 
a-Cut of membership function p;(z) which is defined as 


Ña ={z: u; (z) za, € [0, 1]) 
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B-Cut of nonmembership function v;(z) is as 
Ag = {2:V;(Z) < B,B € (0, 1]} 
Failure rate A is in form of CBFS. Then fuzzy-convexity of CBFS is as follows: 
À, = [AL AS] va e [0,1] 
Also fuzzy-concavity of CBFS is as follows: 
Àg = Ha Ag] VB € [0.1] 


where A and A increasing and decreasing function of a and f, respectively. Aj and 
A increasing and decreasing function of a and £, respectively, with a, B € [0,1]. 
The reliability function R(t) using CBFS defined on X is as follows: 


R(t) = { «x Uro (X) Vro (X) » :x € X) 


Now failure rate follows CBFS, then the reliability R(t) is also in the form of CBFS. 
Now a-Cut for membership function of R(t) is 


where 
r(t)= min(e-“) st A <x<A? (10.3a) 
rR(t)= max(e^*) s.t Ab <x<ar (10.3b) 


B-Cut for nonmembership function of R(t) is given by 


Rat = air, r$ (0)] 
where 
r(t)= min(e™) s.t. A <x <À (10.4a) 
rg(t)= max(e *) s.t. Ap <x < Ap (10.4b) 


For membership and nonmembership function, a-Cut and -Cut of R(t) respectively, 
is 


Ralt) = [e73 eat] (10.5a) 


Rg(t)-|e 8 ,e £] (10.5b) 
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10.4 Fuzzy reliability of series system using 
different types of conflicting bifuzzy 
failure rates 


Consider a system having n components connected in series (Figure 10.1). Let fail- 
ure rates of n components are represented by different CBFN Aj, A, ..., Àn and reliabil- 
ity with conflicting bifuzzy failure rates An, À>, ...,A, are Rj, Ro,...,Rn, respectively, at 
time t. 


€ 


Figure 10.1: Series system. 


Reliability of series system at time t is 


Rs(t) = Ri(t)@R2(t)® - -- @Ra(t) (10.6) 


From eq. (10.2), we get 


R(t) =e", i=1,2,3, ..., n t»0 (10.7) 


The reliability of series system R,(t) is estimated as 


(10.8) 


where Ag = (Aj, Ad, ... An) = XL J is failure rate. 


a-Cut of Às for membership function is 


Àj, = [hl 


jæ 


hE] Yace (0,1), j212, ....n 


We have n intervals Aiq= [hL,, hE], Jo, = [Wb hE] -o Ae (he, MÉ VÀ, 1<j<n, 
respectively, for every a. 
p-Cut of nonmembership function for As is 


Avg = (hig, hh] VB € [0,1], j=1,2, ....n 


Similarly, we have n intervals Ae = [hig Mi], Aog = [hip h3], -- -> Ang = [hip hig] VA), 
1<j<n, respectively, for every £. 
a-Cut of As for membership function is 


Asa = [US UE] Va € [0.1] 
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where 
E ei n RL n 
Wo, = min ` jai Ysa = max ji 
L R pl R L R 
s.t. hig < X1 < hio ha < X2 < eus ha < Xn < Ra (10.9) 


B-Cut of A, for nonmembership function is 
1 L pR 
Asp = [Psg Psp) VB € [0.1] 
where 
. n n 
d om y» d = may 242 


s.t. Nig < xı < hip, hig ox < Nig, --» hig < X1 < hig (10.10) 


After solving (10.9) and (10.10), we find a-Cut and f-Cut of fuzzy system reliability. 
a-Cut and f-Cut of membership function of R,(t) as 


Rsa(t) = [e^ ele fit (10.11a) 
Rep(t) = je ri. e "Tat (10.11b) 


o8- | 
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d --- Nonmembership 
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Figure 10.2: For f= 10, membership and nonmembership function of reliability of radio set. 
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Figure 10.4: For t= 30, membership and nonmembership function of reliability of radio set. 


10.4.1 Numerical problem 


Consider a system consisting of three independent components and all components 
work properly to operate the system successfully. Let failure rates of components 
are in the form of different types of CBFS Ay, Az, and A3, respectively, assumed val- 
ues of Ài, Ao, AG are given in Table 10.1. 


Table 10.1: Conflicting bifuzzy failure rates of components. 


Bifuzzy failure rate of ith component 
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A; = (0.02, 0.05, 0.08; 0.01, 0.05, 0.09) 


Types of bifuzzy number 


Triangular conflicting bifuzzy number 


Az = (0.03, 0.06, 0.08; 0.02, 0.07, 0.09) 


Triangular conflicting bifuzzy number 


A3 = (0.02, 0.05, 0.08, 0.12; 0.01, 0.07, 0.10, 0.13) 


Trapezoidal conflicting bifuzzy number 


a-Cuts of A, Ao, Às are estimated and tabulated in Table 10.2. 


Table 10.2: a-Cut of À;, A> and A for membership. 


a Ata Aza Asa 

0.0 [0.02,0.08] [0.03, 0.08] [0.02, 0.12] 
0.1 [0.023,0.077] [0.033, 0.078] [0.023, 0.116] 
0.2 [0.026,0.074] [0.036, 0.076] [0.026, 0.112] 
0.3 [0.029,0.071] [0.039, 0.074] [0.029, 0.108] 
0.4 [0.032,0.068] [0.042, 0.072] [0.032, 0.104] 
0.5 [0.035,0.065] [0.045, 0.070] [0.035, 0.100] 
0.6 [0.038,0.062] [0.048, 0.068] [0.038, 0.096] 
0.7 [0.041,0.059] [0.051, 0.066] [0.041, 0.092] 
0.8 [0.044,0.056] [0.054, 0.064] [0.044, 0.088] 
0.9 [0.047,0.053] [0.057, 0.062] [0.047, 0.084] 
1.0 [0.05, 0.05] [0.06, 0.06] [0.05, 0.08] 


p-Cuts of A, A; A; are estimated and tabulated in Table 10.3. 


Table 10.3: B-Cut of A, Az, and A3 for nonmembership. 


B Asp Aap Asp 

0.0 [0.05, 0.05] [0.07, 0.07] [0.07, 0.10] 
0.1 [0.046, 0.054] [0.065, 0.072] [0.064, 0.103] 
0.2 [0.042, 0.058] [0.060, 0.074] [0.058, 0.106] 
0.3 [0.038, 0.062] [0.055, 0.076] [0.052, 0.109] 
0.4 [0.034, 0.066] [0.050, 0.078] [0.046, 0.112] 
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Table 10.3 (continued) 


p Ay Ang Asp 

0.5 [0.03, 0.07] [0.045, 0.080] [0.040, 0.115] 
0.6 [0.026, 0.074] [0.040, 0.082] [0.034, 0.118] 
0.7 [0.022, 0.078] [0.035, 0.084] [0.028, 0.121] 
0.8 [0.018, 0.082] [0.030, 0.086] [0.022, 0.124] 
0.9 [0.014, 0.086] [0.025, 0.088] [0.016, 0.127] 
1.0 [0.01, 0.09] [0.02, 0.09] [0.01, 0.13] 


Using formulas (10.9) to (10.11), we get fuzzy reliability of series system 
Rsa(t) = le^ (.28 - .09a)t e — (.07 + .09a)t ] (10.122) 
Rsg(t) 2 C e (09-55 (10.12b) 
The fuzzy system reliability of given system is calculated by eq. (10.12). a-Cut and 
p-Cut of fuzzy reliability are computed for t = 10, 20, 30 shown in Table 10.4 and fig. 


10.2, 10.3 & 10.4. The fuzzy reliability of above system is in trapezoidal CBFN for 
failure rates as conflicting bifuzzy of components for different time t. 


10.5 Fuzzy reliability of parallel system using 
of conflicting bifuzzy failure rates 


Consider a system consisting “n” components are in parallel form (Figure 10.5). Let us 
consider failure rates of n components are in form of different CBFNs Ab Ao, er Àn. 
Let the fuzzy reliability of n components with conflicting bifuzzy failure rates are 
Ab Àa, " d and RB, Ro, ..., Rn respectively, at time t. 

Reliability of parallel system at time t is 


Rp(t)=1- |[a-R(0)-1- [[ 0- ep - 4.5) (10.13) 
e - ~ m n p 
So Rut = Va, Ao, .. Aut) -1- JI (1 - exp( - A0) 
a-Cut of Ai for membership function is 
Nia = (hi, hR] va € [0,1], i=1,2,3,.. n 


We have n intervals Ae - IN, hÈ], Ao, = (RE, We], ..., Ana = (hi, hE] vA; leien 
respectively, for each a. 
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Figure 10.5: Parallel system. 


p-Cut of nonmembership function of A is 
Aig = (hi, hip] VB € [0,1], i=1,2,3,....n. 
Similarly, we have n intervals Ae = (hig, Ai], Aog = [hig Is], - -> Ang = [htp Kal VÀ, 


1<j<n, respectively, for each f. 
a-Cut of Rp(t) for membership function is 


Rpa = AOF Wra(t)] 


where 


= 


j=1 


Wralt) = min H II ae), Dealt) = max(1- a ven (10.14) 


L R pL R L R 
s.t. hig < X1 < hiy hu < X2 < ho e Ana < Xn < Ppa- 


B-Cut of nonmembership function of Rp(t) is as 
Reg = Wiel, Ppp (t)] 


where 


n 


del) = min(1- I] ae), dal = max(1- i nei 


i=1 (10.15) 


s.t. hig < xı < his, hig < X1 < Na, ig < X1 < Pi. 


After solving eqs. (10.14) and (10.15), we find a-Cut and f-Cut of fuzzy reliability of 
parallel system Rp(t). 
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10.5.1 Numerical problem 


Consider a system has three independent components. At least one of them must 
work properly to operate given system properly. Let conflicting bifuzzy failure rates 
of components are in the form of CBFNs A;, À>, and A3, respectively. For numerical 
computation assumed values of A, Ao, and A3 are given in Table 10.5. 


Table 10.5: Conflicting bifuzzy failure rates of components. 


Bifuzzy failure rate Types of bifuzzy number 

A = (0.3, 0.4, 0.5; 0.2, 0.4, 0.6) Triangular conflicting bifuzzy number 
Az = (0.2, 0.5, 0.7; 0.3, 0.6, 0.8) Triangular conflicting bifuzzy number 
A3 = (0.03, 0.05, 0.08, 0.10; 0.02, 0.06, 0.10, 0.12) Trapezoidal conflicting bifuzzy number 


a-Cuts of Aj, Ay, A; are evaluated and given in Table 10.6. 


Table 10.6: For membership, a-Cut of Ai, Az, and A3. 


a Aa Aaa Asa 

0.0 [0.3, 0.5] [0.2, 0.7] [0.03, 0.10] 
0.1 [0.31, 0.49] [0.23, 0.68] [0.032, 0.098] 
0.2 [0.32, 0.48] [0.26, 0.66] [0.034, 0.096] 
0.3 [0.33, 0.47] [0.29, 0.64] [0.036, 0.094] 
0.4 [0.34, 0.46] [0.32, 0.62] [0.038, 0.092] 
0.5 [0.35, 0.45] [0.35, 0.60] [0.04, 0.09] 
0.6 [0.36, 0.44] [0.38, 0.58] [0.042, 0.088] 
0.7 [0.37, 0.43] [0.41, 0.56] [0.044, 0.086] 
0.8 [0.38, 0.42] [0.44, 0.54] [0.046, 0.084] 
0.9 [0.39, 0.41] [0.47, 0.52] [0.048, 0.082] 


1.0 [0.4, 0.4] [0.5, 0.5] [0.05, 0.08] 


198 — Pawan Kumar, S.B. Singh 


p-Cut of Ais Ao, and is are evaluated and given in Table 10.7. 


Table 10.7: For nonmembership, B-Cut of Ai, Ar, and A3. 


B Asp Aap Ase 

0.0 [0.4, 0.4] [0.6, 0.6] [0.06, 0.10] 
0.1 [0.38, 0.42] [0.57, 0.62] [0.056, 0.102] 
0.2 [0.36, 0.44] [0.54, 0.64] [0.052, 0.104] 
0.3 [0.34, 0.46] [0.51, 0.66] [0.048, 0.106] 
0.4 [0.32, 0.48] [0.48, 0.68] [0.044, 0.108] 
0.5 [0.3, 0.5] [0.45, 0.70] [0.04, 0.11] 
0.6 [0.28, 0.52] [0.42, 0.72] [0.036, 0.112] 
0.7 [0.26, 0.54] [0.39, 0.74] [0.032, 0.114] 
0.8 [0.24, 0.56] [0.36, 0.76] [0.028, 0.116] 
0.9 [0.22, 0.58] [0.33, 0.78] [0.024, 0.118] 
1.0 [0.2, 0.6] [0.3, 0.8] [0.02, 0.12] 


Using formulas (10.14) to (10.15), we get 


p, (t -1- (1- e (5 lat) x (1 e-( (7- 200) x (1 — e: 922)t 
R(t) =1- (17 e7 (341M!) x (1— e- (239) x (4 — g- (03.020) 
Lon nmn  p-(A* 2t (64 2B)t - (1«.02)t (10.16) 
Wrp(t) -1- (1- e )x(1-e7 )x(1-e ) 
We g(t) =1- (17 ec 729) x (17 e- (673816) x (1 — g- (06-04 


The system reliability of given system is evaluated using (10.16). a-Cut and fi-Cut of 
fuzzy reliability of given parallel system are computed for t = 10, 20, 30 and are tab- 
ulated in Table 10.8 and Fig. 10.6, 10.7 & 10.8. The fuzzy reliability of system is in 
trapezoidal CBFN. 


10.6 Reliability of parallel-series system when 
failure rates follow conflicting bifuzzy numbers 


Let us consider a parallel-series system consisting “j” branches connected in paral- 
y: 
lel form and each branch contains “” components connected in series form as 
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Figure 10.6: For t= 10, membership function and nonmembership function of reliability of 
achieving orbit. 
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Figure 10.7: For t= 20, membership function and nonmembership function of reliability of 
achieving orbit. 


shown in Figure 10.9. Let us consider failure rates of all components are following 
different CBFNs Ars 

Let failure rate function of the sth component of rth branch (r=1, 2,..., iand 
s=1,2,...,j) is represented by CBFS Aves 

Reliability of the sth component of rth branch is as follows: 


R,s(t) = exp ( - s.t) 


It is well known that the fuzzy reliability of system Rps(t) is given by 
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7 — Membership 
^ --- Nonmembership 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


Figure 10.8: For t = 30, membership function and nonmembership function of reliability of 
achieving orbit. 


Figure 10.9: Parallel-series system. 


Rps(t) =1- Il (1- Un, =1- Il D: Ten An (10.17) 


s-1 r=1 s-1 r=1 


~ H i ~ 
So Brett) = (Ars: = 1, 2.. i and s=1,2,...,5 =1- [T D: Ten Aan) 


a-Cut of Ars for membership function is sah r=l 
Avsa = Dën hE] Ya E€ [0,1],r=1,2,..,iand s=1,2....j 


For each a, we have n? intervals Ara = [hl hE,], Aiza = [Rog hal Anta = [AL RE] 


EA Aja = [hi ba, for each À,, 1<r<iand 1<s<j, respectively. 


ijj? Es 
For nonmembership function, f-Cut of A, is 


Aug Wie Miel VB € [0,1],7=1,2,...,iand s=1,2...j 


Similarly, for every f, we have m^ intervals Aug [hhg hs], Arp = [hhg Mäe), 
Ang = ës, bc, os Ag = [hips hi]; for each A, 1<r<iand 1« s x j, respectively. 
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a-Cut of Rps(t) for membership function is 


Rpsa(t) = Wess deel D) 


S 


r 
R 
uz ff hi, s SXp&z Be hias EZE We se hij a =X NM. 


B-Cut of Rps(t) for nonmembership function is 


Rpsp(t) = [Ups (t), Wisp (£)] 


Wosp(t) = min 
j 


h-IIG- [ew Ac) Voss (f) )- min(1- li Hen -X)] 


(10.18) 


" " (i JJe- int) \, dal? min{ 1- II (-I Hen EU 9)! 


s=1 s=1 
L R L R L R L R 
s.t. ue < X11 < hip hog € < hop hag < Aa < hap — Nig < Xij < his 


(10.19) 


After solving (10.18) and (10.19), we calculate a and f-Cut of fuzzy reliability of 


above system Rps(t). 


10.6.1 Numerical example 


Suppose there is a communication system which receives the input signal and trans- 
mits the output signal. For this, there are two receivers and two transmitters in the 
system connected as is in Figure 10.10. For a successful communication at least one 
receiver and one transmitter connected in series configuration must work properly. 
The failure rates of receivers and transmitters are in the form of different types 


of CBFNs given in Table 10.9. 


The a-Cuts of Aii Ans An, A» are computed and tabulated in Table 10.10. 
The f-Cut of Ai An, An, Ao are computed and tabulated in Table 10.11. 


Using formulation (10.18) and (10.19), we get 
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Ry (Transmitter 1) 


(Receiver 1) 


Input Output 


R- (Receiver 2) 


R5, (Transmitter 2) 


Figure 10.10: Parallel-series system. 


Table 10.9: Conflicting bifuzzy failure rates of components. 


Bifuzzy failure rate 


Ay. = (0.05, 0.10, 0.15; 0.04, 0.12, 0.18) 


Types of bifuzzy number 


Triangular conflicting bifuzzy number 


Ai = (0.06, 0.09, 0.12, 0.16; 0.07, 0.11, 0.15, 0.18) 


Trapezoidal conflicting bifuzzy number 


izn = (0.07, 0.12, 0.17; 0.09, 0.14, 0.18) 


Triangular conflicting bifuzzy number 


Az = (0.03, 0.08, 0.13, 018; 0.05, 0.10, 0.18, 0.23) 


Table 10.10: a-Cut of Ai, Aus, A21, and Ass for membership. 


Trapezoidal conflicting bifuzzy number 


a Ana Anza Ae, Aas, 

0.0 [0.05, 0.15] [0.06, 0.16] [0.07, 0.17] [0.03, 0.18] 
0.1 [0.055, 0.145] [0.063, 0.156] [0.075, 0.165] [0.035, 0.175] 
0.2 [0.06, 0.14] [0.066, 0.152] [0.08, 0.16] (0.04, 0.17] 
0.3 (0.065, 0.135] [0.069, 0.148] [0.085, 0.155] [0.045, 0.165] 
0.4 [0.07, 0.13] [0.072, 0.144] [0.09, 0.15] [0.05, 0.16] 
0.5 [0.075, 0.125] [0.075, 0.140] (0.095, 0.145] [0.055, 0.155] 
0.6 [0.08, 0.12] [0.078, 0.136] [0.10, 0.14] [0.06, 0.15] 
0.7 [0.085, 0.115] [0.081, 0.132] [0.105, 0.135] [0.065, 0.145] 
0.8 [0.09, 0.11] [0.084, 0.128] [0.11, 0.13] [0.07, 0.14] 
0.9 [0.095, 0.105] [0.087, 0.124] [0.115, 0.125] [0.075, 0.135] 
1.0 [0.10, 0.10] [0.09, 0.12] [0.12, 0.12] [0.08, 0.13] 
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Table 10.11: B-Cut of Ans Ay, dar and Lë for nonmembership. 


B Aug Arp Àng Azza 
0.0 (0.12, 0.12] [0.11, 0.15] [0.14, 0.14] [0.10, 0.18] 
0.1 [0.112, 0.126] [0.106, 0.153] [0.135, 0.144] [0.095, 0.185] 
0.2 [0.104, 0.132] [0.102, 0.156] [0.130, 0.148] [0.09, 0.19] 
0.3 [0.096, 0.138] [0.098, 0.159] [0.125, 0.152] [0.085, 0.195] 
0.4 [0.088, 0.144] [0.094, 0.162] [0.120, 0.156] [0.08, 0.20] 
0.5 [0.08, 0.15] [0.090, 0.165] [0.115, 0.160] [0.075, 0.205] 
0.6 [0.072, 0.156] [0.086, 0.168] [0.110, 0.164] [0.07, 0.21] 
0.7 [0.064, 0.162] [0.082, 0.171] [0.105, 0.168] [0.065, 0.215] 
0.8 [0.056, 0.168] [0.078, 0.174] [0.10, 0.172] [0.06, 0.22] 
0.9 [0.048, 0.174] [0.074, 0.177] [0.095, 0.176] [0.055, 0.225] 
1.0 [0.04, 0.18] [0.07, 0.18] [0.09, 0.18] [0.05, 0.23] 
Ger t)=1-{ 1-e- (36- a)t) x(1- e (34-.08a)tyy 
(10.20) 


( ( 

R (t) -1-((1- ern y (1 — g- (09 +.08a)ty 
( (1— e7 (30386) x (1— e~ (33+.088)t) 1 
( ( 


t) =1-{ 1- e7 (26-89) x (ise ty) 


The fuzzy system reliability of given system is evaluated using (10.20). The fuzzy 
reliability of system is trapezoidal CBFN for failure rate of components follow con- 
flicting bifuzzy. The a-Cut and f-Cut of fuzzy reliability Rps(t) are evaluated for t = 
10, 20, 30 and shown in Table 10.12. The value of sum of membership and nonmem- 
bership degree is found to be greater than one (Figures 10.11-10.13) at point A, B, C. 


10.7 Reliability of series—parallel system when 
failure rates follow conflicting bifuzzy numbers 


Let us assume a series— parallel system consisting “” subsystems connected in se- 
ries and each subsystem contains “j” components connected in parallel form as 
shown in Figure 10.14. Let failure rates of all components are assumed to be differ- 
ent CBFNs dee 
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Figure 10.11: For t= 10, membership function and nonmembership function of reliability of 
communication system. 
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Figure 10.12: For t = 20, membership function and nonmembership function of reliability of 


communication system. 


Let failure rate function of the rth component of sth subsystem (r=1, 2, .. ., j and 
s-1,2,.. , i) is represented by CBFS A 
Reliability of the rth component of sth branch is as follows: 


R,s(t) = exp( - s.t) 


It is well known that the fuzzy reliability of system Rsp(t) is given below: 


10 Fuzzy reliability of system using different types of bifuzzy —— 207 
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Figure 10.13: For t = 30, membership function and nonmembership function of reliability of 
communication system. 


Figure 10.14: Series—parallel system. 


(Receiver 1) 


snb, 


R35 (Transmitter 1) 


Figure 10.15: Series—parallel system. 
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Table 10.13: Conflicting bifuzzy failure rates of components. 


Bifuzzy failure rate 


- (0.06 , 0.09, 0.12; 0.07, 0.11, 0.16) 


Types of bifuzzy number 


Triangular conflicting bifuzzy number 


Ay = (0.07, 0.12, 0.17, 0.22; 0.09, 0.14, 0.19, 0.24) 


Trapezoidal conflicting bifuzzy number 


Triangular conflicting bifuzzy number 


( 
= (0.06, 0.11, 0.16; 0.08, 0.13, 0.18) 
( 


i22 = (0.08, 0.13, 0.18, 0.22; 0.10, 0.15, 0.20, 0.25) 


Table 10.14: a-Cut of Au, Aus, Aan, and Aa» for membership. 


Trapezoidal conflicting bifuzzy number 


a Anta Aaa Aan Aga 
0.0 [0.06, 0.12] [0.07, 0.22] [0.06, 0.16] [0.08, 0.22] 
0.1 [0.063, 0.117] [0.075, 0.215] [0.065, 0.155] [0.085, 0.216] 
0.2 [0.066, 0.114] [0.08, 0.21] [0.07, 0.15] [0.090, 0.212] 
0.3 [0.069, 0.111] [0.085, 0.205] [0.075, 0.145] [0.095, 0.208] 
0.4 [0.072, 0.108] [0.09, 0.20] [0.08, 0.14] [0.10, 0.204] 
0.5 [0.075, 0.105] [0.095, 0.195] [0.085, 0.135] [0.105, 0.20] 
0.6 [0.078, 0.102] [0.10, 0.19] [0.09, 0.13] [0.11, 0.196] 
0.7 [0.081, 0.099] [0.105, 0.185] [0.095, 0.125] [0.115, 0.192] 
0.8 [0.084, 0.096] [0.11, 0.18] [0.10, 0.12] [0.12, 0.188] 
0.9 [0.087, 0.093] [0.115, 0.175] [0.105, 0.115] [0.125, 0.184] 
1.0 [0.09, 0.09] [0.12, 0.17] [0.11, 0.11] [0.13, 0.18] 
i j i j 
Rs - || (1- Il (1-R,)) - II (1- (UE emm Au (10.21) 

s=1 r=1 s-1 ral j 

So Rsp(t) = Ps: r=1,2,..,jand s=1,2,..,% t= | [{1- [ [a- exp(-A,))} 
a-Cut of A for membership function is zc) Teal 
Asa = [he RR Ya € [0,1], r=1,2,.. j and s-1,2,... i. 


For every a, we have n? intervals Ang = [hi AR a), Ana = [Aha MÉ, Anta = [Aha RR], 
<, Ajia = LA hil, for each As, 1«rxj and 1<s<i, respectively. 


For nonmembership function, f-Cut of Ars is 


Table 10.15: B-Cut of AnsA12,An1, and A for nonmembership. 
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B Aug Än Aas Àaag 
0.0 [0.11, 0.11] [0.14, 0.19] [0.13, 0.13] [0.15, 0.20] 
0.1 [0.106, 0.115] [0.135, 0.195] [0.125, 0.135] [0.145, 0.205] 
0.2 [0.102, 0.120] [0.13, 0.20] [0.12, 0.14] [0.14, 0.21] 
0.3 [0.098, 0.125] [0.125, 0.205] [0.115, 0.145] [0.135, 0.215] 
0.4 [0.094, 0.130] [0.12, 0.21] [0.11, 0.15] [0.13, 0.22] 
0.5 [0.090, 0.135] [0.115, 0.215] [0.105, 0.155] [0.125, 0.225] 
0.6 [0.086, 0.140] [0.11, 0.22] [0.10, 0.16] [0.12, 0.23] 
0.7 [0.082, 0.145] [0.105, 0.225] [0.095, 0.165] [0.115, 0.235] 
0.8 [0.078, 0.150] [0.10, 0.23] [0.09, 0.17] [0.11, 0.24] 
0.9 [0.074, 0.155] [0.095, 0.235] [0.085, 0.175] [0.105, 0.245] 
1.0 [0.07, 0.16] [0.09, 0.24] [0.08, 0.18] [0.10, 0.25] 


Arsp= (hep, Miel vBe[0.1,r-12,..j and s=1,2....,i 


Similarly, for each f, we have n? intervals Jus = [hiig hug) Ňg= [hip hig) Ang = 
[Wop hg]. -- -> Aus = [hig hig], for any Às, 1<r<jand 1<s<i, respectively. 
a-Cut of Rsp(t) for membership function is 


Rspa(t) = [WSpa(t)> Wspa(€)] 
where 
Wspq(t) = min 
i j i j 7 
{TI (1- LL (1- exo(-4n))) ag = min [TT (1- UD em Al 
s=1 r=1 s-1 1 
s.t. ae <Xu< hee hba < SXp< ss hlas EXAS ns iss hi a SX ij < hÈ (10.22) 
B-Cut of Rsp (t) for nonmembership function is 
Rspp(t) = [WSpp(t), VSpg(£)] 


where 


Wspp(C) = min 
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Figure 10.16: For t= 10, membership function and nonmembership function of reliability of 
communication system. 
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Figure 10.17: For t= 20, membership function and nonmembership function of reliability of 
communication system. 


In (1- Il (1- exp( All \ Wspp(t) = Zo D- Il (1- en Aal 


EEN r=1 EEN r=1 
s.t. hig <Xu < bie, be <Xu < bie bie <Xu < bie — be <Xu < hig (10.23) 


After solving (10.21) and (10.22), we get out a-Cut and fi-Cut of fuzzy reliability of 
series—parallel system Rop(t). 
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Figure 10.18: For t = 30, membership function and nonmembership function of reliability of 
communication system. 


10.7.1 Numerical problem 


Suppose there is a communication system which receives the input signal and 
transmits the output signal. For this there are two receivers and two transmitters in 
the system connected as shown in Figure 10.15. For a successful communication at 
least one receiver and one transmitter must work properly. 

The failure rate of receivers and transmitters follow different types of CBFNs as 
shown in Table 10.13. 

The a-Cut of Ay, An, Àa, An are calculated and tabulated in Table 10.14. 

The f-Cuts of Au, Au, An, Az are calculated and tabulated in Table 10.15. 
Using formulae (10.21) and (10.22), we get 

Wp, (t) ={1- (i= — (0.12 - 0.03a)t )x x (1 — e7 (016 -0.05a)t)} x {1- (1-— e7 (022-0.050)t) 


x(1- pn 
VE, (t) = (1— (1 - e- 09600301) x (1 — e~ (09600500) x (1... (1 — e~ (0.07 + 0.050) 

x (1 — e~ (0-08 + 0.0590) 
Wspp(t) = (1- (17 e~ 03140055) x (1 9 (013+ 0.058)t)) x (4 — (1 — e~ (0.19 + 0.055)t) 

x (1 — e~ (020+ 0.058)t) 1 
Wépp(t) = [1- (17 e~ 021-0048) x (1 — g-(0.13-0.058)t) 1 x £4 _ AAR 

x(1- e~ (0.15 -0.058)t)) 

(10.24) 
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The system reliability of above system is calculated using (10.24). The system reli- 
ability is in trapezoidal CBFN for failure rate of components follow conflicting bi- 
fuzzy. The a-Cut and f-Cut of fuzzy reliability of above system Rsp(t) are evaluated 
for t = 10, 20, 30 and given in Table 10.16. The sum of membership and nonmember- 
ship degree is found to be greater than one (Figures 10.16-10.18) at points A, B, C. 


10.8 Conclusion 


In this study, we have discussed a new method for evaluating fuzzy reliability of 
different systems having components follow different types of conflicting bifuzzy 
failure rates. We also fabricate the membership degree and nonmembership degree 
of fuzzy reliability Kumar and Yadav [9]. Further, using introduced method, mem- 
bership degree and nonmembership degree of fuzzy reliability of different systems 
like series system, parallel system, parallel-series system, and series- parallel sys- 
tem. The components of above systems follow different types of conflicting bifuzzy 
failure rates. We have seen that, the fuzzy reliabilities of above systems are in the 
form of trapezoidal conflicting bifuzzy numbers. The a-Cut and -Cut of fuzzy reli- 
ability are also evaluated at different times (t). 
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